Tome 22 | 


N° 7 _Juicier 1961 


LE JOURNAL DE PHYSIQUE 


ET 


LE RADIUM 


LE GROUPE D'’'INVARIANCE ASSOCIÉ AUX ROTATEURS RELATIVISTES 
ET LA THÉORIE BILOCALE 


Par F. HALBWACHS, 


Institut Henri-Poincaré, Paris. 


J. M. SOURIAU, 


Faculté des Sciences, Marseille. 


J. P. VIGIER, 


Institut Henri-Poincaré, Paris. 


Résumé. — Trois résultats sont obtenus dans cet article : 

19 on montre; par une démonstration très générale, que l’ensemble des constantes fondamentales 
du mouvement qui peuvent être tirées de chaque groupe d’invariance d’un Lagrangien donné sur 
une ligne d’univers, fournit des crochets de Poisson classiques qui satisfont à des relations iden- 
tiques aux relations de commutation de l’algèbre de Lie du groupe ; 

29 toute particule comportant des variables internes peut être représentée par un rotateur 
« bilocal » formé de deux tétrapodes relativistes localisés en deux points différents ; 

30 les transformations de jauge à un et trois paramètres du formalisme spinoriel habituel peuvent 
être traduites dans le formalisme des tétrapodes, avec une interprétation cinématique évidente. 


Abstract. — Three noteworthy results are established in this paper : First,it is shown in the 
frame of a very general treatment, that the set of fundamental constants of the motion which can 
be derived from each invariance group of a given line-Lagrangian gives rise to classical Poisson 
Brackets which fulfil relations identical to the Lie algebraic relations of the group. Secondly, 
any particle endowed with internal variables can be represented by a ‘* bilocal ? rotator composed 
-of two relativistic tetrads localized at two different points. Finally, the usual one-parameter and 
two-parameter gauge transformations in spinor form can be translated into tetrad form with an 
obvious kinematical interpretation. 


Introduetion. — On se propose d'étudier du 
point de vue de la théorie des groupes, la dyna- 
mique de la particule relativiste en mouvement 
dans lunivers physique. Par univers physique, 
_ nous entendons l’espace-temps de Minkowski pos- 
_sédant une orientation déterminée — plus précisé- 
ment, un sens du temps relié à l’existence d’un cône 
de lumière avec une nappe future déterminée, et 
une orientation déterminée de l’espace telle qu’elle 
peut être mise en évidence par les phénomènes qui 

ne conservent pas la parité. En outre, cet univers 
peut être le siège d’un champ extérieur capable 
d'agir sur la particule en chaque point. D'autre 
part, la particule sera, au moins au début, consi- 
dérée comme strictement localisée en un point 
d’univers x, qui décrit une ligne d’univers £. Outre 
sa position, elle comportera des « variables 1n- 


ternes » b® variant avec le temps et sur la variance 
et la nature desquelles nous ne ferons d’abord 
aucune hypothèse restrictive. Toutes les propriétés 
dynamiques de la particule seront déduite d’un 
Lagrangien dépendant de ,, &,, b®, b® (nous dési- 
gnons par un point la dérivation par rapport au 
temps propre + de la particule, c’est-à-dire le long 
de la ligne £). Pour traiter ce problème, démon- 
trons d’abord de façon purement mathématique 
une propriété très générale. 
I. Le théorème de Noether et les variables dyna- 

miques. — Considérons les variables suivantes : 

T (à valeurs réelles) 

Q = (x) (à valeurs dans A?) (1) 


L = jf(+,Q, O) (à valeurs réelles) 
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et l’intégrale : 


T 
4= [Las 


Nous pouvons évidemment 
écrire : 


poser y = Q et 


4= [ue y + PQ — y à 
44 = (PQ + [He Qu — PO— Par. (2 


P étant une fonction quelconque de + à valeurs 
dans le dual de A" (nous supposons que Q et y 
sont des colonnes et que P est une ligne. Leurs 
éléments seront respectivement désignés par Q?, 
y”, P3). Si à représente une variation de la fonc- 
tion ®, nous aurons : 


D4, = [PS0 


de [o— b).50 + CE P) Sy] dr (3) 


Nous pouvons profiter de l'indétermination de P 
pour poser 


_bL 
i 


afin de nous débarasser du coefficient de dy dans (3). 
Nous obtenons alors les équations de Lagrange : 


dL dL . 
6 0 
en tant que conditions nécessaires et suflisantes 
pour que : 


P 


y= Ÿ (4) 


[S Dr) 00) — 0] = [> 4i ==NE (5) 


La forme même de (4) montre que les condi- 
tons (5) ne dépendent pas du choix de +, et +1. 

Considérons, maintenant, une variation quel- 
conque ÔQ fonction de + et de Q, mais qui ne 
s’annule pas nécessairement pour les valeurs 7, 
et +,. On a alors : 


SA — [oz dr. 
To 


S1 la condition (4) est remplie, (2) montre qu’on 
obtiendra : 


AË— [P,801: 


Cette relation nous conduit au théorème de 
Noether [1]: À toute variation 3Q fonction de Q 
et de +, les équations de Lagrange permettent 
d'associer une fonction w de P,Q,r : u — P.5Q 
telle que la différence de ses valeurs en +, et +, 
est égale à la variation de l'intégrale # prise de +, 
AT 


SA = U (0) ee Utruie (6) 


Ce théorème peut être appliqué en particulier au 
cas où l’on considère un groupe de Lie G à p para- 
mètres agissant dans l’espace (Q, +) et conser- 
vant r. Les générateurs infinitésimaux Ô, de ce 


groupe fournissent alors p fonctions fondamentales: 


| UT = PORTO | . (7) 


Appelons phase le couple de variables (o) et 


variable dynamique toute variable qui est fonction 
de la phase et de +. Nous pouvons associer à chaque 


variable dynamique scalaire w une différentielle d, 


définie par : 
du du 
du! P;) eu 30 ‘ 


De plus, si u et v sont deux variables dyna- 
miques, l'expression : 


du(r) = 0 dy'Q?) — 


du dv du dv 
 DP,0QX DQX0P; 


sera le crochet de Poisson de u et v. Les différen- 


[u, v] = dy'v) 


tielles associées à ces variables dynamiques satis- 


font à la relation : 
du dy == dy Dr dur]: (8) 


En appliquant cette relation à une troisième va- KT 


riable dynamique w on obtient l’identité de Jacob : 
Lu, [v, w]] + [v, [w, u]] + [w, [u, v]] = 0. 


Appliquons ces considérations aux variables dy- 


namiques définies par (7). Nous obtenons immédia- 
tement : 


duj(Q) DA ÔJ Q 
dy?) = — PE Gi 0) (9) 
dui(r) = 0. 


D'où nous tirons : 
[uÿ, Ux] =, duÿ(P .0% () Fe dui P.Ôdk Q + -P .duj(dr Q) 
= LPS, 01e O0 + P M0 0 
dQ° e 097 | 
Lu, ux] — P.(d; dx — dx d;) Q. (10) 
Mais les générateurs infinitésimaux 5, d’un 
groupe de Lie satisfont aux relations de commu- 
tation 
D5 dx — 06 D; — The dy 
où les coellicients de structure 7}, sont des cons- 
tantes qui restent invariantes dans tout homo- 


morphisme du groupe. De la relation (10), nous 
trons donc l’identité : 


lé url T} Uy (11) 


N°07 


N° 7 


qui relie au moyen des coeflicients de structure les 
fonctions fondamentales engendrées par les géné- 
rateurs du groupe G. 

La relation (11) peut être considérée comme la 
forme infinitésimale d’une propriété finie très 
simple. En effet, LA ne une transformation 
bi-différentiable o = — F(Q), nous pouvons lui asso- 
cler une transformation corrélative des variables P 
correspondantes en posant : 


P'.4Q' = P.4Q (12) 


ce qui donne : 
= P.4QfdQ’. 


De même, cette transformation qui laisse inva- 
riante la forme P.3Q. laisse évidemment aussi inva- 
riante la forme différentielle dP.5Q — SP.dQ. 
C’est donc une transformation canonique dans l’es- 
pace de phase. 

Si maintenant, nous considérons un groupe G de 
transformations agissant sur @, nous obtenons, en 


appliquant (12) à chaque transformation de C: un 


isomorphisme entre G et un groupe F de transfor- 
mations canoniques. Il est facile de voir que, dans 
cet isomorphisme, chaque transformation infini- 
tésimale à; du groupe G, correspondra à une trans- 
formation canonique infinitésimale du groupe l qui 
coïncide avec la différentielle d,; définie par (9). 
L’invariance des relations de commutation par un 
isomorphisme nous donne alors la relation : 


dui dur — dur duÿ = Tir du. (13) 


Mais, d’après (8), on peut écrire le premier 
membre de (13): du Ce qui donne liden- 
tité (11). 


Signalons que les résultats démontrés ici peuvent 
facilement être étendus à des cas plus généraux, par 


- exemple à des transformations agissant aussi sur le 


_ paramètre + (on exprimera + et Q en fonction d’un 


paramètre arbitraire s)-ou au cas de Laägrangiens 
dépendant de dérivées supérieures (pour le cas des 
dérivées secondes, cf. plus loin, section 4). 

Nous allons appliquer maintenant les résultats 
généraux ci-dessus au problème de la particule. A 
chacune des variables définies plus haut : 

_ À 

Q = (Lu dt )) 
correspondra son moment canoniquement con- 
jugué : 
P=! Guy = LE pa dE” 

dy 2b® 
et à chaque générateur infinitésimal à, d’un groupe 
de Lie G, Hate une fonction : 


Le groupe G ae être formé de deux sortes de 
transformations : les glissements qui agissent sur 
les +,, et les transformations de jauge qui agissent 
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seulement sur les variables internes b®, si bien que 
le groupe G se présente généralement comme un 
produit direct d’un groupe de glissements et d’un 
groupe de jauge. 

En ce qui concerne les glissements, nous pouvons 
considérer que, si le Lagrangien dépend des coor- 
données x,, cette dépendance exprime l’influence du 
champ extérieur A(x), et par conséquent nous pou- 
ons supposer que L ne dépend de x,, que par l’in- 
termédiare de 4,;, qui est une fonction vectorielle 
de æ,. Ainsi tout générateur infinitésimal à, du 
groupe des glissements produit un déplacement 
infinitésimal de la ligne £ dans un champ extérieur 
immobile, si bien que la variation correspondante 
de l’intégrale 4 est : 


NS See dL dus, 


D on 


-et que l’équation (6) devient : 


LS pa dL dAus. 
y To d Au dTy 


soit, sous forme différentielle : 


E dL dAy 
HR dAy 0% CE 


ui(T1) — Wi(To) 4 Tv AT 


Ou enfin, plus simplement : 


(14) 


d; A, représentant la variation subie par le champ 
extérieur À, sous un glissement infinitésimal à, 
dans l’espace-temps. 

Nous voyons, en particulier que les transfor- 
mations infinitésimales à, qui laissent invariante 
l'intégrale d'action #5 sont données par la relation : 

dL 


Il peut exister pour un champ de structure donné 
certaines formes de Lagrangien qui définissent un 
sous-groupe à, de glissements pour lesquels les vec- 
teurs dLfdAy et à, A, sont orthogonaux. C’est le 
cas, par exemple, pour les rotations spatiales dans 
un champ à symétrie (spatiale) sphérique, ou pour 
les translations dans un plan orthogonal à un 
champ uniforme. 

L'étude des glissements infinitésimaux dans un 
champ extérieur fournit les lois de la dynamique 
d’une particule relativiste en rotation soumise à ce 
champ. L’un de nous [2] a utilisé cette méthode 
pour traiter le problème de la précession du spin 
d’une particule chargée en mouvement dans un 
champ électromagnétique. Néanmoins le cas de la 
particule libre, en l’absence de champ extérieur, 
qui parait limité à un « univers physique » très 
particulier, est en fait très important, car c’est 
toujours dans ce cas qu’on se place pour construire 


Ô; Au 10! 


8 
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des modèles du niveau subquantique de la matière, 
niveau où le concept de champ continu perd toute 
signification et doit laisser la place à la notion 
d'interaction, deux particules subquantiques en 
interaction devant être considérées comme for- 
mant un édifice unique pendant la durée de l’inter- 
action. Dans toute la suite, nous allons nous placer 
dans le cas de la particule libre en supposant 
A, = 0. Alors, le Lagrangien ne dépend plus expli- 
citement de x, et reste ineariant sous un glissement 
quelconque. De plus, en ce qui concerne les transfor- 
mations de jauge, nous nous limiterons également à 
celles qui laissent le Lagrangien invariant, si bien 
qu’à partr de maintenant, nous allons seulement 
considérer le groupe d’invariance du Lagrangien 
étudié. 

Dans ces conditions, l’équation (14) devient : 
UN 0 
c’est-à-dire : 

mu = Pi Q= Guy dir + EN 5; bN = constante. 


Chaque générateur infinitésimal du groupe d’in- 
variance donne ainsi naissance à une constante du 
mouvement, et ces constantes fondamentales expri- 
mées en fonction des variables canoniques, obéis- 
sent aux relations algébriques (11) sur les crochets 
de Poisson. 


II. Le groupe des glissements et la théorie bilo- 
cale. — Les Lagrangiens qu’on peut avoir à consi- 
dérer admettent différents types de groupes d’in- 
variance. Mais, dans le cas de la particule libre, le 
groupe d’invariance comporte nécessairement 
(éventuellement comme sous-groupe) le groupe gé- 
néral des glissements dans l’espace-temps, c’est-à- 
dire le groupe inhomogène de Lorentz. Dans cette 
section, nous appliquerons le théorème précédent 
au groupe de Lorentz, sans spécifier les variables 
internes. 

Comme on le sait, les générateurs infinitésimaux 
du groupe inhomogène de Lorentz sont : 

19 les opérateurs 9, des translations infinité- 
simales dans l’espace-temps. Ceux-c1 agissent uni- 
quement sur les x, et sont donnés par : 


Tu Th = Eu 


20 les opérateurs AC,, des rotations infinitési- 
males dans l’espace-temps, qui agissent à la fois 
sur les x, et sur les variables internes. Sur les 
premiers, leur action est exprimée par : 


JA uv Th = Lp Eva — y Bu. 


Sur les variables internes, la forme explicite des 
opérateurs JAC,, dépend de la variance de ces 
variables. 

- Pour former les constantes fondamentales du 
mouvement correspondant à ces opérateurs, nous 
tirerons d’abord du Lagrangien Z les moments 


JOURNAL DE PHYSIQUE 


Ne 


canoniquement conjugués aux variables de confi- 
guration, soit : 


Dans ces conditions, l’opérateur %, engendre le 
vecteur constant : 


GNT 2 = ON En Gn 


c’est-à-dire simplement le vecteur impulsion. D'autre 
part, l’opérateur AC,, engendre le tenseur anti- 
symétrique constant : 


Muv == GX ACuy TA + Ba AC pb) = Ty Gy + AY Gy + Suv. 


C’est le moment angulaire total qui se compose de 
deux termes dont la signification physique est évi- 
dente : æ, G, — x, G, est le moment angulaire 
orbital relatif à l’origine. Le tenseur 

AE 55 pa My pa) 
dont la forme dépend de la variance des b® est le ” 
moment angulaire propre qui caractérise la rotation à 
interne de la particule, et dont le rôle canonique a 
déjà été mis en évidence par deux d’entre nous [3]. 

Pour nous résumer, l’invariance du Lagrangien 
par rapport au groupe des glissements nous fournit 
deux relations fondamentales de conservation : 


Gy = 0 (15) 


Ms = 0. 


Remarquons que si on explicite ,, dans la 


seconde équation, celle-ci devient : 


SES = Gy Ty ES G Zu (16) 


en tenant compte de (15). Ainsi nous aboutissons au 
système (15), (16) où l’on reconnaît le système de 
Weyssenhoff [4]. : 

D'autre part, on sait que l’algèbre de Lie du 
groupe inhomogène de Lorentz est construite sur les 
relations de commutation : 


[Tu sh] —= 0 


j [ACuv, Ta] SERA Ty — Lav Tu 


[Gus Ag] = Sue Ag + V8 Aux -— gup Ave — va Aug. 


A ces relations, vont correspondre des relations 
classiques exprimant les crochets de Poisson de nos 
constantes fondamentales, conformément au théo- 
rème général établi ci-dessus, soit : 


Re ( 
LMuv, Gal = gou Gv — gav Gu 


Mouv, Mag] = Eux Mog + 5,8 Mux 


Lu Mu — Eva Mug: 


(17) 


Remarquons que le tenseur M, n’est défini que 


Po Yu = —c 


No.7 


:ÉLael à une origine arbitraire. Mais, comme 

’écrit M. L. de Broglie [5] : « Dans chaque système 
Fe référence, l’origine des coordonnées est un point 
arbitraire, et le moment cinétique par rapport à ce 
point arbitraire n’a pas, en général, de signification 
physique particulière. Ce qui a une signification 
physique intéressante, c’est le moment cinétique 
par rapport à un centre doué de propriétés phy- 
siques. » Nous allons montrer, à partir des relations 
de conservation ci-dessus qu'un tel centre physique 
apparait très naturellement en liaison avec l’étude 
générale du mouvement. Ce centre est lié à la parti- 
eule et, par conséquent, est en moüvement, mais 
on verra que nous pouvons cependant le prendre 


comme origine pour évaluer le moment angulaire 


orbital. 
Nous pouvons écrire la grandeur constante de 
impulsion sous la forme : G, G, = — Mic, M, 


ayant les dimensions d’une masse. Si nous consi- 
dérons alors le moment angulaire propre S,, et si 


Rey x D nv-Cll 


aura les dimensions d’un rayon-vecteur dans 
Pespace-temps. Appliquons alors l’origine du vec- 
teur — À, au centre x, de la | particule. Son extré- 
mité définira un point ; — R,. Ce point 
dont le rôle a été pressent: par MO [4] et 
qui a été mis en évidence par l’un de nous [6], nous 
pouvons l’appeler le «centre de gravité », il possède 
la particularité de se mouvoir d’un mouvement uni- 
forme parallèlement au vecteur impulsion, comme 
c’est le cas pour le centre de gravité classique d’un 
système mécanique non relativiste en l’absence de 
champ extérieur. En effet, dérivons Ÿ, par rapport 
au temps propre + relatif au point x, : 


nous formons le vecteur 


. 2 ° 4 ° 
Yu = Ly Ry = Ty M ze Suv Gy 
1 : ve! 
= Lu Mic (Gu vi Cy du Gy 
ou encore, en posant G,d, — — A, €? (la quan- 


tité AC,, généralement variable, a encore les dimen- 
sions d’une masse), nous obtenons : 


AC 


Taie M? 


2 Gp. 


Comme Ge — 0, le mouvement est rectiligne. Sa 
Vitesse LETTRES constante “ est obtenue en faiient 
intervenir le temps propre d'relatif au point Y,: 


— ÿy Yuldr/d0)? 


— (AIME) (ME ce?) (dx [d0)° 
c’est-à-dire : 
drfd0 = MIA 
et 
Ye d+/d0 — GylM, = constante, 


74 
Yu = 
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Maintenant, nous pouvons exprimer le moment 
angulaire total relativement à ce centre de gravité, 
soit : 


AC = By (er EE Ry Guy +- iSyLv. 


Ce moment aura une signification physique in- 
trinsèque indépendamment du choix de l’origine 
des axes, si bien qu’on aura avantage à l’introduire 
à la place de M,,. D'autre part, il sera encore une 
constante du mouvement et ainsi vérifiera encore 
les relations de Lie données ci-dessus. En fait, i! 
n’est rien d’autre que le moment angulaire relati- 
vement à une origine fixe, mais choisie sur la droite 
d’univers déterminée parcourue par le centre de 


gravité, si bien qu'à un instant t le vecteur Y, 


eompté à partir de cette origine, est colinéaire 
à G,. I s'ensuit qu'on peut ajouter au moment 
angulaire total M, — x, G, — x, G, + S,,, le 
tenseur Ÿ, G,-—.Y, G, qui est nul. On obtient 
ainsi : 


Ms = Tu RS TN GENE SET ANS 


Ainsi à toute particule régie par un Lagrangien 
invariant par glissement, et comportant des va- 
riables internes localisées en un point déterminé x, 
nous sommes conduits à associer un autre point Ÿ, 
en mouvement uniforme, autour duquel le point x, 
exécute un mouvement plus où moins compliqué. 
Nous retrouvons les deux points « canoniquement 
conjugués »récemment mis en évidence par T. Taka- 
bayasi [7]. De plus, le mouvement de la particule 
nous fournit au point Y, un système déterminé de 
vecteurs que nous pouvons valablement choisir 
comme référentiel relativement fixe par rapport 
auquel on étudiera le mouvement du point x, et 
des « vecteurs internes ». En effet, considérons le 
vecteur G, et le tenseur AC,, comme localisés au 
point Y,. Nous aurons ainsi un premier axe coli- 
néaire à G, et qui va glisser sur lui-même d’un 
mouvement uniforme. Puis nous obtiendrons un 
second axe orthogonal au précédent en considérant 
la « polarisation » [8] (aussi appelée « spin de 
Moller [6]»), à savoir : 


L 

Au = 2M,c Euvap Gv A af 

(Euves est le tenseur complètement antisymétrique 
de Lévi-Civita). Ce vecteur, évidemment constant, 
va se déplacer parallèlement à lui-même. Finale- 
ment, nous choisirons, dans l’hyperplan orthogonal 
à Gy, et À,, qui reste parallèle à lui-même, deux 
autres axes perpendiculaires, qu'il est possi- 
ble — et tout indiqué -— de prendre encore de 
directions fixes, et qui seront alors définis dans cet 
hyperplan à une rotation constante arbitraire près 
(qui jouera le rôle de jauge). Ainsi nous définissons 
un «tétrapode » orthonormé formé de quatre vec- 
teurs unitaires 4 te est un indice non tensoriel 
prenant les valeurs 1, 2,3, 4) dont le quatrième 46 


Pis 
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est du genre temps. Ces vecteurs satisfont aux deux 
conditions équivalentes d’orthonormalité : 
AD 40 2 ps AS AUD 2 GE 

Nous appellerons ce tétrapode le système fixe. 

A cette étape, il apparaît que notre formalisme 
très général, qui paraissait initialement décrire une 
particule ponctuelle, engendre nécessairement un 
modèle « bilocal ». Il se rattache ainsi à de très 
nombreuses théories proposées récemment, qui con- 
sidèrent d'emblée, soit deux points distincts [9], 
[10], soit tout un édifice, par exemple une goutte- 
lette fluide [6], [11]. En particulier d’après cette 
dernière conception, le fluide est caractérisé loca- 
lement par un vecteur conservatif « densité de 
matière » 7,, et par un tenseur symétrique conser- 
vatif «densité d’énergie-impulsion » {,,. L’impulsion 
globale G, est obtenue en intégrant sur une section 
quelconque d’espace Z la composante t{,, du ten- 
seur {,, relative au référentiel correspondant à 2 : 


1 
Gu = [ru dv. 


C’est un vecteur constant. Le rôle du point x, 
considéré ici est alors joué par le barycentre X, 
relatif à la densité de matière 7, dans une section 
d'espace particulière 2, orthogonale à G, ; c’est le 
«centre de matière » : 


A. [ja du = fe ün Ja AD. 


Le rôle du point Y, est joué par le barycentre 
relatif à la densité d’énergie t,, dans la même sec- 
tion d’espace ZX, ; c’est le «centre de gravité » : 


Yu [ aa AU = L Ty taa A. 
0 0 


On montre alors facilement que le moment angu- 
laire relatif au point Y, calculé par intégration sur 
une section d’espace quelconque, soit : 


1 
ACuy = = U Lau — Vu) tva — (xs — Yy) tua] do 


1e 

est un tenseur constant, et que le moment angulaire 
relatif au centre de matière, 

1 

S uv ec 2 
1C JE 
est aussi un tenseur satisfaisant à la relation : 

S uv —= Gy x, ee Gy Xp. 

Ainsi tout le formalisme de la gouttelette est 
entièrement identique à celui que nous présentons 
ici, On arriverait à la même conclusion pour la 
théorie développée il y a longtemps par Matthisson 
et ses élèves [12], [13] qui étudient à différents 
degrés d’approximation l’évolution d’une particule 
considérée comme une singularité du champ de 
gravitation régi par les équations d’Einstein. 

Il nous paraît digne de remarque que le forma- 


Lau — Zu) tva — (tv — À) tu] dv 


lisme très général proposé ici, qui à l’origine ne 
semblait pas attribuer à la particule une extension 


dans l’espace, nous a conduit à un modèle bilocal 


qui concorde exactement avec les résultats obtenus 
par la démarche opposée, à savoir en schématisant 
des structures étendues dans l’espace. Cette re- 
marque peut être considérée comme un argument 
en faveur de l’opinion générale suivant laquelle le 
fait d'attribuer à une particule soi-disant ponctuelle 
d’autres propriétés qu’une pure et simpleimpulsion 
— par exemple, de la doter d’un «spin » — revient 
en réalité à lui reconnaître implicitement une cer- 
taine extension dans l’espace. Selon cette concep- 


tion, on serait toujours, en fait, en présence d’une. 
structure étendue dans l’espace, schématisée par. 


différents modèles correspondant à des approxi- 
mations différentes. Ainsi la première approxima- 
tion représenterait l’ensemble de la particule par 
un point doté d’une impulsion, celle-ci étant néces- 
sairement colinéaire à la vitesse. Le formalisme 
général développé ci-dessus montre alors immédia- 
tement que le moment angulaire total se réduit à sa 
partie orbitale sans moment angulaire propre. Le 
rayon-vecteur À, s’annule, si bien que les points x, 
et Y, coïncident ; le modèle devient unilocal et le 
choix de l’origine sur la droite d’univers décrite par 
la particule annule le moment angulaire total A0, 
Le point comportant en outre des variables internes 


tel que nous l’avons étudié ci-dessus, représenterait . 


alors une seconde approximation qui se traduit par 
le modèle bilocal. 


I semble aussi que cette remarque puisse être le. 
point de départ d’un examen critique de la concep- 


tion de la particule à spin telle qu’elle est utilisée 
par la mécanique quantique : cette particule est 
considérée comme rigoureusement ponctuelle, ce 
qui se traduit par l’apparition de la fonction 3 de 


Dirac dans le formalisme, et en même temps elle est - 


dotée de spin et d’autres propriétés spécifiques — 


par exemple, de moment magnéto électrique — qui . 


du point de vue classique sont en fait des « va- 
riables internes ». Mais, d’autre part, une telle par- 
ticule quantique est à la fois strictement ponctuelle 
et sans localisation précise, étant donné les rela- 
tions d'incertitude. Peut-être est-il permis, au con- 
traire, de rechercher un modèle « déterministe » qui 
pourrait être à la fois strictement localisé et non 
ponctuel, ou, pour parler plus précisément « multi- 
ponctuel ». 


IIT. Le tétrapode fondamental. — Après ces 


considérations générales, nous allons préciser notre : 


modèle pour retrouver les représentations phy- 


siques ordinaires de la particule tournante rela- 


tiviste. Ainsi nous devons introduire nos variables 
internes de telle façon qu’elles puissent décrire 
convenablement une rotation relativiste, c’est pour- 
quoi nous allons généraliser de façon naturelle la 
représentation de la toupie non-relativiste. 


Ke 
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On sait que celle-ci est représentée par un trièdre 
mobile de vecteurs orthonormés bf”, qui nous per- 
met de définir un tenseur antisymétrique «rotation 
instantanée » &;; — b® b® (le point représente la 
dérivée par rapport au temps et les indices r, à, j 
prennent les valeurs 1, 2, 3). Outre son caractère 
orthonormé, le trièdre doit être orienté par le choix 
d'un certain sens de rotation. Si on se donne deux 
des vectèurs et l'orientation, le troisième vecteur 
est complètement déterminé par la relation : 


(r) rst ,(s l 
d— se "bi EE 


(x et <"* sont les symboles complètement anti- 
symétriques du troisième ordre de Lévi-Civita). 
Cette relation entraîne l’orthonormalité qu’on peut 
exprimer par les relations équivalentes : 
HU = guy D D = 

De plus, chaque vecteur est déterminé sans ambi- 
guité à partir des deux autres, dès lors qu’on a 
choisi des permutations positives pour 2.7.k et 
pour r.s.t, soit, par exemple &.7.k = r.s.t — 1, 2, 
3. Remarquons que £;;4 est un pseudo-tenseur com- 
plètement antisymétrique, tandis que les e"* sont 
des quantités égales à + 1 ou — 1 et évidemment 
invariantes sous une rotation « propre », c’est-à-dire 
ne comportant pas d’inversion d’espace. Mais nous 
avons besoin d’une convention supplémentaire si 
nous désirons attribuer une variance aux e'# sous 
une inversion d’espace, et cette convention déter- 
minera en même temps la variance des vecteurs 
formant le trièdre. Si les e’* sont invariants sous 
une inversion d'espace, alors, comme le produit 
b® b® est évidemment invariant, tandis que le 
pseudotenseur £;;, change de signe, les compo- 
santes b® sur les axes inversés deviennent — b®,. 
En d’autres termes, les b® sont de « vrais » vecteurs 
(vecteurs polaires). Au contraire, si nous supposons 
que les £"* sont des pseudoscalaires qui changent 
de signe sous une inversion d'espace, alors les b® 
garderont le même signe sur les axes inversés et le 
trièdre ainsi défini sera formé de pseudovecteurs 
(vecteurs axiaux). 

La généralisation relativiste nous conduit de 
même à considérer un trièdre de quadrivecteurs 
contenus dans l’espace propre de la particule et à 
le compléter par un quatrième quadrivecteur du 
genre temps colinéaire à sa vitesse unitaire x, (doré- 
navant le point représente la dérivation par rapport 
au temps propre + de la particule). Pour exprimer 
l’orthonormalité et l’orientation d’un tel « tétra- 
pode », nous devons écrire deux conditions : 


pÉPE 2 o7e (18) 
ic 


‘ C rst j(r) (8) (1) / 
Lu = TE Euvap € by bee bp (19) 


De ces relations résultent, d’une part, les deux 
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relations équivalentes d’orthonormalité relativiste : 
DE D = guv De D = g°1 (20) 


(cette fois, les indices grecs &, n désignant les vec- 
teurs, prennent, comme les indices tensoriels y, v, 
les valeurs 1, 2, 3, 4, bŸ représentant à, Jic, alors 
que les indices latins r, s. t, prenant les valeurs 1, 2, 
3, désignent spécialement les trois vecteurs du 
genre espace) ; d’autre part, une orientation déter- 
minée du trièdre du genre espace et du vecteur du 
genre temps, ces orientations étant fixées sans 
ambiguïté par le choix d’une permutation positive 
des tvaf (par exemple p, v, «, B = 1, 2, 3, 4) et 
des r.s.t ; (par exemple, r, s, t — 1, 2, 3). Comme 
dans le cas précédent, les e'* sont invariants sous 


le groupe connexe de Lorentz (transformations 


orthochrones propres). Mais nous devons fixer con- 
ventionnellement leur variance sous un renver- 
sement du temps (transformation antichrone 7) 
et sous une inversion de l’espace (transformation 
impropre ou parité P). Il y a ainsi quatre conven- 
tions possibles qui définiront quatre types diffé- 
rents de variance pour le tétrapode et qui cor- 
respondant naturellement à quatre cas physiques 
différents : 

10 e*# est un « vrai » scalaire, il est invariant 
aussi bien sous T' que sous P : on a un tétrapode 
«propre » ; 

20 e# change de signe sous 7 et est invariant 
sous P. On a un tétrapode «antichrone » ; 

39 e# change de signe sous P et est invariant 
sous 7. On a un tétrapode «impropre » ; 

40 et change de signe aussi bien sous P que 
sous 7”. On a un tétrapode «invariant de charge » 
(parce qu’il est invariant sous PT qui est la conju- 
gaison de charge). Naturellement, dans chacun de 
ces cas, la forme du Lagrangien doit être choisie de 
façon qu’il soit invariant sous une transformation 
de groupe complet de Lorentz, compte tenu des 
transformations subies par les tétrapodes du type 
considéré. 

Ceci étant dit, et la possibilité des différentes 
variances étant réservée par l'éventualité de con- 
ventions correspondantes sur les Z*%#, disons tout 
de suite que nous laisserons de côté, dans cet 
article, l’étude des transformations par inversion, 
problème très important que nous nous proposons 
d'aborder dans un travail ultérieur. Dans tout ce 
qui suit, nous allons donc restreindre les glisse- 
ments aux transformations orthochrones propres de 
Lorentz. La rotation instantanée relativiste du 
tétrapode s'exprime évidemment par le tenseur 
Ov — L® b® avec sommation sur £. Ce tenseur 
généralise de façon naturelle le tenseur non-relati- 
viste de rotation instantanée : relativement au sys- 
tème propre, nous devons ajouter aux composantes 
d'espace : 3; — 0% b®) les composantes de temps 
ou — 0® bŸ = bŸ qui représentent l'accélération 
conformément aux conceptions relativistes. Aïnsi 


PS 
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l’évolution dans le temps de notre tétrapode fonda- 
mental représente en même temps la rotation et 
l'accélération du point matériel considéré, ces deux 
grandeurs étant unies en un tout indissociable dans 
tout référentiel autre que le système propre. 
Comme nous l’avons dit plus haut, ce modèle 
ponctuel doit être nécessairement complété par un 
centre de gravité auquel est lié un autre tétra- 
pode qui se déplace parallèlement à lui-même avec 
une vitesse constante. En effet, une fois qu’on a 
choisi un Lagrangien déterminant l’évolution 


res: dL 
des b®, nous pouvons en tirer l'impulsion G, = — 
Ty 
et le moment angulaire propre 
(iQ ere b 
Suv = by 2bE } dbE 
d’où le rayon-vecteur 
| 1 
Ru = Fo a Gy 


qui nous fournit le contre de gravité Y, = x, — A. 
De plus, le moment angulaire total (relativement 
au centre de gravité) est 


M y —= Ro Gy + re Gy +- iS av. 


On peut alors définir les deux premiers vecteurs 
du système « fixe » : 


(4 TAC. G MARS 
D) = GulV— A Gr er D de 
avec 


(4 
Au = = Euvop Us Ag. 
Enfin les deux derniers sont déterminés, à une 
jauge constante arbitraire près, par les équations : 


(1) (2) 
du = dy 


1 ne 
sr) COS rst (Tr) ,(s) ag) 


On IE nval 2 RUN C 
! A 010) v Co 


(£) (t) St 


Cu  =£. 


Ainsi nous aboutissons à un édifice formé de deux 
tétrapodes orthonormés localisés en deux points 
différents, liés Pun à l’autre. Cet édifice bilocal, 
décrit en premier lieu par Takabayasi et ses colla- 
borateurs [7], [20], et qu’on conviendra d’appeler 
un rotateur relativiste représente par un schéma 
unique, d’une part, les modèles présentés explici- 
tement comme bilocaux [9], [10], d’autre part, les 
modèles pseudo-ponctuels tels que, par exemple, 
la particule à spin ordinaire, et aussi les modèles de 
particules étendues convenablement  schéma- 
tisés [11]. Mais nous pensons que cette représen- 
tation par deux tétrapodes est la plus suggestive et 
aussi la plus maniable. Cette remarque permet de 
donner un très vaste champ de validité à la théorie 
récemment proposée par deux d’entre nous [14] 
qui décrit la dynamique d’un tel rotateur par le 
mouvement « global » ou « externe » du centre de 
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gravité, et par le mouvement «interne » consistant 
dans l’évolution du vecteur À, et de l'orientation 
relative du tétrapode « mobile » par rapport au 
tétrapode « fixe », cette orientation étant exprimée 
au moyen de la généralisation relativiste des angles 
d’Euler [15], ce qui permet de former des opéra- 
teurs de rotation susceptibles d’être quantifiés [16]. 

Nous pouvons immédiatement appliquer ces 
considérations sur le tétrapode fondamental au cas 
très simple du « Lagrangien linéaire » introduit par 
Takabayasi et deux d’entre nous et qu’on a fré- 
quemment utilisé par la suite [17], [18], [19], [21. 
Un premier terme représentant l’énergie cinétique 
de rotation propre est formé de la façon suivante : 
plaçons-nous dans le cas dit de Frenkel [21] où le 
moment angulaire propre de la particule est un 
tenseur d'espace propre : 


S pv Ly Ce 0. $ 


Alors, si nous formons le spin 
L . 
Su — 2e Euvag Év S ap; 


la relation précédente entraîne la possibilité d’ex- 
primer #,, entièrement en fonction de #, ets, [6] 


; " - 
SN ; Euvap La SR: 


Prenant alors le vecteur b® sur à, et le vecteur 
DESutiS 
RE ict En SH 0 
(s, est le module du spin), nous obtenons : 
Spy = So cuve De ee 


c’est-à-dire, d’après (18) et (19) : 


s (1) ,(2) (1) ,(2 
CE CM TE 


L'énergie cinétique peut alors être formée au 
moyen de l'expression classique : 


1 (1) }(2) (1) 3@)) 3 1) 
= soft Es, —Ù by, hrs 


2 


T — 


DORE 


S'y Ouv — 
(1) ; (2 j (I 2 1) ; (2 
> UE bi SE ba bi) 50 be be 


Les autres termes s’annulent d’après les condi- 
tions d’orthonormalité (20). En complétant ce | 
terme par les conditions de Lagrange discutées dans 
cette section, le Lagrangien devient : 


Si HMS ER RATE 8 
Set re Me 7 LA Re 
jap ie | mst Cr) j(s) z(4) 
+ Au (du EYE EuvaB = Ly la te 


PE : : | 


, et N° étant des multiplicateurs de Lagrange 
(A est évidemment symétrique en r et s). 
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Un tel Lagrangien est invariant sous un glis- 
sement général (transformation inhomogène de 
Lorentz). De plus, il est également invariant sous 
une rotation de b® et b® autour du plan 09, b®, 
si bien que le groupe d’invariance est le produit 
direct du groupe inhomogène de Lorentz et du 
groupe des rotations planes orthogonales à b®, b@, 

Les générateurs infinitésimaux seront : 

— l'opérateur f, des translations infinitésimales 
agissant Sur à, ; 

:— lopérateur A0, des rotations infinitésimales 
dans l’espace-temps, agissant sur x, et sur les b® ; 

— l'opérateur # des rotations planes infinité- 
simales orthogonales à b®, b® agissant sur D et 
:L 

Pour exprimer les constantes du mouvement 
correspondantes, constituons d’abord le formalisme 
canonique. Les seules variables dont les dérivées 
figurent dans le Lagrangien sont x, et b®. Il faut 
donc introduire deux moments conjugués : 


Ca = DEP = FL = LE = 50 bp. 


Remarquons que, comme le moment conjugué 
de b® n’est autre que b® il est possible de former 
un Hamiltonien, quoique le Lagrangien soit li- 
néaire en b®. On obtient : 


He: a CV pis) 7, . 6x rst pl” 15) bp. 


Euvap € 


Ceci étant une expression abrégée dans la- 
quelle b® doit être remplacé chaque fois qu’il appa- 
rait par b®/s,. Un cas analogue a été étudié par 
l’un de nous [19]. 

Le vecteur constant correspondant à l’opéra- 
teur 9, est l'impulsion G,, comme on la vu plus 
haut. L'opérateur A, agit sur un vecteur quel- 
conque {Ah Suivant : 


ACuv 9x — Qu En — D Eux. 


Le tenseur constant correspondant est donc : 


* (1) (1) 
Mys = Gi AGuy TX + Ex Cu UN 


(2) ,(4) 
QUES (an RE 


2 1 
— Gy ru — Gu 2 +Ss a . )), 


Le dernier terme est donc simplement le moment 
angulaire propre qui avait été utilisé initialement 
pour former le Lagrangien. Ainsi M,, est, comme 
il se doit, le moment angulaire total utilisé dans la 
représentation physique, soit la somme des mouve- 
ments angulaires orbital et propre. 

Enfin, l'opérateur #% agit sur bY d’après 
= 27 00, Zus étant la matrice des 
rotations in hnitésumales dans le plan 4,90 


0 1 0 
D == — Il 0 0 
() 0 0, 
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Le scalaire constant qui lui correspond est : 


pou ADD — _ Gt (2) so 0 e2 


— Sp: 
C’est simplement le module du spin. 
Les relations de conservation sont donc : 
Gy = 0 Mas 0 Sp 0 
Nous pouvons aussi écrire les relations de Lie en 
crochets de Poisson : 
LG, G] = 0 
[M uv, Cal =" Bou Gy — Gay Gu 
LMos, Mg] — gux Mig + gug Mio — 808 Mis — 854 Mug. 
Comme l’opérateur Æ commute avec T, et A, 
les crochets de Poisson de s, avec G, et M,, sont 
nuls. Toutes ces relations de commutation se véri- 


fient immédiatement à partir des expressions cano- 
niques : 


Cr = Gr My = tn Co an Gue CN = ECIREIUE 
se = Fa 
IV. Le rotateur quadratique. — Nous allons 


maintenant appliquer notre théorie à un autre 
Lagrangien, proposé initialement par Nakano [22] 
et dont la forme quadratique permet des dévelop- 
pements plus intéressants [23], [14]. 

L'expression de la rotation instantanée 


du be be 
sert encore de point de départ, et on généralise 
expression classique de l’énergie cinétique de rota- 


tion physique élémentaire T =; lu? (la cons- 


tante Z est le moment d'inertie). La généralisation 
relativiste fournit : 


1 
ARE % Tale] Oap Ouy, 


le tenseur d'inertie constant Zsemw étant symé- 
trique relativement aux couples «6 et uv et anti- 
symétrique à la fois en «, 8 et en x, v. Une restric- 
tion intéressante consiste à supposer que la parti- 
cule possède une structure de symétrie hyper- 
sphérique dans l’espace-temps, ce qui donne au 
tenseur d'inertie la forme 


T(gou LBv — Gav ge) 


il 
Luv nn: 


avec une seule constante Z. On obtient done : 


1 L - Th ARTE) NA AU 


4 - ñ Lou Opv ne 


ou, en tenant compte des relations d’orthonor- 
malité : 


1 
r 11e 5e 


PE 
fs 
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(ces conditions seront naturellement introduites par 
la suite par des conditions de Lagrange conve- 
nables). 

Une particularité de ce Lagrangien est la pré- 
sence d’un terme en 0%, c’est-à-dire contenant à. 
Nous avons ainsi à faire appel à un formalisme 
canonique élargi au cas où il apparaît des dérivées 
secondes, formalisme proposé par Ostrogradski [24] 
[25]. La démonstration est analogue à celle donnée 
à la section I. Soit : 


= f(x, Q, 0, 0) 


une fonction des variables Q et de leurs premières 
et secondes dérivées. On pose 


=D =:0 


et on étudie l'intégrale d’action : 


re MO, 
P et IT étant des lignes P4, I, à valeurs dans le 
dual de À”. Si 3Q est une variation quelconque de 


la fonction @, et 5Q une variation quelconque de la 
fonction @, la variation correspondante de #4 est : 


343 — [P.8Qn + [301 


ES 


SL >L 
A: es ES ô — .Ôz T e 
> = P ri). Dee Ê 11) ] d 


r,Q, y, 4] + P(Q — y) + IT(Q — z)] d= 


r,Q, y, 2) — PQ — Py — TIQ — Tz] dr 


en tenant compte de l'hypothèse 
T= 0; Sy = 5Q. 
La condition nécessaire et suffisante pour que : 
[Qt = IQ — 0, DQ(x); Qt) = 0] > [D 4 — 0] 


est fournie par les équations de Lagrange généra- 
lisées : 


d’où : 


dL/dQ = P (dL[dy) — P — IT dLfd: = II 


Si ces équations sont vérifiées, toute varia- 
tion 3,0, 80, fonction de + et de O, qui laisse le 
Lagrangien invariant (transformation infinitési- 
male du groupe d’invariance) donnera : 


83 A3 = 0 = [P.8; Ole + [LS Of 
c’est-à-dire, puisque 7, et +, sont arbitraires : 
ui = P.0; Q +- 11.5; Q — constante. 


Aïinsi, nous pouvons étendre à ce cas le théorème 
de Noether en considérant un espace de configu- 
ration construit sur les variables g* = (Q%, Où), 


avec les moments conjugués généralisés 


dL dL 
"No === Po =—— 
œ 200 a dQ% 


a IT el 


Poe = 1P00 ln) 


Pour tout groupe de Lie G laissant le Lagrangien | 
invariant, agissant dans cet espace de configu- » 


ration, sans changer +, chaque transformation infi- 
nitésimale à, de G engendre une constante du mou- 
vement 
Uj = p.04. 
Maintenant, nous pouvons comme ci-dessus, 


considérer la phase E 


miques w, fonctions de la phase et de +. Nous pou- 


vons attacher à chacune d’elles une différentielle 


dE 
du 


dut) = — dulPa) = 


dult) = 0 dPa dq® 


et de même définir le crochet de Poisson de deux 
variables dynamiques 


GE Un du dv du dv 
VASE dPa 0% d4* dPa 
autrement dit : 
DES OU du DV 
[uv] = | = 
d Pc 0Q® DIT, 30 


du 0 du dv 
20% 2Po 30 do) 


De ces formules généralisées, nous pouvons tirer 


sans difficulté le théorème suivant : Si les généra- 
teurs infinitésimaux 0; du groupe d’invariance satis- 
font aux relations de Lie : 


Dj 0h — dx d; = Tir de, 


les crochets de Poisson des constantes fondamen- 
tales du mouvement qui leur correspondent obéis- 
sent à la relation de même forme : 


ture Tr Uy. 

Dans notre Lagrangien quadratique, la seule 
dérivée seconde qui figure est #,. Il faut donc 
considérer d’un côté les vecteurs du genre espace b® 
dL 


avec les moments conjugués usuels 8? = VA) , et 


de Pautre le vecteur du genre temps #9 qu’on 


en prenant les conditions d’ortho- 


normalité sous la forme (18), (19). Il faut alors 
considérer la variable x, (qui ne figure pas expli- 
citement) et son moment conjugué qui prend la 
forme particulière 


écrira b® = à 
ic 


dL, d dL 


dËy. dr dËy 


* 
Gu = 


et les variables dyna- 
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et de plus la variable #, considérée comme 
indépendante, et qui à pour moment conjugué 

dL 
D 

dË 
plémentaire, la constante du mouvement associée 
à chaque transformation infinitésimale à; de groupe 
d'invariance s'écrit : 


. En tenant compte de cette variable sup- 


Uÿ — Gu d i Lu = Nu 5; à, . eo 5; CU 


En particularité les générateurs infinitésimaux 
du groupe inhomogène de Lorentz engendreront 
les constantes : 


TRE d£s de 9L 
ji, dr ou 
x re 
Ms = Tu Gy Sen Guy ne CAT dËy 
AL FRET Te in dL 
Éarr by = (r) Ÿ 23 
dy db, db) 


et le moment angulaire propre prendra la forme 


dL 2 Ô0L uw) dL ) dL 


ER b 
ee en db? 


déjà obtenue par deux d’entre nous [3]. 
Le moment angulaire total par rapport au centre 
de gravité sera : 


Au = Ru G A Ry Gy - uv 


et les crochets de Poisson de G, et A0, exprimés au 
moyen des variables canoniques, ,, G,, à,, N,, 
b®, 8% obéiront aux relations de Lie [17]. 

“Nous appliquerons la théorie prise sous cette 
forme au Lagrangien [26]. 


At ie 


1 : (r) 5 { SES : 
L = — F2 Æ (a bu LT: Ty äy) ie À Re 5) 
4 C 


É ) 
= : 1C st j(r) (5) LE 
+ Du fé, —& Euvas © Re Ce Je 


Les moments canoniques sont : 


: 1 : se (r) 1 ,;6n 
Ny == 2e T&y, Gy = y nr VEN bre = 9 Tbie 
On forme facilement S,, : 
Lee : (£) 1 (Ë 
Syuy = 1(u? pr RE à) = FE pile Toy: 


Ce n’est pas autre chose que l'expression homo- 
logue au moment cinétique classique S — Zo, qui 


: d 
correspond à l’expression 7 — G Lw? pour l’éner- 


gie cinétique que nous avons prise pour point de 
départ. 

Le Lagrangien proposé montre deux particu- 
larités intéressantes : d’une part, le moment angu- 
laire propre n’est pas contenu dans l’espace propre, 
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contrairement à ce qui était posé par hypothèse 
pour le Lagrangien linéaire de la section précé- 
dente. Le produit contracté S,,#, — que nous 
avons appelé le «balourd » [6] — devient : 


TE PE opt = po 


d’où : 
S'y ŒY es Ty. 


Le balourd est colinéaire à l'accélération d'univers. 
Cette relation qui remplace, dans le cas considéré, 
la condition de Frenkel $,, #, — 0 du cas précé- 
dent, forme avec les relations de conservation du 
système de treize équations sur treize variables 
indépendantes, qui déterminent entièrement le 
mouvement (aux conditions initiales près). L’un 
d’entre nous [27] a intégré ce système dans le cas 
de conditions initiales assez particulières. 

D'autre part, L se trouve être invariant par une 
permutation cireulaire des trois vecteurs du genre 
espace bŸ. Autrement dit, le groupe de jauge est, 
cette fois, le groupe des rotations tridimensionnelles. 
On a donc à introduire les générateurs infinité- 
simaux du groupe des rotations agissant sur les 
indices supérieurs, soit : 


JC EU 2 9 PO 


gi (à) 
by 
ce qui nous fournit les constantes du mouvement : 


Nous ee ainsi trois st Pa. 
dans le temps qu’on peut écrire : 


EC ük hp 20) OO 
GUL—= Fos ©” = V5 Lu by — r 9 Te by Le (22) 
avec 6° — 0. On sait que l'opérateur de spin qu’on 


peut écrire S — 5 ik AC obéit aux relations de 
commutation du groupe des rotations tridimen- 
sionnelles : 


[SÈ, 85] — ei Sk, (23) 


On en déduira donc les relations correspondantes 
sur les crochets de Poisson des fonctions o/, soit : 


[oi of] — eûk ot 


qui peuvent être aisément vérifiées à partir de (22). 
Ces relations (23) doivent être jointes aux relations 
de Lie (17) pour définir complètement l'algèbre de 
Lie du groupe d’invariance pris dans son ensemble 
(bien entendu les crochets des 6° avec les G, et 
les A, sont nuls). 

Pour mettre en évidence la signification phy- 
sique de ces trois nouvelles constantes, tirons le 
vecteur relativiste de spin de l’expression du mo- 
ment angulaire propre : 


t l (4) 3CE) jtE) 
Su — de Epivap Ty BP == 5 TEyvog b\ by be 


LÉ 
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ou encore, le terme en b# b# s’annulant par anti- 
symétrie : 
Teuvap ES bE ER, 


Su 


ES 


Si nous exprimons b en fonction des trois 
autres vecteurs au moyen de (19) et contractons s, 
par l’un de ceux-ci, soit 9, il vient : 


(à) 


1 rst 0) 50) 1% 70) (8) 30) 
Su bu — = Teuves Eve et bn 0e BON D NES 


Cali 


ce qui fournit après quelques calculs : 


; 3 o RER : 
Su ED jure Ier A te = 6". 


Ainsi la projection du spin sur chacun des vec- 


teurs du genre espace du tétrapode est la constante 
correspondante sf — oi. Si nous nous souvenons 
que le spin est orthogonal à la vitesse unitaire nous 
voyons que le-spin est un vecteur d'espace propre 
rigidement lié au trièdre d'espace du tétrapode 
mobile. Naturellement sa grandeur : 


So = VSyu Su 


est aussi constante. 

Nous pouvons d’ailleurs particulariser la jauge en 
choisissant l’un des vecteurs, par exemple db 
sur $,, Ce qui est permis parce que la rotation néces- 
saire pour cela ést indépendante du temps (il res- 
tera encore un degré de liberté dans le plan des 
nouveaux b® et b®). Cela revient à choisir les 
deux constantes 61 et 6? nulles : 


5) RU) (3) 30) 

Comme D B% — 0 cela signifie que 0® est 
colinéaire au quatrième vecteur b®. En d’autres 
termes, en raisonnant dans le système propre, le 


d : 
vecteur d'espace “pe b% est nul, autrement dit : 


le spin est axe instantané de rotation dans le systéme 
propre. Les vecteurs d'espace b® et b% qui restent 
orthogonaux à l’axe fixe b(%, tournent autour de 
lui avec une vitesse angulaire Q, et on a 


in DM — (be) 
dr 


et 


— — (ht 
d 


Ainsi 


de (£ bu) bO 


(2 
n’est autre que le scalaire constant 


° _ (F4 po) = 6? IL. 


Aiïnsi le mouvement de rotation dans le système. 


propre est uniforme, le spin prend l’expression 
su = (by lu) be 
et sa grandeur est 
So io 
ce qui fournit comme conséquence du choix de la 
jauge, une relation importante (et tout à fait clas- 


sique) entre les constantes du mouvement [28]. 
Enfin, en utilisant l'expression du spin 


+ (3) 
RU 
nous pouvons exprimer le moment angulaire propre 


en fonction du spin et du balourd £,,, par la formule 


bien connue [6] : 
RE à Euvop du Sp — (lu À — I du) 
SOIÉ : 
SAS en envon re ED PE AUO DU De UE 
soit, en tenant compte des relations (19) : 
Guo = SE EEE MO AE RS 


On voit que c’est une généralisation de l’expres- 
sion considérée dans la section précédente. Du 
reste, si on contracte par ©,, — b® b®, on obtient : 


(1 6 1 : A 
T = 7 Su ous = 2 (so de due Ton bu). 


Le premier terme est bien le Lagrangien linéaire 
de la section précédente, le second est une consé- 
quence de l’existence d’un balourd. Nous voyons 
donc, en conclusion, quelle relation étroite il y a 
entre les deux formes proposées par le Lagrangien. 


V. Forme spinorielle des conditions de jauge. — 
Le modèle du tétrapode a été récemment relié par 
T. Takabayasi [18], [29] à la représentation spino- 
rielle habituelle des fonctions d’onde quantique. 
Soit 4 un spineur de Dirac, normé en posant : 


(po)? — (oys D)? = 1 avec Y— d+ y 


(+ est le conjugué hermitien de 4). On forme le 
vecteur du genre temps et le pseudovecteur du 
genre espace habituels, tous deux normés et ortho- 
gonaux l’un à l’autre : 


(4) _ 
ln — Ÿ Yu Ÿ 


(8) _ :T 
bu = 14 Ys Yu Ÿ. 
Puis on introduit les spineurs conjugués de 
charge d° en utilisant la matrice de conjugaison de 


charge C (C+ C—1, C+ = — C) qui transforme 
chaque matrice y, en sa transposée Mu 


T . 
Yaf= — C1 y, C 


No7 0 


N°7 
d’où 
= +- Cest VE Ce 
Les spineurs conjugués de charge sont alors : 
ge cuil C4, 
d’où : 
(G} = — gr ot 


On peut ensuite former deux nouveaux vec- 
teurs : 


f 


DD = (UT vu Ÿ — Y ru 
(24) 
(= CD vu Y + Dvu PT. 


On peut montrer facilement que ces nouveaux 
vecteurs sont tous deux du genre espace, qu'ils 
sont réels, normés, orthogonaux l’un à l’autre et 
aux deux précédents, si bien que nous avons ici un 
tétrapode semblable à celui que nous avons défini 


* dans la section III. Si une équation d’onde (ou 


plus généralement un Lagrangien) est donnée en 
fonction du spineur Ÿ, elle peut aussi être exprimée 
en fonction du tétrapode correspondant, ce qui 
peut mettre en évidence une interprétation phy- 
sique suggestive. 

Nous nous bornerons ici à montrer le lien des 
transformations de jauge mises en évidence dans 
le cas des deux Lagrangiens étudiés ci-dessus, avec 
celles qu’on considère généralement dans les équa- 
tions d'onde spinorielles. Une transposition sem- 
blable de la forme spinorielle dans celle du tétra- 
pode a déjà été proposée par l’un de nous [30] en 
associant des bases spinorielles et vectorielles parti- 
culières. En premier lieu, considérons la jauge spi- 
norielle la plus fréquemment utilisée, conservant 
par exemple l’équation de Dirac et l’ensemble des 
seize grandeurs tensorielles attachées au forma- 
lisme de Dirac, soit : 


be eir , 
r étant un nombre réel. bŸ et b® restent évidem- 


ment invariants. Mais on voit immédiatement que 
les spineurs conjugués de charge deviennent : 


ge — e—ir We y" = eir d° 
d’où on tire : 
MERE taire elie 
Ya He FEENC dan 117 4 
T1 /e [HE DAT EE ET 
Yu Y = 67 Ÿ Yu Y. 


Nos deux premiers vecteurs deviennent alors : 


‘ 2 5 

pe —= po cos 27 — AL sin 2r 
(9 : 2 É 
pie _ FA sin 2r + be cos 2r. 


Donc une telle transformation de jauge fait tour- 
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ner les vecteurs b®, b® autour du vecteur fixe b& 
dans l’espace propre. Nous retrouvons ainsi le 
groupe de jauge du Lagrangien linéaire de notre 
section IT. 

Considérons maintenant une transformation de 
jauge plus compliquée mettant en jeu la conju- 
gaison de charge, qui a été proposée récemment 
par Pauli [31] comme une transformation cano- 
nique linéaire d’un grand intérêt : 

g= ad + Pys ÿ° (25) 
d’où 
7” 2 k *k IC 
RE RE PEU UE 
où & et 8 sont deux nombres complexes quelconques 
vérifiant la relation : 
CA RO EEE 
ce qui ne laisse que trois degrés de liberté. Si on 


porte les équations (25) dans l’expression de D, 
celle-ci devient : 


el = (a 


Le calcul s’exécute sans peine en exprimant les 
spineurs conjugués de charge au moyen de la ma- 
trice de conjugaison C. Celle-ci peut sauter par 
dessus chaque matrice y en la transformant en sa 
transposée. Il faut remarquer que toute expression 
bilinéaire 4 À 4 est égale à sa transposée (4 À p)T 
qui symbolise un ordre différent des opérations : 


da de Cle 


Grâce à cette remarque, on peut montrer que : 


Ÿ ES B* 4° Y5) ua d HER BYS ÿ°). 


p7 AT QT. 


Pve P = Dre = rev 
q Ys Yu Ÿ — as 5 A g)T = 2h Ys Yu Ÿ = 0 
et de même 
Dove vu = (D) 5 vatU) = — ÿ vs vu Ÿ = 0. 
Finalement, il ne reste que : 
Eu = (ava + B*B) Ÿ ya Ÿ = de. 


Nous voyons que le vecteur D est invariant 
sous la transformation de jauge considérée. Au 
contraire, si on fait le même calcul sur bŸ, il vient 


D = (aa — B*B) LE + ifata — BYa) Li? 


+ (8 + BYa) bu”. 


Ce vecteur est évidemment orthogonal à bE et on 
voit facilement que : bf D = 1. 

Pour appliquer la transformation (25) aux vec- 
teurs (24), il faut utiliser une nouvelle identité 
qu’on démontre comme les précédentes : 


CT OT TT 


"Fe gere Nes AT 
Ÿ Yu Ys Ÿ = Ÿ Yu Y6 Ÿ = — V0 Yu Ys ÿ. 


x 


T * 
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Nous avons alors les transformations suivantes : 
a? ÿ° Yu Ÿ — By y Se 2aP À Vs Yu (} 

Et 1e Fe — LS 

D Yu D = XD yu D — BE D vu d + 2atB* DYs Yu Ÿ 
ce qui donne pour les vecteurs : 


RES 
S 
— 
F 
ne 
ll 


D = LL + pr + at +"Bt) ED 


Hg ler + port pa 1 — iap — a*p#) D 
ML Lun ge jme — 59 
(O2 — BA — ak G2) LD — (ap + a*p*) LU 


_ Sur ces expressions, on vérifie aisément les autres 


relations d’orthonormalité sur D, EP, D. 


Pour nous résumer, nous avons montré que la 


transformation de jaune spinorielle de Pauli expri- 
mée en fonction du tétrapode de Takabayasi est 


équivalente à une rotation tridimensionnelle des vec- 
teurs du genre espace b?, D, bP qui laisse b 
invariant. Il nous paraît intéressant de faire remar- 
quer que c’est là justement le groupe de jauge que 


nous avons mis en évidence dans le Lagrangien que 
nous avons étudié dans notre section IV, et qui a 


été utilisé par la suite pour la théorie des particules 
de Bohm-Hillion-Vigier. 
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VISCOSITÉ ET BIRÉFRINGENCE DYNAMIQUE DES MACROMOLÉCULES SOUPLES 


Par ANTON PETERLIN, 


Physikalisches Institut, Technische Hochschule, München, Allemagne. 


Résumé. — Si l’on introduit les valeurs moyennes de distances inframoléculaires, correspondant 
à une valeur donnée du gradient de vitesse, on peut calculer pour le modèle de l’haltère élastique 
la viscosité intrinsèque et la biréfringence dynamique dans tout le domaine du gradient de 0 à oo. 
La viscosité intrinsèque décroît d’abord jusqu’à un minimum puis croît lentement mais cons- 
tamment par suite de l’extension progressive de la pelote macromoléculaire, supposée parfai- 
tement souple. La biréfringence d'écoulement change moins dans ce cas que selon les résultats 
obtenus par la théorie simple qui ne tient pas compte de la variation de l’interaction hydrody- 
namique avec l’extension de la macromolécule. Les résultats présentés ne sont valables que pour 
des molécules infiniment longues et parfaitement souples. La longueur finie et la rigidité partielle 
des macromolécules réelles provoquent des effets de saturation qui vont être traités séparément. 


Abstract. — Intrinsic viscosity and dynamic double refraction can be calculated for all values 
of velocity gradient between 0 and co, in the case of an elastic dumbbell, when, for a given velo- 
city gradient, the mean values of intramolecular distances are introduced. Then the intrinsic 
viscosity first decreases to a minimum value and then slowly increases with velocity gradient, 
because of progressive stretching of the supposed perfectly flexible macromolecule. Dynamic 


- double refraction is less varying in this theory than in the simpler one which does not take into 


account variation of hydrodynamic interaction with stretching of the macromolecule. 
The results obtained are valid only for infinitely long and perfectly flexible molecules. 
Effects of saturation due to finite length and partial stiffness of macromolecules will be repor- 
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ted elsewhere. 


1. Pour l’étude des macromolécules linéaires en 
solution idéale (à la température 0 de Flory) le 
modèle de Rouse [1] est le plus approprié. La 
chaîne à Z chaînons, statistiquement indépendants, 
possède Z modes à l’aide desquels le mouvement 
de là chaîne et de tous ses chaînons dans l’écou- 
lement laminaire peut être déterminé. 

L’interaction hydrodynamique modifie les va- 
leurs propres de tous les modes, la modification 
dépendant essentiellement des distances actuelles 
entre les chaînons. Puisque celles-c1 elles-mêmes 
sont déterminées par les valeurs propres de tous 
les modes, le problème devient d'autant plus com- 
pliqué que jusqu'alors on n’a pas encore réussi de 
trouver une solution exacte et complète dans tout 
le domaine du gradient de vitesse G de 0 à co. 

Pour simplifier le problème on a d’abord rem- 
placé la distance actuelle entre les chaînons 7 et k 
par la valeur moyenne ;r (2,3) 


< Ar > =V6/r/hr 
= <> = bi — | (1) 


b — longueur du chaînon que l’on trouve en solu- 
tion au repos. Avec cette hypothèse la viscosité 
intrinsèque devient indépendante du gradient de 
vitesse [3]. Ce résultat contraire à l'évidence expéri- 
mentale est une conséquence de l’invariabilité des 
distances intramoléculaires moyennes (éq. 1). Il 


s'ensuit que ce ne sont que les valeurs initiales : 


de [1], An/G'et ctg 27 qui sont représentées correc- 
tement par le modèle de Rouse à interaction hydro- 
dynamique constante [3, 4]. 

La situation change si l’on se rend compte de 
l’extension de la pelote macromoléculaire dans 
l’écoulement. L’interaction hydrodynamique entre 
deux chaînons décroit avec G croissant par suite 
de l’accroissement de la distance mutuelle moyenne. 
En remplaçant < 1/r > par une fonction linéaire 


Ar SAR RE PENSE A APCE) 


et introduisant pour < r? > le carré moyen de la 
distance correspondant à la valeur actuelle du gra- 
dient de vitesse, Zimm [5] a pu montrer que la 
viscosité intrinsèque décroissait avec G dans tout 
le domaine de G 


[nla = [nlo(1 — a G? ee) (3) 


ce qui correspond, du moins qualitativement, à 
l'expérience. 

Le choix d’une valeur moyenne X en fonction du 
gradient néglige le fait que les distances intramo- 
léculaires changent deux fois pendant une rotation 
du maximum au minimum dans les directions res- 
pectives de dilatation et de compression maxi- 
males de l’écoulement laminaire. En outre il est 
difficile d’évaluer précisément l’influence de l’appro- 
ximation linéaire (éq. 2) qui sûrement empêche 
l'application de la méthode aux grandes valeurs du 
gradient. Il s’ensuit que de cette manière on ne 
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peut pas calculer avec une précision suffisante la 
viscosité intrinsèque et la biréfringence dynamique 
dans tout le domaine de G, de 0 à co, bien que les 
déviations initiales des valeurs constantes, que 
l’on avait obtenu auparavant par le modèle à l’in- 
teraction hydrodynamique constante, aient le sens 
correct. 

La situation est plus favorable avec le modèle de 
l’haltère élastique qui ne rend compte que du pre- 
mier mode du modèle de Rouse. Mais puisque ce 
mode représente exactement le mouvement de base 
de la pelote macromoléculaire dans l’écoulement 
laminaire — en moyenne les deux extrémités ont 
des vitesses exactement opposées et la même situa- 
tion se retrouve dans le mouvement des deux 
. moitiés de la macromolécule (fig. 1) — généra- 


F1G. 1. — Mouvement moyen des chaînons de la macro- 
molécule pelotonnée dans l’écoulement laminaire Hnéaire. 
Les vitesses des chaînons correspondants dans les deux 
moitiés ont même valeur mais sens opposé. 


lement les valeurs initiales de [1], An/G et ctg 2y 
calculées d’après le modèle de l’haltère ne diffèrent 
de celles du modèle de Rouse que par les cons- 
tantes multiplicatives. C’est pourquoi on s’attend 
à obtenir bien plus aisément à l’aide de ce modèle 
les traits caractéristiques du comportement de la 
pelote macromoléculaire dans lécoulement lami- 
naire sans avoir besoin d'introduire des approxi- 
mations trop vagues n’ayant aucune raison phy- 
sique et dues pour la plupart aux difficultés du 
traitement mathématique du modèle. 

2. En ce qui suit on adoptera les mêmes suppo- 
sitions qui ont été introduites par: Kirkwood-Rise- 
man [2] et Zimm [3, 4] sauf la substitution de 
< 1/r > par une fonction linéaire de h (éq. 2). 
C’est-à-dire, pour chaque valeur du gradient on 
introduira la valeur correspondante moyenne de 
distance réciproque < 1/r > sans se rendre compte 
de sa variation considérable pendant la rotation de 
la pelote. Pour les extrémités cette valeur s’obtient 
de la fonction de répartition, qui dans un champ de 
vitesse 

vi CLEO) (4) 
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s'écrit [6] 


Pre er, e— ul + 2B)a8—2 pay + ye+ (+ B9z2]/(1+ 62) 
ue — 3/2 18 = 3/2 Z (5) 


B = GJa uD = Min] nG/NAkT 


D = kT/n,W constante de diffusion des extré- 
mités, 10 viscosité du solvant, M masse molé- 
culaire de la macromolécule, VN; nombre d’Avo- 
gadro, k constante de Boltzmann, T température 
absolue. On a [7] 


: 1 = V2: F{(k, ©) (6) 
“ho Vvr e 
F'intégrale elliptique du premier ordre avec 
1 5 
REC EE (7 
| 1| Vi x) | 
Sin? ® — 6 


RVT EE 


En se servant de la fonction de répartition des 
chaînons (8) on obtient pour les distances intra- 
moléculaires réciproques moyennes 


1 CS: F(k, ©) 1 + 22/3 
ee be tr bV V2BVI + 2 m — 1 26? 
m m 84/1 + p2 ) 
Se TZ + 1 38 
(8). 


m = |] —k|. 


Si l’on compare les distances 
moyennes à différents m, on voit que l’extension 
de la pelote macromoléculaire n’est pas uniforme, 
les parties éloignées (m grand) s’étirant beaucoup 
plus que les parties voisines (m petit). 

L’interaction hydrodynamique modifie le coef- 
ficient de résistance moyenne de translation W [7] 

raZ 


We 


| AUeE 


* LUZ, 6 aVZ 
ms Tr b: 


E — [0,598 


ZI2 12 Z 

m a 

se | E (7) 22 
+ 1 


m=1l 


- 1 + 2623 
|) a — 
PTT tee _ 


réciproques " 


Nos 


a — rayon hydrodynamique du chaînon. Le rem- 
placement du coefficient 1,196 [9] dans e par le 
facteur semi-empirique 2/3 [9, 10] tient compte des 
approximations dans le calcul des moyennes. La 
distribution spatiale des chaînons est considérée 
par le facteur Z/6 [11] dans W qui remplace le 
facteur original Z/4 de Kuhn [12]. En accord avec 
les suppositions, à 6 constant, on ne se tient pas 
compte des variations de W pendant la rotation de 
la macromolécule. La valeur moyenne d’après 
équation 9 détermine le paramètre £ et par consé- 
quence la fonction de répartition ®(8). 

La variation de W en fonction de B est déter- 
minée par les facteurs S et f. En cas d'extension 
uniforme de la pelote macromoléculaire le facteur 


commun en S deviendrait 1/V/1 + 282/3 et la 
somme serait égale à 0,598. Il s’ensuivrait que le 
produit Ë — 5f/0,598 décroïîtrait et la viscosité 
intrinsèque, proportionnelle à W, croîtrait cons- 
tamment avec le gradient en désaccord complet 
avec l’expérience. Ce n’est que la non-uniformité 
de l’extension qui donne lieu à un caractère plus 
réel. Puisque la sommation en $ contient le fac- 
teur s(m/Z) de signe différent pour m petit et grand 


respectivement (fig. 2), le facteur commun, fonc- . 


10\ s(m/2) 


05 
(0) 2e Al m/Z 
SO 


4 


Fig. 2. — La fonction s(m/Z) de l’équation 9. 


tion de (m — 1)/(Z — 1) fait disparaître plus vite 
les termes négatifs, correspondant à m grand. En 
conséquence le produit & en fonction de $ com- 
mence par croître pour B petit jusqu’à un maxi- 
mum, décroît ensuite très lentement et finit par 
disparaître pour $ — co. 

3. La viscosité intrinsèque [7] peut être calculée 
à partir de £ (éq. 5) 


Nah°W _ Nabt raz 
EM = çm 6M, 1 +eë 


(10) 


M, masse moléculaire du chaînon. Pour une molé- 
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cule assez grande, Z — oo, la viscosité relative 
hunlessul etre edit le 
[nlret LE [nlo A 1 sr eË Fe a dE 1 Je (11) 


s'approche de 1 /E 
[nfrel, Z=00 = 1/6 = 1 — 0,1524 PB? + 0,0593-B4 ... (12) 
Normalement on trace [nlra sur G ou B, 


Bo = Mfn]0 0 ae g Lalo 


Na kT M]: 


Pour une comparaison avec l’expérience on doit 
alors remplacer 8 par B,. On obtient 


[nlrel, Z=co — 1 — 0,1524 B2 + 0,1058 B4 — ... 


(13) 


(14) 


La viscosité relative décroit d’abord jusqu’à un 
minimum puis croît assez lentement mais cons- 
tamment avec le gradient (fig. 3). La courbe com- 


3 15 l nl 


10 
05 
(0) 
OM 20 RU 0ER. 
Fic. 3. — Viscosité intrinsèque relative ([n]/[n]0)z=c0 


en fonction de f. 


mence par une tangente horizontale (éq. 14). Cette 
partie initiale disparait pratiquement dans la 
figure 3 à cause de l’échelle trop petite pour l’abs- 
cisse Bo. 

L’accroissement final de la viscosité intrinsèque 
résulte de l’extension de la pelote macromoléculaire 
dans l’écoulement. Si la molécule n’est pas parfai- 
tement souple et résiste au changement de forme 
(viscosité interne de Kuhn [13] et Cerf [4]) ou bien 
si la molécule à une longueur finie, ce qui est 
toujours le cas, l'extension devient plus petite que 
ne l’exige l’équation 5 et par conséquent F’accrois- 
sement de [nlra est fortement réduit et peut même 
disparaître complètement [14]. Tout de même il y a 
quelques mesures de viscosité, bien entendu en 


. solutions trop concentrées pour que l’extrapolation 


à dilution infinie soit encore possible, qui montrent 


27 


h: 


&* 
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une variation semblable à celle de la figure 3. Ce 
phénomène de « dilatation » a été trouvé par 
Selby et Hunstad [15] sur des solutions extrême- 
ment visqueuses (19 — 104 P) des polysiloxanes 
dans des huiles de graissage et par Merrill [16] sur 
des solutions de polyisobutylène dans le benzène. 
On doit attendre que des expériences ultérieures à 
concentration plus basse, permettant une extrapo- 
lation à la viscosité intrinsèque, confirment ou non 
la conclusion que cet effet est, du moins partiel- 
lement, dû à l’extension des pelotes macromo- 
léculaires dans un milieu d’autant visqueux que la 
résistance interne n'intervient plus sensiblement. 
4. La biréfringence dynamique résulte de l’aniso- 
tropie de l’indice de réfraction dans le champ 
d'écoulement laminaire. On a généralement 


FR RAR RE 
An — Vi — nèx) ? + Endy 


SEE + JE (ay — à) BVT + BE (15) 


RSA 3 M 


ARS 
Noy 

æ, et «, sont les polarisabilités optiques du chaînon 

dans la direction parallèle et perpendiculaire à l’axe 


du chaînon. La variation de l'interaction hydro- 


41 An 
(LS Eee / 


2 
Ur) 


F1G. 4. — Biréfringence dynamique 


An]C = (Andc) : [(4xf4a5n) (n° + 2)? Nales — «) [M] 
en fonction de $, pour le modèle à interaction hydro- 
dynamique constante (— —— — —) et variable( 
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. dynamique avec £ n’influe que sur le rapport B/8,, 


c’est-à-dire sur l’échelle des.abscisses. Les résultats 
d’après l’équation 15 ne diffèrent de celles du 
modèle de Rouse que par l’angle d’extinction où $ 
est à replacer par 0,4 8 [3]. Sur la figure 4et 5,ona 
tracé 


An 


: se ( + 2) at — à») C = An | 
or 3 M C 


et y en fonction de 8, et, à la fin de comparaison, 
les courbes respectives à l’interaction constante 


0 ; 50 100 5 pà 
F16. 5. — Angle d'extinction y en fonction de $, pour le 
modèle à interaction 


hydrodynamique constante 
(— — — —) et variable Je 


(8 — B,) qui correspondent aussi au cas sans inter- 
action hydrodynamique [6]. On voit que les courbes 
pour lesquelles linteraction est fonction de 8 
restent moins inclinées sur l’axe des abscisses, du 
moins au commencement, coupent ensuite les 
courbes à l’interaction constante à la valeur 8, 
pour laquelle [nl = [n]lo et à partir de là conti- 
nuent avec une pente plus grande. Sûrement on ne 
doit pas oublier que la longueur finie des macro- 
molécules réelles imite l'extension de la pelote et 
par conséquent la valeur de la biréfringence ne 
peut pas surpasser celle correspondant aux molé- 
cules complètement étendues. Ces effets de satu- 
ration vont être traités prochainement [17]. 

La valeur initiale de la biréfringence dynamique 
spécifique (constante spécifique de Maxwell) 

2] 


M 4 An = ef FE 
PP NC) em Goo 00 3n 


2 £ 
] (eu — a) El (16) 
ne diffère pas de celle qui l’on obtient sans consi- 
dérer la variation de l’interaction hydrodynamique 
avec l'extension de la pelote. On trouve la même 
situation pour l’angle d'extinction. La pente initiale 


— (5%) s- (5) A2 
TOfaen ANG / ice nee kT 10 


(17) 


Ne:7 


est proportionnelle à la viscosité du solvant. Il 
s’ensuit qu’à l’aide de ce modèle on ne peut pas 
reproduire l’effet spécial découvert par Cerf [18] et 


. étudié ensuite surtout par Leray [19], qui fait que 


la valeur limite de (dy /d Geo) n’est pas proportion- 

nelle à la viscosité du solvant 
dy A 

(le A Pre (18) 

Pour expliquer cet effet il faut introduire la 

notion de viscosité interne de la macromolécule 

comme l’a fait Cerf [4]. Mais il reste à voir comment 

ses résultats doivent être modifiés pour obtenir un 

accord satisfaisant en viscosité intrinsèque où la 


variation de l’interaction hydrodynamique, com- 


plètement négligée par Cerf, suffit pratiquement 
pour expliquer quantitativement la diminution ini- 
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tiale de la viscosité relative. A l’exception de la 
déviation initiale (premier terme dans la série en 
G?) les calculs présentés jusqu'alors par Cerf ne 


_peuvent pas être comparés aux courbes expéri- 


mentales. Même dans le premier terme conte- 
nant G? la variation des dimensions linéaires de la 
pelote influe sensiblement sur le coefficient : l’équa- 
tion 14 exige 0,1524, le calcul de Cerf une valeur 
entre 0,72 et 1,71 pour des molécules rigides 
[Kubhn [13] obtient la valeur 2] et entre 0,113 et 
1,436 fois FJZan,s pour des molécules souples. 
FJZan, est le rapport entre la résistance de la 
macromolécule à la déformation et la force de frot- 
tement par laquelle le solvant tend à déformer la 
pelote. 
Manuscrit reçu le 23 février 1961. 
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MESURE DIRECTE DE LA STRUCTURE HYPERFINE DU NIVEAU 6 ‘P, 
DE L’ISOMÈRE NUCLÉAIRE 1°’Hg* (1) 


Par C. BROT, 


Laboratoire de Chimie Physique de la Faculté des Sciences de Paris 


Résumé. — On décrit la mesure, par une méthode de double résonance récemment développée, 
de l’intervalle hyperfin F = 13/2 — F — 11/2 de l’état 6 °P, de l’isomère nucléaire 1°*Hg*, Cette 
méthode est basée sur l’observation du changement de la longueur d’onde de la lumière réémise 
après transition de radiofréquence. La précision obtenue est de l’ordre de 6.106, 


Abstract. — The hyperfine interval F — 13/2 — F — 11/2 in the 6 %P, state of the nuclear 
isomer 1°’Hg* has been measured using a recently developed double resonance method. This 
method depends on the observation of the wavelength change of the light re-emitted after the 
RF transition. The result obtained has a probable error of about 6 X 106, 


I. Principe de la mesure. — Le noyau 1°7Hg 
possède, en plus de son niveau fondamental de 
durée de vie 65 heures et de spin 1/2,.un niveau 
isomérique, habituellement dénoté 1?7Hg*, dont le 
spin est 13/2 et la durée de vie 24 heures. Dans 
l’état excité 63 P, de l’atome, cet isomère nucléaire 
présente donc trois composantes hyperfines dont 
les moments angulaires totaux sont Æ# — 11/2, 
F = 13/2 et F — 15/2. Des mesures purement 
optiques [1] ont montré que les séparations hyper- 
fines sont voisines de 14 KMHz pour 


F=MF = 13/2 


et de 18 KMHz pour F = 13/2 — F — 15/2. 
Plus récemment des valeurs plus approchées ont 
été obtenues par H. Hirsh [2] en utilisant des tech- 
niques de croisement de niveaux et de double réso- 
nance dans des champs magnétiques « intermé- 
diaires ». 

Parallèlement au travail de Hirsch, nous avons 
employé une méthode de mesure directe en champ 
nul ou faible, qui nous a permis d’obtenir une 


largeur Doppler: 


dune aie 


204 199 201 


700 201 798 
Fe FA F:Tfe F3/2 


-15000"1h3 0 


valeur précise de l'intervalle F = 13/2 — F = 11/2. 
Cette méthode, due à Kohler [3] fait appel elle aussi 
aux techniques si fructueuses de double résonance 
optique-radiofréquence proposées, et expérimen- 
tées par Kastler, Brossel et Bitter [4], [5]. 

Elle consiste dans son principe à détecter la 
transition radiofréquence au moyen du change- 
ment de longueur d’onde (et non de polarisation) 
de la lumière réémise par l’échantillon. C’est dire 
que son applicabilité nécessite : 

1) L’illumination sélective d’un seul des niveaux 
composants de la structure. 


À 


+ 


2) La détection sélective de la lumière, de lon-. 4 


gueur d’onde différente, réémise par celui des 
autres niveaux vers lequel on à induit la transi- 
tion de radiofréquence. 

Il faut donc disposer d’une radiation lumineuse 


excitatrice qui coïncide, dans les limites de sa - 


largeur Doppler, avec l’une et une seule des compo- 
santes de la raie étudiée. Il faut aussi que les sépa- 
rations des composantes soient supérieures aux 
largeurs Doppler. Cette dernière condition est réa- 


196 |199 207 
F3 F7 


197 
F-3/2 


+ 15000Mh g 


HTre de 


(1) Le présent travail a été effectué au Massachusetts 
Institute of Technology, Research Laboratory of Elec- 
tronics. 


lisée dans le cas des composantes hyperfines de la 
raie de résonance 19, — 3P, à 2 537 À des difté- 
rents isotopes du mercure. “ 


Na 


Pour réaliser l’illumination sélective de l’un des 
niveaux hyperfins, on pourrait songer à construire 
une source lumineuse de longueur d’onde ajustable 
en sélectionnant, à l’aide d’un quart d’onde et 
d’un polariseur, lune des composantes © de la 
structure Zeeman d'émission d’un isotope pair ; 
cela nécessiterait que la source soit soumise à un 
champ magnétique élevé réglable dans de larges 
limites. Une solution plus élégante [3] consiste, 
quand cela est possible, à profiter des circonstances 
suivantes : Le mercure possède de nombreux iso- 
topes et les déplacements isotopiques de leurs 
raies sont tels qu’il existe plusieurs coïncidences 
fortuites entre les composantes hyperfines pour les 
différents isotopes (fig. 1). 

C’est ainsi que, dans notre cas, nous avons pu 
profiter de l’assez bonne coïncidence qui existe 
entre la raie de l’isotope 196 et la composante 
F = 13/2 de 1°7Hg* : nous avons utilisé une lampe 
enrichie en Hg. 

Pour la détection sélective de la lumière réémise 
après transition de radiofréquence, il suffit d’absor- 


Re [Hg 1% } 


BOBINES 
D HELMHOLTZ 
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ber la lumière primaire directement diffusée : une 
cellule d'absorption contenant le même isotope 196 
remplit ce rôle. La raie réémise par le sous-niveau 
F = 11/2, se trouvant fort éloignée de toute autre 
raie des isotopes ses n’est absolument pas 
absorbée. 


IT. Techniques expérimentales. — L’échantillon 
de mercure 197* était préparé par bombardement 
d’une feuille d’or avec des deutons dans le cyclo- 
tron du M. I. T. (réaction 17Au, d, 2n). Le mercure 
produit était ensuite distillé sous vide dans une 
cellule de quartz en prenant les précautions décrites 
par Bitter et coll. [5]. Cette cellule était cylindrique, 
de dimensions voisines de 1 em x 4 cm. Le nombre 
d’atomes de mercure transférés dans la cellule, 
était, suivant notre estimation, de l’ordre de 
1013 — 1014, La méthode de préparation utilisée 
donne pour le mercure obtenu la composition sui- 
vante : 
17H o%* 17 9 : 1Ho 54 0 ; 


1%6Hg 17% ; 18Hg 42 %. 


CELLULE D'ABSORPTION 


(_Hg_196 


ff 


CELLULE DEF RESONANCE 


IG 0? 


Montage de mesure. — Le montage expérimental 
est schématisé sur la figure 2. La source lumineuse 
est constituée par une lampe sans électrode remplie 


d’un mélange de différents isotopes du mercure 


(fourni par Oak Ridge National Laboratory), 
enrichi à 8,4% de 1%Hg. La lampe est excitée par 
un magnétron fournissant 100 watts stabilisés 
vers 3 000 Mes. La lumière de la lampe est foca- 
lisée par un système de lentilles de quartz sur la 
cellule de résonance contenant l’échantillon. Celle-ci 
est placée au centre d’une cavité électromagnétique 
cylindrique percée suivant son axe d’un trou pour 
lPillumination. Des fentes de sorties de la lumière 
sont ménagées sur la paroi latérale de la cavité, 
parallèlement à ses bases, La cavité résonne dans 


FILTRE UV. 
PHOTOMULTIPLICATEUR \ 


Le 


un mode TE pour lequel, au centre, le champ ma- 
gnétique est axial et maximum. Elle est accordable 
grâce à un piston constituant l’une de ses bases. 
Son Q est de l’ordre de 6 000. Elle est alimentée 
par un klystron VA 92 C (délivrant à peu près 
200 mw) par l’intermédiaire d’une ligne en guide 
d’onde (standard KUÜ) comportant les organes de 
réglage habituels. Le champ magnétique hyper- 
fréquence au centre de la cavité est, à puissance 
maximum, de l’ordre de 1 gauss. La lumière dif- 
fusée perpendiculairement par la cellule est foca- 
lisée à travers les fentes de la cavité sur un photo- 
multiplicateur IP 28, avec interposition de la cel- 
lule d'absorption contenant elle aussi le mélange 
enrichi à 8,4 % de Hg 196, De plus un filtre inter- 


‘+ 
LS 
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férométrique UV, absorbant le visible, contribue à 
améliorer le rapport signal/bruit. 

L’énergie hyperfréquence est modulée en ampli- 
tude par un atténuateur à lame semi-circulaire 
tournant dans une section de guide fendu. Cet 
atténuateur est actionné par un moteur synchrone 
alimenté par le secteur. La modulation s’effectue 
ainsi à 30 Hz. A la sortie du photomultiplicateur 
la composante à 30 Hz est extraite par un ampli- 
ficateur sélectif et détectée dans un détecteur de 
phase (lock in) avec une bande passante de 0,2 Hz. 

La fréquence exacte du klystron est mesurée par 
battement dans un cristal avec l’harmonique con- 
-venable d’un fréquencemètre Guertsch FM 4 A 
lui-même accroché à un générateur à quartz 
Guertsch FM 6. L’oscillateur BF de celui-ci est 
lui-même étalonné par un fréquencemètre comp- 
teur Hewlett-Packard. La précision ultime de cet 
ensemble est de l’ordre de 107, Il à montré que 
linstabilité à court terme du klystron était de 
l’ordre de 35 kHz. 

Pour la commodité de la mesure le montage est 
complété par un balayage en champ magnétique 
faible permettant de mettre en évidence l'effet 
Zeeman sur la transition hyperfine (cf. $ suivant). 
Pour ce faire une paire de bobines d’Helmoltz est 
montée autour de la cavité de mesure, avec leur 
axe parallèle à la direction de la lumière incidente 
et au champ magnétique hyperfréquence. Ces 
bobines sont alimentées par une grosse batterie 
d’accumulateurs à travers un rhéostat dont le cur- 
seur est müû par un moteur électrique. Le champ 
magnétostatique varie ainsi de 0 à 50 gauss en une 
dizaine de minutes. Le courant dans les bobines 
est mesuré par une résistance étalon dont la diffé- 
rence de potentiel aux bornes est appliqué à 
l’entrée X d’un enregistreur XY présentant une 
bonne linéarité. Le signal de double résonance, 
prélevé à la sortie du détecteur de phase, est appli- 
qué à l’entrée Y de cet enregistreur. La différence 
de potentiel aux bornes de la résistance étalon est 
aussi mesurée sur un potentiomètre Rubicon : 
lorsque ce potentiel passe par une valeur prédé- 
terminée, l’opérateur envoie un top de marquage 
sur l’enregistreur XY, ce qui permet ensuite l’éta- 
lonnage des abscisses dans les enregistrements 
obtenus. 

Pour chaque valeur de la fréquence du klystron 
on prenait quatre enregistrements en inversant et 
le sens du courant dans les bobines (pour éliminer 
l'influence du champ terrestre), et le sens de par- 
cours du balayage (pour annuler l’effet de la cons- 
tante de temps du détecteur de phase). 


Résultats. — 4) MESURE EN CHAMP NUL. — Après 
avoir décelé la résonance autour de 14 235 MHz, 
nous avons relevé l'intensité du signal en fonction 
de la fréquence ; pour cela nous faisions varier 
point par point la fréquence du klystron en retou- 


chant à chaque fois l’accord de la cavité. Il fallait 
prendre soin évidemment de travailler à niveau de 
puissance constant, en réglant le klystron au milieu 
de son mode et en agissant au besoin sur un atté- 
nuateur. La courbe de résonance obtenue a une 
largeur voisine de la largeur naturelle de l’état 3P.. 
Son maximum f, pouvait être pointé à 0,2 Mes près. 


b) MESURE EN CHAMP MAGNÉTIQUE FAIBLE. — 
Les 14 sous-niveaux Zeeman du niveau F — 13/2, 
et les 12 sous-niveaux Zeeman du niveau F = 11/2 
sont sensiblement linéaires en champ pour les 
champs utilisés. 


Les transitions AF — 1 étaient, dans notre expé- : 


rience, du type Amy — 0, du fait que, dans le 
volume de l’échantillon (au centre de la cavité) le 
champ magnétique hyperfréquence était sensi- 


blement parallèle au champ statique. On a donc ” 


(fig. 3) à champ fixe 12 fréquences de transition 
(Mr = — 13/2 ... + 13/2). Lorsqu’à fréquence 
fixe on balaie le champ magnétique, il apparaît 
6 résonances qui se séparent ou non suivant que la 
fréquence f utilisée est éloignée ou non de f4,. 
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Pour faciliter l'identification rapide des raies 
enregistrées au cours d’une expérience nécessai- 
rement de courte durée, nous avons précalculé 
théoriquement, pour quelques valeurs de f — f,, le 


No7 LA 


| 
| 


N°7 


spectre attendu, c’est-à-dire la courbe théorique du 
signal en fonction du champ. L’Hamiltonien de 
Patome dans le champ 4€ s’écrit avec les notations 
usuelles : 


R = AI.J + BQo + gr us J.H + SP TH, (1) 


Les éléments de matrice de cet opérateur sont 
bien connus [6]. Dans la représentation [ZJ Fmy >, 
en négligeant pour l’instant la faible perturbation 
due aux éléments non diagonaux, l’énergie des 
sous-niveaux 7» est proportionnelle au champ. On 
obtient bien alors, par soustraction, pour la posi- 
tion des raies le diagramme linéaire et symétrique 
par rapport à f, tracée sur la figure 3 ; nous avons 
pour ce calcul adopté les valeurs numériques sui- 
vantes [2], les énergies étant exprimées par les 
fréquences correspondantes : 


ur — 1,399686 MHz/gauss 


gr = 1,4861 (au lieu de 1,5 si le couplage était pur Russel- 
Saunders) 


Wior — — 1,03 magnétons nucléaires. 


Pour calculer la forme d’un spectre il faut de plus 
faire intervenir la largeur des raies et calculer l’in- 
tensité relative de chaque raie. Celle-ci s’obtient 
de la manière suivante : on calcule d’abord le peu- 
plement relatif de chacun des sous-niveaux m de 
3P,, F = 13/2 produit par l’excitation lumineuse ; 
la lumière utilisée est sensiblement « blanche » dans 
le faible intervalle de la structure Zeeman, qui est 
très inférieur à sa largeur Doppler ; sa polarisation 
est o* + ©. Les probabilités de transition optique 
(transition du dipôle électrique P) correspondantes 
sont données par [<< «Fm|Plx' Fm + 1 >{?œet x’ 
symbolisant les états électrique LS, et P; respec- 


. tivement, F étant le même (13/2) dans les deux 


états. En explicitant ces éléments de matrice [6] 
on trouve que les probabilités de transition 
6° + o— et finalement les populations des sous- 
niveaux m»m sont proportionnelles à F(F +1) — m2. 

Les probabilités de transitions magnétiques de 
radiofréquence vers chacun des sous-niveaux My 
de F — 11/2 sont, en tenant compte de la restric- 
tion expérimentale Am — 0, proportionnelles à 
[<F, mJ|F—1,m >|?c'est-à-dire à F?— m?[6]. 
Nous avons donc pris pour intensité relative de 
chaque raie du spectre 


Jon = [F(F + 1) — m°] (F2? — mi) 


et calculé le spectre en admettant pour chaque raie 
une forme Lorentzienne en fréquence et une lar- 
geur égale à 2,3 MHz voisine de la largeur natu- 
relle : l’intensité du signal en fonction du champ 77 
est proportionnel à 


m= +11/2 


Jo m 
D à 
m=—11/2 1 | A (* en Pr =) 


2,9 


J = où Af = f—f0 
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et où les p,, sont les pentes des droites de la figure 3 
(MHz/gauss). Nous avons calculé ces fonctions pour 
trois valeurs de A. Celle correspondant à 


Af — 30 MHz est représentée sur la figure 4a. 

Les enregistrements expérimentaux étaient en 
bon accord avec les spectres théoriques. La figure 4b 
en donne un exemple. 


ÿ 
| 


RAT Ut 


Top de marquage A 


Fic. 4b. 


Le rapport signal-bruit, de l’ordre de 20, est très 
inférieur à celui obtenu par Kohler [3] pour Hg 
par la même méthode. Ceci n’est guère surprenant 
étant donné le nombre beaucoup plus faible 
d’atomes de l’échantillon qui participent à la 
double résonance, et les conditions d’excitation 
optique sélective bien moins favorables (8,4 % seu- 


te 
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lement de mercure utile dans la lampe). Malgré 
cela la mesure avec balayage magnétique reste 
légèrement. plus précise que la mesure en champ 
nul, car dans cette dernière ‘on ne peut être tout à 
fait certain que le changement de la fréquence du 
klystron et de l’accord de la cavité n’influe pas un 
peu sur le niveau HF dans celle-ci, | 
Nous avons pris des enregistrements tels que 
celui de la figure 4b à des fréquences différant de fo 
de + 10, + 20 et + 30 MHz à peu près. Pour 
chacune de ces fréquences nous avons pointé, par 
interpolation, pour celles des raies qui apparaissent 
bien séparées, le voltage V correspondant aux bornes 
de la résistance étalon, puis pris la moyenne V,(f) 
entre les 4 enregistrements correspondant à cette 
fréquence. Si les positions des raies étaient exacte- 
ment linéaires en champ, il suffirait, pour une raie 
donnée à différentes fréquences, d’extrapoler linéai- 
rement vers Vn — 0 pour obtenir /,. En fait même 
aux champs faibles utilisés, il est nécessaire d’intro- 
duire, pour le calcul définitif de f,, les termes en A7? 
dans les formules donnant les positions des raies en 
fonction du champ. Ces termes s’obtiennent en dia- 
gonalisant l’'Hamiltonien de l’atome par un calcul 
de perturbation. Leurs coefficients sont de l’ordre 
de 104 MHz/gauss?. Un étalonnage relativement 
grossier de V en fonction de Æ permettait de cal- 
culer ces termes avec une précision suffisante, étant 
donnée leur petitesse. Après soustraction de ces 
termes correctifs, nous avons calculé f, par extra- 
pelation linéaire à partir de différentes raies my. 


6) RésuLTars. — Les deux méthodes a) et b) 
concordent pour donner 


fo = Wr=i3>pr=u2 = 14 234,86 + 0,09 MHz. 


Il aurait été intéressant de mesurer avec la même 
précision le deuxième intervalle hyperfin, 


F = 13/2 > F — 15/2 
vers 18 250 MHz. L’auteur et H. H. Stroke ont 
récemment tenté, sans succès, de déceler cette tran- 
sition en utilisant le même montage alimenté cette 
fois par un klystron V 40 B. Cet échec était proba- 
blement dû : à la puissance insuffisante de ce klys- 
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tron ; aux condition très défavorables dans les- 
quelles était utilisée la cellule d’absorption : conte- 
nant, outre 1%Hp9, plusieurs isotopes stables, elle 
absorbait beaucoup trop la raie 15/2 réémise (fig. 1). 
Nous avons essayé également, sans plus de succès, 
de faire la transition dans le sens 


F=152-F= 132 


en illuminant le niveau F — 15/2 avec une lampe 
contenant un mélange 22Hg, *24Hg sous pression 
de gaz inerte. 

En combinant notre résultat ci-dessus avec ceux 
de Hirsch, on calcule pour 1°7Hg* : 


F — 13/2 => F = 15/2 = 18 248 + 11 MHz 
A = — 2 328,89 + 0,84 MHz 
B = — 902,9 + 5,4 MHz. 


Il est inutile, pour calculer le moment nucléaire 
à partir de À, de connaître exactement la fonction 
d'onde électronique de 3P, : On peut obtenir u97* 
par comparaison avec un autre isotope du mercure 
dont le moment a été mesuré directement, comme 
199Ho [7] : 


* ET A E*X 
U. 
Hi97  , 197 : 197 


H199 Lio A199 


À:99 étant connu [8], on obtient 
Uo7* = — 1,0316 + 0,008 


magnétons nucléaires. Ceci suppose évidemment 
que l’anomalie de la structure hyperfine soit négli- 
geable devant l’erreur probable sur À, ce qui est le 
cas présentement. Cette anomalie hyperfine, qui 
reflète la répartition du magnétisme nucléaire, 
aurait d’ailleurs été très intéressante à connaître 
dans un noyau de spin élevé comme celui qui nous 
intéresse. Malheureusement il aurait fallu obtenir À 
avec une précision supérieure, c’est-à-dire qu’il 
aurait fallu connaître le deuxième intervalle hyper- 
fin avec la même précision que le premier. 
L'auteur tient à remercier le P*' F. Bitter qui l’a 
accueilli dans son laboratoire et a dirigé son travail 
ainsi que le D' H. Stroke qui l’a constamment aidé 
de ses conseils. 
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DÉTERMINATION SPECTROSCOPIQUE DU SPIN DE 1’6Lu 
ET DES MOMENTS NUCLÉAIRES MAGNÉTIQUE ET QUADRUPOLAIRE DE !’5Lu ET !’Lu 


Par JEAN BLAISE, Jacques BAUCHE, Simon GERSTENKORN (*#) 
et 


Frank S. TOMKINS (*#), 
Laboratoire Aimé-Cotton, CG. N. R. $., Bellevue (Seine-et-Oise). 


Résumé. — Les structures hyperfines de raies de Lu Z et de Lu ZI émises par un échantillon de 
lutécium contenant 74,5 % de ’6Lu et 25,5 % de !’5Lu ont été enregistrées avec un spectromètre 
Fabry-Perot. Du rapport des intensités des composantes hyperfines, on a déduit que 1176 — 7, 
ce qui confirme le résultat de Kaliteevsky et Tchaika. Les rapports des moments des deux isotopes 
sont : 

175 [0176 — 0,7109 + 0.0050 et Q175/Q176 — 0.71 + 0.01. 


La valeur proposée par Steudel, 175 = + 2,0 + 0.2 un est confirmée, maïs la valeur du moment 


quadrupolaire doit être révisée : 


QU5 = (4.0 + 0.5).1/{(1 — R}.10—%4 


= (5.6 + 0.6).1/(1 — R’).10—% cm?, 


où À et R’sont les corrections de Sternheimer pour les électrons 6p et 5d. 
Les résultats sont comparés aux prévisions du modèle de Mottelson et Nillson, et le rapport 
des moments quadrupolaires indique que les déformations des deux noyaux sont identiques. 


Abstract. — The hyperfine structure of several lines of Lu Z and Lu 11, emitted by a sample 
of lutetium containing 74.5 % 476Lu and 25.5 % 175Lu, have been recorded using a Fabry-Perot 


spectrometer. 


that 2176 = 7, confirming the result of Kaliteevskü and Chaïka. 


of the two isotopes are : 


From the ratio of intensities of the hyperfine structure components, it is deduced 


The ratios of the moments 


Do ET 0 74090 0050 and 0 ir )OME=0 7422001" 
The value proposed by Steudel, u — + 2.0 + 0.2 unis confirmed, but the value of the quadru- 


pole moment must be revised : 


Q175 — (4.0 + 0.5).1f(1 — R).140—% = (5.6 + 0.6).1/(1 — R’) 10—% cm?, 


where R and R’ are the Sternheimer corrections for the 6p and 54 electrons, respectively. The 
results are compared with the predictions of the model of Mottelson and Nillson, the ratios of the 
quadrupole moments indicating that the deformations of the two nucleïi'are identical. 


I. Introduction. — Schüler et Gollnow [1], en 
1939, ont étudié à l’aide d’un interféromètre Fabry- 
Perot la structure hyperfine des raies 6 463 À 

(bd6s 3D, — 6s6p *P5) 
et 5 984 À 

(bd6s D, — 6s6p *P?) 
du Lu 1. Ils ne disposaient que de lutécium natu- 
rel contenant seulement 2,6 % de 176Lu et les inten- 
sités des composantes hyperfines étaient déduites 
d’enregistrements effectués sur plaque photo- 
graphique. Ils ont pu néanmoins établir que le spin 
de 176Lu était supérieur ou égal à 7 et inférieur à 20, 
que son moment magnétique p était à peu près 
indépendant de la valeur choisie pour le spin et 
égal à + 3,8 + 0,7 ux et que son moment quadru- 
polaire Q décroissait de + 8.0 à + 6.0.10—% cm?, 
suivant qu’ils adoptaient pour le spin les valeurs 7 
ou 13. 


(*) Ingénieur au Centre d'Études Nucléaires de Saclay. 
(**) Guggenheim Fellow (1960-1961), adresse perma- 


nente : Argonne National Laboratory, Argonne, Illinois. 


En 1957, Steudel [2] a repris l’étude de la struc- 
ture hyperfine de raies de Lu Z et Lu 77, mais cette 
fois à l’aide du spectromètre Fabry-Perot photo- 
électrique [3]. De cette étude faite aussi avec du 
lutécium naturel, il a conclu que la valeur la plus 
probable du spin était 6, les valeurs 5 et 7 n’étant 
pas exclues. Il recaleula les moments nucléaires de 
MELrLt 

175 = 2,0 + 0,2 un 
et 
OSEO EE O0 6) CM 
et à partir des rapports mesurés des moments des 
deux isotopes, 1l conclut que, pour 
RTE CORPS EROTSUN 
et | 
Q176 — (8 4 1).10—% em?, 
Disposant d’un échantillon de Lu,0O, conte- 


nant 30 % de 176Lu, Kaliteevsky et Tchaïka [4] ont 
refait en 1959 l’expérience de Schüler et Gollnow 


EX 
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sur la raie 6 463 À avec la même épaisseur d’étalon, 
mais en utilisant un spectromètre Fabry-Perot pho- 
toélectrique. La mesure du rapport des intensités 
des deux composantes extrêmes de la structure 
hyperfine de 1*6Lu leur permit de conclure que la 
valeur du spin de cet isotope était 7, tout en préci- 
sant que «il serait souhaitable d’effectuer des me- 


sures sur d’autres raies du lutécium pour éliminer 


complètement les valeurs du spin différentes de 7 ». 

L’un d’entre nous ayant obtenu d’Oak Ridge un 
échantillon de Lu,0, contenant 74,5 % de 176[u 
et 25,5 % de 15Lu, il devait être facile de réaliser 
ce souhait. 


IT. Dispositif expérimental. — «) SOURCE DE 
LUMIÈRE. — Nous avons utilisé une cathode creuse 
refroidie à l’azote liquide. Le gaz porteur était du 
néon et l’intensité du courant à varié de 10 à 
50 mA, selon la raie étudiée. La quantité totale 
d’oxyde de lutécium introduite dans la lampe a été 
d’environ 6 mg. 


b) SPECTROMÈTRE FABRY-PEROT PHOTOÉLEC- 
TRIQUE ENREGISTREUR. — Ce type d’appareil, en 
service au Laboratoire depuis 1948, a été décrit à 
plusieurs reprises [3], [5]. Le modèle utilisé com- 
portait un monochromateur à réseau du type Ebert- 
Fastie, un étalon Fabry-Perot à sept couches ZnS- 
cryolithe et un photomultiplicateur EMI du type 
9 558 B, sensible dans le rouge, couplé à un poten- 
tiomètre électronique du type Speedomax. On a 
exploré les anneaux d’interférence en laissant ren- 
trer l’air dans la cloche de l’étalon, vidée au préa- 
lablé, à l’aide d’un robinet à pointeau à réglage fin. 
Cette manière d’opérer complique un peu la mesure 
des positions des composantes de structure hyper- 
fine sur les enregistrements, car le temps mis pour 
explorer un ordre d’interférence augmente réguliè- 
rement par suite de la diminution de la différence 
de pression entre la cloche de l’étalon et l’atmo- 
sphère, mais elle laisse inchangée la hauteur des 
composantes et ne gêne donc en rien les mesures 
d'intensité dont dépend la détermination du spin. 


ITT. Détermination du spin de 176Lu. -— Comme 
Pont déjà remarqué nos prédécesseurs, en raison 
de la valeur élevée prévue pour le spin et pour le 
moment quadrupolaire, il est impossible de déter- 
miner le spin par application de la règle de Landé 
aux intervalles hyperfins ou d’après la multiplicité 
des termes hyperfins, cette dernière étant toujours 
limitée par la faible valeur de J dans les spectres du 
lutécium. La seule ressource est donc de mesurer 
les intensités relatives des composantes hyperfines. 
Cette méthode conduit à un résultat sûr quand le 
rapport théorique des intensités de deux compo- 
santes varie suffisamment avec le spin pour que, 
malgré les inévitables erreurs, le rapport mesuré ne 
donne qu’une solution. Ayant l’avantage de dis- 
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poser d’un échantillon considérablement enrichi | 
en 176Lu, il nous était facile de limiter nos mesures | 
d'intensité à des composantes choisies d’après le . 


critère précédent, parfaitement résolues et assez 


intenses pour réduire l'incertitude sur le fond pro- | 


venant des composantes voisines. Il fallait en outre 
s'assurer au préalable que le spectromètre Fabry- 
Perot permettait d'enregistrer correctement les in- 
tensités de composantes hyperfines. 


a) VÉRIFICATION DES INTENSITÉS DES COMPO- 


SANTES HYPERFINES FOURNIES PAR LE SPECTRO- 


MÈTRE. — Dans notre montage, l’étalon Fabry- 


Perot se trouve à la sortie du monochromateur. Si * 


les fentes d’entrée et de sortie de ce dernier sont 
prises égales, sa fonction d’appareil est triangu- 


Fe | 


Ines Bee 
EH 2 |: 
fente d'entrée) (fehte de sortie) Fit 
4h Fa Verre 


FC Te 


A 


laire (fig. 1). Le flux lumineux transmis au Fabry- 
Perot n’est maximum que pour la radiation 6, sur 
laquelle est réglé le monochromateur. Il en résulte 
que pour mesurer l'intensité de radiations 6,, 62. 


tion d'apparzil (triangle) 4 


voisines de 6,, il faut effectuer un balayage syn- 
chrone du monochromateur et de l’étalon. Le « 


balayage synchrone peut être réalisé de plusieurs 
manières [6], mais nous avons préféré adopter une 
solution très simple, bien qu’entraïînant une légère 
perte de luminosité. Les fentes d’entrée et de sortie 
du monochromateur sont prises inégales, ce qui 
conduit à une fonction d’appareil trapézoïdale 


2 = 
0 
F L_Fe 
(fente d'entrée) (fente de sortie): 


rc 2: 


(fig. 2) : toutes les radiations comprises dans le 
domaine 69 + 90, sont correctement transmises 
par le monochromateur à l’étalon. La figure 3 repré- 


4 
L 
| 


ee 


a  — 
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Fi. 3. — Enregistrement de la fonction d’appareil du monochromateur : 
b et c représentent les positions des composantes hyperfines de la raie À 6 463 À, 


a, 
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sente un enregistrement de la fonction d'appareil 
du monochromateur seul pour les raies 6 463 À 
de Lu ZI et 6 477 À de Lu J, réalisé, avec du luté- 
cium naturel, en faisant tourner le réseau. On a 
indiqué en pointillé la fonction triangulaire qui 
aurait été obtenue avec des fentes égales : on cons- 
tate que l’emploi d’une fonction trapézoïdale est 
nécessaire si l’on ne veut pas fausser les rapports 
d'intensité entre composantes de la structure de la 
raie 6 463 À. 

Pour nous assurer que nous ne commettions pas 
d'erreurs systématiques, nous avons enregistré les 
raies 6 463 À de Lu ZI et 4 905 À de Lu 1 avec du 
lutécium naturel (97,40 % de Lu) dont le spin 
est 7/2 et nous avons vérifié que les rapports d’in- 
tensité mesurés pour différents couples de compo- 
santes hyperfines étaient bien les rapports théo- 
riques (fig. 4a et 5a). 


b) MESURE Du SPIN DE 176Lu. — En substituant 
l’échantillon enrichi en 176Lu au lutécium naturel, 
les mêmes raies ont été enregistrées dans les mêmes 
conditions (mêmes étalons, mêmes distances entre 
les lames). Les figures 4b et 5b montrent des enre- 
gistrements types. Nous avons conservé la dési- 
enation des composantes adoptée par Steudel [2] : 
a, b, c, ... pour les composantes de Lu, g’, b’, 
c', ... pour les composantes correspondantes de 
176Lu et les figures 4c, 4d, 5c et 54 montrent la 
structure des raies et les schémas de niveaux hyper- 
fins correspondants. Au début et à la fin de chaque 
expérience, nous avons eu soin d'enregistrer le cou- 
rant d’obscurité (zéro de l’enregistrement). Les 
corrections de fond êt de pieds ont été soigneu- 
sement effectuées, soit par symétrisation de profils 
purs, soit par le calcul, en tenant compte des diffé- 
rents paramètres (effet Doppler, fonction d’Airy) 
qui contribuent à la formation du profil enre- 
gistré [6], [8]. Ces corrections s’avèrent peu impor- 
tantes, car nous avons choisi des couples de raies 
bien résolues, dont les intensités sont du même 


ordre de grandeur (les rapports mesurés varient 


de 1,3 à 2,6) et ces raies sont souvent perturbées de 
façon équivalente par les pieds des composantes 
correspondantes de Lu. A titre d'exemple, pour 
la raie 4 905 À, on trouve pour le rapport a’/f' les 
valeurs 1,49 sans corrections et 1,50 après correc- 
tions. Les résultats obtenus sont indiqués dans le 
tableau I où seul le dernier rapport mesuré nécessite 
la connaissance du rapport d’abondance des deux 
isotopes, 2,922 — 0,016 déterminé au spectro- 
graphe de masse. 

L'accord excellent entre les valeurs mesurées et 
les valeurs prévues théoriquement pour 7? —7 nous 
conduit à confirmer la valeur obtenue par Kali- 
teevsky et Tehaika [4] sur la seule raie 6 463 À. 
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TABLEAU I 


LoNGUuEUrR RAPPORT D'INTENSITÉS 
D’ONDE COMPOSANTES : 
(À) CALCULÉ MESURÉ | | | | k; 
— —- _ — = me ET 
: 6 1,364 D D 2 M 
6 463 _ 7 1,308 1,32 & 0,02 SL : 
OIL il ll 
8 1,267 DIE 
LE 6 1,600 ser EUR HA 
7 7 4,500 1,50 + 0,02 HT 1) PONE 
8 1,429 Ne EEE TE RE AE 
: Gore 1897. Le CLR 
4 905 D 7 41,808 1,82 + 0,02 ERA D a et 2. Le 
SAT as a 
Ê CRU rypiqe TT DR NE 
“ 7 2,550 2,57 # 0,03 DRM 
8 2,500 A Le : 
TE SO NEC ‘ 
Fa) l L $ | . 
JS | Nil i | 
RENE | Le 
eue (Ut 
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IV. Rapports des moments magnétiques et des 
moments quadrupolaires de !’’Lu et !’6Lu et véri- 
flcation du spin de !76Lu. — Outre les raies utilisées 
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pour la détermination du spin de {’6Lu, nous avons 
enregistré deux autres raies raies de Lu Z : 5 001 À, 
déjà étudiée par Steudel [2] et 6 004 À (fig. 6 et 7). 


O 20 28 


168,6 256 264 c— 


(10 êm") 
ei D: 


La raie 5 001 À (65? 6p 2P9,, — 652 7s 2S,%) per- 
met de calculer immédiatement le rapport des fac- 
teurs d’intervalles magnétiques A15/A176, 

(Dans tout ce qui suit les notations utilisées sont 
celles de Kopfermann [7].) 

On a en effet : 


GS7s S, 236,0 
3 
1 D ce 
POSE 2e API 
5 
6#6p 674 
è . RG 
37 Do 
Fic.7c. 


D'où l’on tire : 
A5 [A176 — 4,4917 + 0,0020. 


Pour la raie 6 004 À (652 Gp 2P9,, — 65? 7s 29,2), 
nous indiquons les résultats pour lisotope 175 seu- 
lement, l’enregistrement faitavecle mélange enrichi 
étant difficile à interpréter en raison de la proxi- 
mité des composantes &, b, cet a’, b’, c’ d’une part 
et des composantes e, f et e’, f’ d'autre part. 


TABLEAU II 


d— a 8 AUS 
FER UE 
SPECTRE TERME I URI PA) 


AUS CES) 


6sAGD 2 P1}, 42,48 + 0,08 
I 65 6p.2P2;. 7,43 + 0,30 
6sè7s 2, 59,00 + 0,10 
6nA6SED,, — 35,67 + 0,08 
IL 5d6s D, — 68,23 + 0,10 


BAIE ce) 


OUT) 
Ame) AU Cr) 


_- 29,91 + 0,08 — 
70,39 + 1,50 = Le 

pris 41,47 + 0,10 3 
6,39 + 0,40 —95,19 + 0,08 8,98 + 1,00 
28,64 + 0,10 — 48,05 + 0,08 40,51 + 0,10 
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Le tableau IT indique les valeurs des facteurs 
d’intervalles magnétiques À et des facteurs de cou- 
plage quadrupolaire B déduites de nos mesures. 

Il y a très bon accord avec les valeurs de Steudel 
pour lisotope 175 et avec celles de Schüler et 
Gollnow pour le terme 5d6s D; pour les deux 
isotopes. 

Dans le tableau ITT nous indiquons les valeurs des 
rapports A1%5/A17%6 et B17%5]Ba16, 


TABLEAU III 


SPECTRE TERME AO) ATIS DECIBES 
6s? 6p ?P?ys 1,420 + 0,006 
I 65? 7s Se 1,423 + 0,006 
6p?6s4Pyje 1,420 + 0,007 0,712 + 0,120 
I 5d6s D, 1,420 + 0,004 0,707 + 0,007 
Moyennes 4,421 + 0,003 0,71 + 0,01 


Les valeurs de A 1%/A17%6 pour les différents 
termes spectraux sont en accord avec la valeur 
1,4217 + 0,0020 calculée plus haut. En négligeant 
Panomalie de structure hyperfine, ce qui peut 
introduire, d’après Steudel, une erreur supplémen- 
taire de 0,5 %,, on en déduit le rapport des moments 


JOURNAL DE PHYSIQUE 


magnétiques des deux isotopes : 


le 1375 A175 


né = TR — 0,7109 + 0,0050. 


Le fait que le spin de l’isotope 175 est passé de 6. 
à 7 n’a réduit la valeur calculée par Steudel que . 
de 1/91. 
Pour le rapport des moments quadrupolaires des 
deux isotopes, on a : 


ONE D175 


QE puis = 0,71 + 0,01. 


L'étude de la variation du rapport B1®/B1%6 en w 
fonction de la valeur du spin choisi pour l’iso- 
tope 176 permet de vérifier la valeur déduite des 
mesures d'intensité. Ce procédé, déjà utilisé par : 
Speck et Jenkins [9] dans le cas de l’hafnium a été ” 
appliqué aux trois termes spectraux avec un J 
supérieur à 1/2, dont la structure hyperfine a été … 
mesurée pour les deux isotopes par Schüler et 
Gollnow, par Steudel ou par nous. On doit trouver » 
une valeur de 7176 pour laquelle la valeur de 
B175]B176 est la même quel que soit le terme spec- 
tral considéré. ki | 

Le tableau IV indique les valeurs de B175/B176 
calculées à partir des structures observées en attri- 
buant au spin de 176Lu les valeurs 6, 7 et 8. 


TABLEAU IV 


Rarporrs B175/B17,6 


176 1216 
> SCHÜLER- 6 SCHÜLER- 
TERME tn ee STEUDEL BELLEVUE Coins 
GPA GS EPS: 0,45 0,447 
2080 7/4 0.060 
5d6s D, 0,684 0,70 0,665 0,729 
+ 0,03 + 0,010 
6s 6p%P9 0,688 0,70 0,743 
ce 0,01 


On voit que pour les termes 3), et #P%, la varia- 
tion de B1%5/B1%6 en fonction du spin est lente et 
qu’il est impossible de choisir en raison de l’impré- 
cision des mesures. Pour le terme 4P,, par contre, 
la valeur du rapport varie très rapidement et bien 
que les erreurs de mesure possibles soient relati- 
vement grandes, il ressort clairement du tableau 
que c’est la valeur 7 qui donne les valeurs du rap- 
port les plus voisines pour les trois termes. Nos 
mesures ayant confirmé l’interprétation de la struc- 
ture hyperfine du terme 4P,, donnée par Steudel, 
il convient de réduire l’incertitude dont il a affecté 
. la valeur de B15/B176 dans l’hypothèse F = 6 
(0,45 + 0,30). À moins d'admettre que le terme 
#Py2 d’une part, ou les termes °D, et *P? d’autre 


F7 T8 

: SCHÜLER- : Si 

STEUDEL BE£ELLEVUE Riva es STEUDEL BELLEVUE 

0,697 0,712 1,30 1,39 

+ 0,120 + 0,40 
0,750 0,707 0,770 0,796 0,746 

+ 0,010 + 0,010 
0,760 0,795 0,817 


part, soient soumis du fait de termes voisins à des 
perturbations telles que le calcul des facteurs À 
et B n'ait aucun sens, nous pouvons conclure sans 
ambiguité à la valeur 7 pour le spin de l’isotope 176. 


V. Moments magnétique et quadrupolaire de 
FSLu. — Depuis que Schüler et Schmidt [10] ont 
confirmé que les déviations de la règle d'intervalle 
de Landé observées dans les structures hyperfines 
du lutécium s’expliquaient comme dans le cas de 
l’europium par l’existence d’un moment quadru- 
polaire, les moments magnétique et quadrupolaire 
de cet isotope ont été calculés à plusieurs reprises. 

Casimir [11] a utilisé les résultats expérimentaux 
de Schüler et Schmidt pour les termes 3D,,%,3 et 


No 7 


3P? 2 de Lu ZI ; Gollnow [12] a repris ces calculs à 
she des mêmes résultats expérimentaux auxquels 
il à ajouté le terme ?2D,, de Lu Z ; récemment, 
Steudel [2] a mesuré deux nouveaux “termes spec- 
traux et révisé les valeurs des données spectro- 
Scopiques entrant dans les formules de Goudsmit- 
Fermi-Segrè ; enfin le terme 6s? 6p 2P3,, de Lu 
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mesuré au cours de ce travail donne une nouvelle 
possibilité de calculer les moments. 


a) MOMENT MAGNÉTIQUE. — Nous indiquons 
dans le tableau V les différentes façons dont le 
moment magnétique a été calculé par Gollnow, par 
Steudel et par nous-mêmes. 


TABLEAU V 


FACTEUR 
AUTEURS TERME SPECTRAL D'INTERVALLE 
(TOC) 
{ BUS GS ED) Gdsfo — ‘ 4,62 
Case a 0) 
Lu 17 5d6s DRE 
3 ds b=—= 520 
GOLLNOW 
Apalo — 722) 
Gs6 3p à pal2 » 
A = 860 
Lu 7 Ed6s? Das. Aa 3 2 ra 6,41 
Do Ada/2 — 9,9 
GS 6p .2PUy dp so 42:30 
STEUDEL Gay à =, 10,6 
Lu ZI 5564 D, | 
RUE ge as = 320 
Ag —23024 
66 S 3p9 2 
| # Apilz — 94,3 
2P9)2 dpilz = 42,48 
BELLEVUE Lu I 6s?6p 
2PÈo Apsla — 7,43 


On peut constater l’excellent accord entre les 
résultats expérimentaux des différents auteurs, à * 
l’exception du terme 5d6s? 2D;,, dont la structure 
est très difficile à mesurer. 

Les facteurs d'intervalle des électrons @»1,2, Aya 
dasjo Ct Gage Sont égaux aux facteurs d’intervalle 
À des termes dans le cas du spectre d’arc ; dans le 
cas des termes des configurations 5d6s et 6s6p de 
Lu ZL,ils se determinent par la méthode du cou- 
plage intermédiaire de Breitt et Wills [13] en sup- 
posant que le rapport Gala est. égal au 
rapport théorique 


LH 3/2 Fruspe, I 0) À — E)}r-19 
t— 1/2 Fogyige [(1 — à) (1 — €)]1+a7a 


où les corrections à et e sont respectivement égales 
à 0,010 et 0,004 pour l’électron p,,: et négligeables 
dans les autres cas. 

Cette condition n’est pas respectée dans le dou- 


p.145 © 
CONSTANTES UTILISÉES VALEUR 
OPA Cie) ADOPTÉE 
Zi= "607; 0 Pi=004,9 3520 
Zi 0 OI 08 2,66 
na [Za®{AnalAn) = 2,20 2,28 2,6 
i À 
na [Za?(dna]An) — 2,79 2,85 + 0,9 
Zii= 610 = 4, 851,9 SRE 
Ë 1,84 
Zi — 60; SW = 1 993.9: 
î Fe 3,85 
Zi = 65005) 340,35 2,07 
1— 5 — 0,990 ; 1 — € — 0,996 
Zi = 60: SW — 2 440,35 2.48 2,0 
| na®|Za2{dna/dn) — 1,80 1,99 + 0,2 
EEE CE TC 
“ = 65; SV — 4 851,9 1,83 
152 8 = 0,990 :; 1 — = 0,996 
1 — D — 0,990; 1 —- € — 0,996 2,08 
PROMOS 00 
d © — (0 2,65 


blet 6$? 6p 2P9 et encore moins dans 5d6s? 21), où 
les facteurs d'intervalle ont été mesurés. direc- 
tement. Comme dans ces deux cas, le rapport me- 
suré est trop petit, 1l en résulte que les valeurs du 
moment magnétique déduites de @py7o OU Gaspa SOnt 
probablement trop petites, la conclusion inverse 
étant vraie pour les valeurs déduites de üp3,9 @t daspo 
par la même formule. 

Les deux termes de la configuration 6s6p ne per- 
mettent pas de vérifier la validité de la méthode de 
Breit et Wills. Par contre, à partir des termes 32,32 
combinés deux à deux, Gollnow a calculé pour a, les 
valeurs 313,4, 321,0 et 324,8 tandis que Steudel a 
trouvé 310,8, 318,9 et 322,9.10—3 cm1, On voit 
donc que le grand facteur d'intervalle a, est peu 
sensible aux variations possibles des petits facteurs 
d'intervalle &, quelles qu’en soient les raisons et on 
est ainsi conduit à préférer la valeur du moment 
magnétique déduite de ce facteur d'intervalle. Sa 


28 


Pol 
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détermination exige en outre la connaissance de 
l’expression rà(6s)/Z£(dn/dn) que Steudel a évalué 
à 1,80 par deux méthodes d’interpolation entre les 
spectres d’étincelle d'éléments voisins. Cette esti- 
mation est confirmée par la valeur 1,84 que nous 
avons déduite du schéma de niveaux de Lu 77 
récemment publié par Bovey et Pearse [14], en 
supposant que le potentiel d’ionisation de ce 
spectre est 115 000 cm1. Nous conclurons donc, 
comme Steudel, que la valeur la plus probable du 
moment magnétique de l’isotope 175 est 


+ 2,0 Ge 0,2 UN 
et celle de l’isotope 176, 
SE 2,8 + 0,3 UN: 


b) MOMENT QUADRUPOLAIRE. — Les différentes 
déterminations du moment quadrupolaire Q1%5 ne 
sont pas plus concordantes. Dans le tableau VI, 
nous comparons les valeurs déduites des divers 
termes spectraux par les mêmes auteurs que précé- 
demment. 


TABLEAU VI 


Q (10—?4 cm?) 
TERME SPECTRAL 


GOLLNOW  STEUDEL  BELLEVUE 
2 LA 
5d6s? D 312 9,9 
Lu Z Dj: 5,98 5,0 
6s?6p è 3/2 3,99 
SD, 5,60 5,49 5,49 
5d6s  3D, 5,83 5,30 
3 ECC 1 a. 
Lu ZI D, 5,98 6,16 
eee hLh 4,27 
PR. 4,0 3,85 


Gollnow et Steudel considèrent la configuration 
Gs6p comme perturbée et ils ont adopté comme 
valeur du moment quadrupolaire respectivement 
9,9 et (5,6 + 0,6).10-% cm1. Le fait que de 
6s? Gp ?P3,, nous déduisions une valeur de Q encore 
inférieure à celles calculées par ces auteurs à partir 
des termes 6s6p *P?,,, nous conduit à penser que 
dans les trois cas on se trouve en présence d’un 
effet de « shielding » de Sternheimer [15] dû à 
l’électron 6p. Que les termes des configurations 
od6s? et 5d6s donnent systématiquement des va- 
leurs plus élevées prouve simplement qu’on a un 
effet « antishielding » dû à l’électron 5d, car la 
dispersion des valeurs n’est pas moins grande et 
rien ne permet de les préférer au premier groupe. 

Le spectre d’arc étant un spectre à un électron, 
on peut, en supposant connu le moment magné- 
tique, déduire le moment quadrupolaire du rapport 


N°7 


des facteurs d'intervalle B/A des termes au moyen 
de la formule suivante [16] : 


De een) 


TRE 
Pr ae Noel 


P Rr 


Si l’on admet que & — 2,0 ux, on obtient les 
résultats indiqués dans le tableau VIT, où les fac- 
teurs d'intervalle des termes ?D sont empruntés à 
Steudel. 


TABLEAU VII 


5d6s? 
TERME SPECTRAL 6s?6p? Pi» 
3Dype De 
: 7,80 9,64 9,47 
€ 40% 3,185 1,300 1,355 
Q (16—% cm!) 6% 2,6 2,7 


Ce procédé de calcul du moment quadrupolaire 
qui est généralement considéré comme le plus 
correct puisqu'il permet d’éliminer l’incertitude 
provenant de l’évaluation de 1/r°, nous conduit à 
des résultats inconciliables avec les valeurs trou- 
vées précédemment et l’accord entre les valeurs 
déduites des termes 2D,7, et 2P3,, doit être consi- 
déré comme fortuit. 

Dans l’état actuel de nos connaissances sur 


l’effet Sternheimer, il est impossible de préciser da- . 


vantage la valeur du moment quadrupolaire et 
nous adopterons provisoirement les valeurs sui- 
vantes : 


OMS = (4,0 Ge 0,5) 4022 


1—R 


= (3,6 4 0,6) 


où À et À’ sont respectivement les corrections de 
shielding des électrons 6p et 5d. 


VI. Comparaison avec les prévisions théoriques. 
— Le modèle de Mottelson et Nillson [17] permet 
de prédire les spins et les moments des deux iso- 
topes du lutécium. 

Pour *#Lu, l'orbite [404 7/2] fournit l’état fon- 
damental Q — 7/2 + pour le 71€ proton, et, par 
suite, le spin 7/2. En supposant que l’excentricité 
du noyau à — 0,28, Nillson prévoit un moment 
magnétique de + 1,4 ux, soit 30 % trop faible. 

S1 Q est le moment quadrupolaire intrinséque du 
noyau, On à : 

a 
Q ER 


J 


5223 (1 es _ 2) 


No: 
où le rayon du noyau 
Ro = 1,2.10—1 A3 em. 
D'où : 
GITES 08: LOG Mes. 


Cette valeur est en excellent accord avec les va- 
leurs 7,6 et 8,5.10—% cm? déduites de mesures 
d’excitation coulombienne par Heydenburg et 
Temmer[18]et avec le comportement systématique 
des moments quadrupolaires que Townes [19] a 
décrit en traçant la courbe Q,/R$ en fonction du 
nombre impair de nucléons : pour Z — 71, cette 
courbe permet de prévoir Q,/R5 — 18, alors 
qu'avec la valeur 0,28 admise pour à, 


GUSIR — 18,1. 


Le moment quadrupolaire se déduit du moment 
intrinsèque par la relation : 


I(21 — 1) 


Qo- 
QI 088 75.16% 0m, 
15 
ce qui confirmerait l’effet d’antishielding de l’élec- 


tron 5d. À 
Pour 1*6Lu, le modèle de Mottelson et Nillson 
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prévoit que le niveau du 1052 neutron est l’orbite 
[514 7/21, soit Q — 7/2. — Il résulte de la règle du 
couplage sphéroïdal de Gallagher et Mosz- . 
kowski [20] que Z —Q, + Q, —7. Pourle moment 
magnétique, ces auteurs prévoient + 3,06 x, 
tandis que Hooke [21] prédit + 2,93 ux avec la 
théorie du couplage fort des noyaux impair- 
impair. L'accord de ces deux prédictions avec la 
valeur expérimentale 176 = + 2,8 + 0,3 ux est 
tout à fait satisfaisant. Il n’y a pas eu, à notre 
connaissance, de détermination directe du moment 
quadrupolaire intrinsèque Q675. Comme 


; 91 7 hrs 
QU — a 178 et Q175 = 1 Q'°, 
OA 
(D RE EME 
BEA 195 QE 


Avec la valeur mesurée Q176/Q17%5 — 170,71, il 
vient : O6 = 0,98 055. Ce résultat permet de 
conclure que les deux noyaux ont pratiquement la 
même déformation. 

Les auteurs sont heureux d’exprimer au Profes- 
seur Jacquinot leur gratitude pour avoir mis à leur 
disposition les ressources du Laboratoire Aimé- 
Cotton et leur reconnaissance pour lintérêt qu'il 
a pris à ce travail. 


Manuscrit reçu le 16 mars 1961. 
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SUR LES SPECTRES D'ÉMISSION DES SULFURES MIXTES 
DE ZINC ET DE CADMIUM ACTIVÉS A L'ARGENT : 
EXISTENCE DE DEUX BANDES D'ÉMISSION 


Par M. Hucues PAYEN DE LA GARANDERIE, 


Laboratoire de Luminescence, Faculté des Sciences de Paris. 


Résumé. 


Les ZnS et CdS activés au cuivre et à l’or présentent deux bandes principales 


d'émission. Il en est de même des CdS activés à l'argent. Seuls les ZnS(Ag) paraissaient faire 
exception. Le présent travail confirme des indications antérieures d’Henderson et collaborateurs 
sur l’existence des deux bandes également dans les ZnS(Ag). Ces bandes ont été suivies dans les 
sulfures mixtes ZnCdS(Ag), ce qui en a permis l'identification précise (1). 


Abstract. — Two main emission bands are observed in the spectra of the Cu- and Au- activated 


zinc and cadmium sulphides. 


The same is true for CdS(Ag). 


But ZnS(Ag) appears as an excep- 


tion. . The present work gives evidence for the occurence of both bands in ZnS(Ag) also, supporting 


previous indications by Henderson and co-workers. 


The displacement of these bands has been 


studied in the mixed crystals ZnCd$S(Ag), leading to their accurate identification. 


- L. Introduction. — La plupart des sulfures lumi- 
nescents photoconducteurs présentent deux bandes 
prificipales d'émission. C’est le cas du sulfure de 
zinc activé au cuivre ZnS(Cu). Généralement la 
plus intense de ces bandes à son maximum dans le 
vert, à une longueur d’onde voisine de 0,530 u. 
Mais à basse température on peut en observer une 
deuxième dans le bleu, dont le maximum est à 
0,450 y environ ; cette bande est souvent éteinte 
dès la température ordinaire et en tous cas dispa- 
rait 100 ou 1500 plus bas que la bande verte. 


Banpgs D’ÉMISsION DE ZnS(Cu). 
TEMPÉRATURE ORDINAIRE 
BANDE VERTE BANDE BLEUE 


Zn$S wurtzite 5 230 À & 450 À 
(« gap» = 3,70 eV) 
Zn$ blende 5 350 À a 600 À 


(«gap — 3,64 eV) 


Dans certains produits (haute concentration en 
activateur et absence d'oxygène dans l’atmosphère 
de calcination) on peut rendre la bande bleue pré- 
dominante [1]. 

Il y à quelques années encore, certains auteurs 
se demandaient si la bande bleue des ZnS(Cu) était 
réellement due à l’activateur cuivre et sil ne 
s'agissait pas de la bande bleue des ZnS « self- 
activés », attribuée à des lacunes [Zn]. Les lon- 
gueurs d'onde de ces deux bandes à température 
ordinaire sont en effet voisines ; mais lorsqu'on 


(4) Certains des résultats contenus dans ce travail ont été 
déjà présentés au Colloque sur les centres colorés et la 
luminescence cristalline, organisé à Turin par le Pr G. Bonfi- 
glioli (5-7 septembre 1960). 


abaisse la température elles se déplacent en sens « 
inverses [2], [3]. | 


TEMPÉRATURE AIR LIQUIDE 
ORDINAIRE 
Bande bleue 4 450 À 4 380 À 
du Zn$(Cu) 
Lacunes [Zn] 4 520 À 4 560 À 


Deux bandes ont été aussi observées pour le sul- 
fure de cadmium activé au cuivre CdS(Cu) et pour 
le sulfure de zinc activé à l’or ZnS(Au). E. Grillot 
en 1950 a montré l'existence d’une deuxième bande 
pour le sulfure de cadmium activé à l’argent 
CdS(Ag) [4]. 

Rien d’analogue n’avait pu être obtenu pour le 
sulfure de zinc activé à l’argent ZnS(Ag), dont le 
spectre d'émission montre un seul maximum dans 
le bleu. Toutefois Henderson, Ranby et Halstead [5] 
ont obtenu une indication en faveur de l’existence 
d’une deuxième bande de plus grande longueur 
d'onde. Admettant que s’il y avait un seul type de 
centres luminogènes introduits par l’argent dans le 
réseau cristallin on aurait une répartition gaus- 
sienne de l’intensité lumineuse /(v) en fonction de 
la fréquence v, ils trouvent que l'expression 


I(v) = 1,e—xv—v} 


ne représente pas bien le spectre observé. Ils y 
parviennent par contre avec une formule du type : 


I(») = À, eau) LI, eat, 


Le premier terme représente la bande principale 
avec un maximum à la longueur d’onde À;, le deu- 


NOT 


xième une bande moins intense avec un maximum 
pour la longueur d’onde X,. 


BANDES D’ÉMISSION DE ZnS(Ag) 
D'APRÈS HENDERSON [5]. 
TEMPÉRATURE ORDINAIRE 

BANDE PRINCIPALE BANDE SECONDAIRE 


4 760 À 
5 000 À 


4 370 À 
4 520 À 


Zn$S wurtzite 
ZnS blende 


Contrairement à ce qui se passe pour les ZnS(Cu),. 
la bande bleue-violette des ZnS(Ag) wurtzite se 
distingue aisément de la blande bleu clair des 
lacunes [Zn] ; d’ailleurs elle se déplace vers les 
courtes longueurs d’onde par refroidissement, con- 
trairement à celle des lacunes [Zn] [3]. 

Même en modifiant la température de calcination 
ainsi que les concentrations en activateur, ou en 
fondant (présent ou absent selon les essais), nous ne 
sommes pas parvenus à modifier appréciablement 
amplitude de la bande secondaire : dans les 
ZnS(Ag), celle-ci semble n’apparaître que comme 
un faible élargissement de la bande principale du 
côté des grandes longueurs d’onde. 

Toutefois, nous pouvons, semble-t-il, confirmer 
existence de cette deuxième bande : 


19 Nous avons fait une étude expérimentale des 


spectres d'émission des sulfures mixtes ZnCdS(Ag), 
avec des proportions variables de ZnS et de Cds. 
Nous avons ainsi pu suivre de proche en proche les 
deux bandes À, et à, depuis ZnS(Ag) exempt de 
CdS jusqu’à CdS(Ag) exempt de ZnS, où nous re- 
trouvons les résultats de E. Grillot [4]. Or pour 
CdS(Apg) l’existence de deux centres Ag distincts ne 
saurait faire de doute : l’émission de grande lon- 
gueur d’onde 2, est prépondérante dans les sulfures 
calcinés à haute température 900 ou 1 000 °C, 
lémission de plus courte longueur d’onde à, est 
prépondérante dans les sulfures caleinés à basse 
température 600 ou 700 0C [4]. 

20 Chacun des sulfures mixtes ZnCdS(Ag) a été 
étudié à température ordinaire 20 0C et à plus 
haute température, généralement 105 °C. Dans ce 
dernier cas la bande À, est relativement moins 
éteinte que la bande à, de même que dans les 
ZnS(Cu) la bande verte s’éteint moins rapidement 
que la bleue. Nous avons effectivement constaté 
que les spectres se déforment très peu du côté des 
courtes longueurs d’onde quand la température 
s'élève, et beaucoup du côté des grandes longueurs 
d'onde. | 


II. Spectres dus aux niveaux localisés de l’acti- 
vateur dans la bande interdite du réseau cristallin. 
— Un sulfure luminescent photoconducteur pré- 
sente deux bandes d’énergie permise, la bande de 
valence et la bande de conductibilité, séparés par 
une bande interdite (fig. 1). | 
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L’activateur introduit des niveaux localisés dans 
la bande interdite. Par analogie avec ce qui se 
passe dans les ZnS(Cu), nous admettrons que le 
niveau fondamental du centre Ag est situé à envi- 
ron ! eV au-dessus de la bande de valence, tandis 


Bande de conductibilité 


Activateur 


Bande de valence 


Fi@. 1. — Schèma d’un centre luminogène. 


que le niveau excité est extrêmement proche du bas 
de la bande de conductibilité. D’autres centres d’un 
modèle différent peuvent exister également (modèle 
de Lambe et Klick [6]) ; nous y reviendrons. 

Sous l’action de la lumière ultra-violette, le 
centre Ag passe dans l’état excité, ce qui se traduit 
en définitive par son ionisation. [’émission lumi- 
neuse a lieu lorsqu'un électron de la bande de con- 
ductibilité se trouve trappé par le niveau excité, 
puis retombe sur le niveau fondamental. La chute 
directe de l’électron dans le niveau fondamental à 
partir de la bande de conductibilité est vraisem- 
blablement un événement rare : sinon les spectres 
d'émission dépendraient beaucoup du mode d’exei- 
tation employé, ce qui n’est pas le cas. 


ÉQUATION DU SPECTRE D'ÉMISSION. — D’après ce 
qui précède, nous pouvons accepter pour décrire les 
spectres d’émission dans un sulfure photoconduc- 
teur la même formulation que celle utilisée pour les 
KCI(TI). L’interaction entre le centre et le réseau 
est à l’origine de l’élargissement du spectre en une 
bande d’émission. 

Soit [u}? le carré de l’élément de matrice du 
moment dipolaire et Z(v) l'intensité lumineuse pour 
la fréquence v. On a : 


la sommation étant à effectuer sur tous les états 
excités vibrationnels du centre luminogène à partir 
desquels l'émission d’un photon Av est possible, Ces 
états sont décrits par des fonctions d’onde d’oscil- 
lateur harmonique d'(r —r,) et d'(r) : 


+00 
= Const. sk de(r— r5) Y(r) dr 
—00 


r décrit la configuration des ions en vibration et r 
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le décalage de la position d’équilibre des ions entre 
l’état fondamental et l’état excité du centre. Dans 
les ZnS, il s’agit vraisemblablement, en première 
approximation, de la distance au centre des ions 
soufre ou du coactivateur voisins du centre lumi- 
nogène. 

Un calcul classique (voir par exemple [7]) montre 
que la sommation 


G(v) = 2 ul? 
conduit à un facteur sensiblement Gaussien : 


4 ke hk2(v — v,)?. 2 hy 
G{(\) = exp [5 e { 0) , FRE re . 


tanh 
2 k? rè 


Dans cette formule, v. représente la fréquence de 
vibration du centre dans l’état excité et ke est 
relié à v. par 


27Ve =V ke] M 


M masse des ions en vibration. De même y est la 
fréquence de vibration du centre dans l’état fonda- 
mental et l’on a 


27rwy = Vk/[M. 


Enfin T est la température absolue et X la 
constante de Boltzmann. 
Nous poserons pour simplifier : 


G{\) =— ga v— Vo)? 
et 
Ga vi 
I(») Tes G{»). 


Bien que le facteur v{ varie de 1 à 16 d’une 
extrémité à l’autre du spectre visible, souvent les 
spectres expérimentaux se laissent bien représenter 
par une formule du type : 


I(s) &: VIF e—a(v— vi)? 


à condition qu'ils ne comprennent qu’une seule 
bande. 

Si le spectre comprend deux bandes d'émission, 
on doit avoir pour /(v) l'expression : 


I(V) = Aj vert 2 4, vhenalv=u)? 
ou la relation approchée : 
FAST EME TS Er AE 


L’approximation gaussienne pour G(v) peut 
n’être pas valable ; c’est le cas quand l’état excité 
du centre et l’état fondamental correspondent à des 
configurations des ions extrêmement voisines (r, 
petit). Dans ces conditions, à un centre donné 
correspond une bande unique dissymétrique et la 
résolution de cette bande en composantes gaus- 
siennes n’est pas permise. Un exemple bien connu 
est celui de la willémite Zn,Si0,(Mn), étudié par 
Klick et Schulman [8]. Mais le cas de l’activateur 
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manganèse (retournement du spin d’un électron) est 
très particulier et il ne semble pas qu’il y ait lieu de 
mettre en doute la validité de l’approximation. 
gaussienne dans les sulfures phosphorescents du. 
type ZnS(Cu) ou ZnS(Ag) [9]. 


PARAMÈTRES DES SPECTRES D'ÉMISSION. — Les 
spectres d'émission construits expérimentalement 
donnent l'intensité lumineuse Z(à) en fonction de 
la longueur d’onde À. On a évidemment : 

I(x) dx = 1(v) dv 
“d’où 
I(à) = 1(v) cp? 

Par suite, pour une bande d’émission unique, 

nous utiliserons expression : 


I) = A(1/X) se GEAR 
La longueur d’onde 2 du maximum du spectre : 
AIX) [AR = 0 
est plus courte que À : 
RARES LORIE 


Mais en fait l’écart est faible et varie suivant An 


de 5 mu à 15 mu. 


IIT. Spectres expérimentaux des sulfures mixtes | 
ZnCdS (Ag). — Nous avons réalisé une gamme de 
sulfures du type 


æ CdS, (1 — x) ZnS (Ag) 


x proportion (en poids) de CdS variant de 0 à 
100 %. Nous donnons ci-dessous les spectres obte- 
nus pour quelques-uns d’entre eux. 

Tous sont activés par 2.101 Ag et calcinés à 
650 °C pendant une demi-heure, en courant d’azote. 
Cette basse température de calcination favorise 
d’après E. Grillot [4] l'apparition de la bande de 
courte longueur d’onde, s'étendant jusqu’au visible, 
des CdS(Ag) à luminescence normalement infra- 
rouge. 

Le Zn$ était donc sous forme blende, tandis que 
les ZnCdS sont sous la forme wurtzite : on sait en 
effet depuis les premiers travaux de A. Guntz [10] 
que dans les sulfures mixtes le CdS impose la eris- 
tallisation dans le système hexagonal. Aucune dis- 
continuité ne s’observe d’ailleurs entre la blende et 
la wurtzite, aussi bien dans le déplacement des pics 
des bandes d'émission que dans la variation de la 
largeur de la bande interdite (Krôger [11]). 

L'introduction du fondant (3 % NaCI) ne modi- 
fiait pas les positions des bandes d'émission, mais 
par contre faisait varier en général l'intensité de la 
luminescence émise. Pour ZnS, surtout lorsqu'il est 
calciné à une température aussi basse que 650 0C, 
le fondant joue un rôle favorable. Pour CdS par 
contre, nous trouvons comme l’avait déjà indiqué 
Grillot que son introduction est plutôt néfaste, 
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150 
nes 
100 
5 
0,400 0450 0,500 0,550. 0,600 0,H00 0,450 500 0,550 0,600 0650 


À {micran) x (micran) 


0,450 0,500 0,550 0,600 0,650 (micran) 0,500 0,550 0,600 0,650 0,700 


0,500 0,550 0,600 0,650 000 {micron) 
Fc. 2. — Spectres d'émission de ZnS(Ag) et des sulfures mixtes ZnCdS(Ag) 
renfermant respectivement 0,20, 40, 60 et 80 % de Cds. 


Courbe I: Spectre à 20 °C. É 

- Courbe II: Spectre à 105 °C, amené à la même ordonnée maximum. 
Courbe III : Courbe IL — courbe I, indiquant le profil approximatif de la bande secondaire, 
Concentration en activateur 2.10—4 g Ag/g sulfure. Calcination à 650 °C, 
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La figure 3 indique le déplacement de la longueur 
d’onde du maximum des deux bandes Any et Amo 
lorsque x varie. 

Pour la bande principale vers le bleu, l’on re- 
trouve en gros les résultats déjà connus. Plus préci- 
sément, nos résultats sont très voisins de ceux de 
Rothschild [12] tant que la proportion de CdS ne 
dépasse pas 60 à 70 % ; au delà, les auteurs qui 
n’indiquent qu’un seul pic donnent évidemment à 
celui-ci une position intermédiaire entre An et Amo- 

Pour obtenir simplement la position du maxi- 
mum de la bande secondaire de grande longueur 
d’onde, nous avons tracé les spectres Z(à) à diffé- 
rentes tempéiatures, par exemple à la température 
ordinaire 20 0C et àlatempérature 105 0Coù l’extinc- 
tion thermique est déjà importante, surtout pour la 
bande principale. Les deux spectres ont été alors 
ramenés à la même hauteur maximum, comme sur 
les figures 2, et nous avons admis que la différence 
des ordonnées des deux courbes donnait le profil 
de la bande secondaire. Cette façon de faire néglige 
les variations de position et de largeur de la bande 
principale lorsque la température s’élève ; mais elle 
parait approximativement justifiée si l’on remarque 
la coïncidence des deux spectres du côté des courtes 
longueurs d’onde. 


0,800 


Âmp 
À (micron) D 
4 4 Am, 
0,700 7 


0,400 RS re En De 0 ee 
0 10 20 30 40 50 60 70 80 90 100 


F1G. 3. — Déplacement des maxima des bandes d'émission 
pour une gamme de ZnCdS(Ag). 
æ pourcentage en poids de Cds. ; 
Concentration en activateur 2.16—1 g Ag/g sulfure. 
Calcination à 650 0C. Am, bande principale ; Am bande 
secondaire. 


Une analyse plus fine des spectres nécessite une 
précision que ne permettent guère la plupart de nos 
mesures : dès que la proportion de CdS dans le 
sulfure devient notable, une grande partie du 
spectre est émise dans le rouge, où la dispersion du 
spectrographe utilisé devient faible. Toutefois nous 
avons effectué cette étude sur le ZnS(Ag) à 20 °C, 
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D’après ce qu’on a vu précédemment, nous 


ASCEA || 


N9:7 


avons cherché à représenter la répartition de lin-. 


tensité lumineuse émise par une expression de la 
forme : 


IAE AN 
1 = At pe e GTR) 4 no e7GTR) 


es) 
150. 


0,50 0,500 


F1G. 4. — Spectre de ZnS(Ag). 
En trait plein : spectre expérimental. 
En pointillé : analyse du spectre en deux bandes, 
d'équations respectives ! | 


0,400 


5 LP 2 
TO) = (AIN) UC Es) 
(courbe TI, bande principale) et : 
LR iNS 
L(n = (4:1X) 5 CG 0,500/ 


(courbe IT, bande secondaire). 


Concentration en activateur 2.10—4g Ag/o ZnS. Tem- 


pérature de calcination 650 °C, 


Nous y sommes parvenus avec les valeurs sui- 
vantes des paramètres : 


A1 = 0,458 pu CARE 
As = 0,500 p to — 64,5 u? 
ACTA OMES 


Ces valeurs de à, et à, correspondent aux maxi- 
mums à, des éléments de matrice G(v). Les posi- 
tions des pics An des bandes d'émission, que l’on en 
déduit au moyen de la formule 


3 9 2 
A0 — Âm — 3 Xm/( 0 RS 3m) 
sont respectivement 


Ami = 0,450 4: Ame = 0,495 y. 


On comparera aux valeurs obtenues par la mé- 
thode simple précédente et figurant sur le gra- 
phique de la figure 3 : 


Am 1 = 0,450 te Am 2 = 0,485 HL. 


Pour la bande secondaire la. précision est faible 
10 my, mais néanmoins suffisante pour affirmer 
5 SAIe 
qu'il ne s’agit pas d’un effet de contamination par 
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le cuivre (position de la bande verte du cuivre dans 
ZnS blende : 5 350 À). A 105 0C l'émission étant 
fortement réduite du fait de l’extinction thermique 
il est difficile de trouver les paramètres relatifs à la 
deuxième bande. Pour la bande principale nous 
donnerons : 


A1 = 0,465 


On observe donc un élargissement de la bande, 
« diminuant avec la température, ce qui est con- 
forme à la théorie. Cet effet est accompagné d’un 
léger décalage vers le rouge. Ces deux phénomènes 
s'ajoutent pour déformer la partie du spectre située 
vers les grandes longueurs d’onde, tandis qu’ils se 
compensent du côté du violet. C’est pourquoi la 
méthode simple utilisée ci-dessus (fig. 3) ne donne 
que des résultats approchés. Toutefois l'amplitude 
de ces phénomènes reste très faible. 


DISPOSITIF EXPÉRIMENTAL UTILISÉ. — Les spec- 
tres ont été obtenus à l’aide d’un monochromateur 
Desvignes fonctionnant dans le visible, et d’un 
photomultiplicateur d’électrons A. Lallemand à 
vingt étages, débitant dans un galvanomètre Zer- 
nicke. L’étalonnage a été fait en utilisant une lampe 
Philips à filament de tungstène dont la température 
de couleur a été préalablement mesurée. 


IV. Résultats et discussion. — Nos résultats 
sont en accord avec l’existence des deux bandes À: 
et À, données par Henderson, Ranby et Hals- 
tead [5] et les positions de ces bandes dans les 
sulfures mixtes æ CdS (1 — x) ZnS(Ag) se raccor- 
dent progressivement aux bandes dans l’infra- 
rouge proche décrites par E. Grillot [4]. Nous don- 
nerons pour position de ces bandes : 


Am 1 Âm 2 
ZnS(Ag) blende & 500 À 4 950 À 
CdS(Ag) wurtzite 7160 À 7 800 À 


Les différences avec les valeurs données par les 
précédents auteurs peuvent être imputées aux 
erreurs d'expérience. 

En accord avec E. Grillot [13], nous considé- 
rerons ces bandes comme les homologues des bandes 
bleue et verte des ZnS(Cu), pour lesquelles un mo- 
dèle a été récemment proposé [14]. Ce modèle, en 
accord avec le schéma de niveaux de la figure 1, 
considère le niveau fondamental des centres comme 
située à une hauteur AE, de l’ordre d’un électron- 
volt au-dessus de la bande de valence. Le niveau 
excité du centre est situé à une distance AF, au- 
dessous de la bande de conductibilité ; AF% est 
faible, surtout pour le centre bleu. Pour le centre 
vert, ce niveau est légèrement plus profond et on 
peut le considérer comme un niveau donneur asso- 
cié au niveau accepteur fondamental [15]. Il faut 
admettre que le cuivre et l’or s’associent facilement 
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à un tel donneur, l’argent plus difficilement, du 
moins dans ZnS et dans les sulfures mixtes conte- 
nant une forte proportion de ZnS. Dans CdS(Ag) 
par contre, les bandes À, et à, retrouvent des inten- 
sités comparables, comme dans ZnS(Cu). 

Dans ZnS(Cu), la bande bleue est favorisée par 
les concentrations élevées en activateur, la bande 
verte par l’incorporation de donneurs (oxygène ou 
halogène) en concentration égale à celle de lacti- 
vateur. De la même façon, dans CdS(Ag), l’émission 
visible (queue de la bande X,) est défavorisée par 
l’incorporation du fondant NaCI ; van Gool [16] a 
montré que les phosphores avec des concentrations 
égales d’activateur et de coactivateur ne présentent. 
que l’émission infra-rouge ; il ne distingue pas deux 
pics dans cette émission car il a calciné ses sulfures 
à température élevée 850 ou 950 0C : ses résultats 
s’appliquent donc en réalité à la bande À. 

Enfin, du point de vue de la persistance de la 
phosphorescence, la bande bleue du cuivre est peu 
persistante, tandis que les centres verts sont liés 
aux pièges profonds. Les ZnS(Ag) qui présentent la 
bande bleue très prépondérante sont peu persis- 
tants. De ce point de vue aussi, 1l paraît done 
normal d’assimiler la bande fondamentale des 
ZnS(Ag) à la bande bleue des ZnS(Cu). 


PosITIONS ET ÉNERGIES 
DES DIFFÉRENTES BANDES D'ÉMISSION 
DUES AUX ACTIVATEURS Cu, Ag ET Au DANS ZnS ET GdS 


À hv Egap — hv 
cs ÂEËy + ÂEe 
ZnS(Cu) wurtzite & 450 À 2,79 eV 0,91 eV 
pê 5 230 À 2,37 eV 1,33 eV 
ZnS(Cu) blende & 600 À 2,70 eV 0,94 eV 
5 350 À 2,32 eV 1,32 eV 
ZnS(Ag) wurtzite 4 370 À 2,84 eV 0,86 eV 
k 760 À 2,61 eV 1,09 eV 
ZnS(Ag) blende & 500 À 2,76 eV 0,88 eV 
k 950 À 2,51 eV 14%ev 
ZnS(Au) Wurtzite & 700 À 2,64 eV 1,06 eV 
5 300 À 2,34 eV 1,36 eV 
ZnS(Au) blende & 800 À 2,59 eV 4,05 eV 
5 500 À 2,26 eV 1,38 eV 
CdS(Cu) wurtzite 8 200 À 1,51 eV 0,92 eV 
10200 À, 4,22 eV 1,21eV 
CdS(Ag) wurtzite 7160 À 1,73 eV 0,70 eV 
7 800 À 1,59 eV 0,84 eV 
CdS(Au) wurtzite [20] 8 000 À 14,55 eV 0,88 eV 
11 500 À 1,08 eV 1,35 eV 
BANDES A COURTE LONGUEUR D’ONDE : 
ZnS 3 900 À 3,18 eV 0,52 eV 
Cds 6 200 À 2,00 éV 0,43 eV 


BANDE A COURTE LONGUEUR D'ONDE DANS ZnS 
ET CdS. — Fonda et Klick ont montré, indépen- 
damment, que les CdS(Ag) à forte concentration en 
argent possèdent une émission rouge, avec pic vers 
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6 200 À [17], [18]. Cette bande n’apparaît qu’à 
basse température ; elle est souvent étudiée à l’air 
liquide, alors qu’elle est éteinte à température ordi- 
naire. 

‘Généralement cette émission est attribuée, soit à 
des ions Ag interstitiels, soit aux lacunes S. Mela- 
med [19] a donné des arguments en faveur de la 
seconde interprétation. Van Gool a suggéré [16] de 
la décrire par des centres du modèle de Lambe et 
Klick. 


La bande homologue dans les ZnS est située vers 
3 900 À 

Le tableau ci-dessous montre que l’énergie du 
photon lumineux émis est beaucoup plus voisine de 
la largeur de la bande interdite que pour les autres 
bandes d'émission dues aux activateurs Cu, Ag ou 
Au. En accord avec les auteurs cités ci-dessus, nous 
pensons donc qu’il ne faut pas considérer ‘cette 
bande comme homologue des précédentes. 

Manuscrit reçu le 4 mars 1961. 
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MOUVEMENTS INTERNES D'UNE PARTICULE RELATIVISTE 


Par M. KLEMAN, 


Théories Physiques, Institut Henri-Poincaré. 


Résumé. — Une particule relativiste sans champ extérieur, obéissant aux lois de conservation 
de l’impulsion et du moment cinétique, est d'ordinaire décrite par un tenseur antisymétrique Au, 


et sa dérivée covariante C4, ainsi que par une impulsion constante p,. Pour concrétiser le fait 
physique de la rotation interne, on introduit dans la théorie une deuxième particule ponctuelle, 
déduite de la première sans hypothèses supplémentaires : tout mouvement de point matériel peut 
être décrit comme un mouvement relatif de deux particules, lequel est analysé par un tenseur 
antisymétrique w4g, analogue au pseudo-vecteur rotation du corps solide de la mécanique clas- 
sique. On trouve ainsi une formule qui lie les vitesses unitaires des deux points, mais où og dépend 
de 6 paramètres arbitraires. Des considérations physiques simples de réciprocité entre ces deux 
points (réciprocité des impulsions et des vitesses) permettent de fixer ces degrés de liberté, et 
d'obtenir, dans le cas le plus général, une expression simple de la rotation. On passe ensuite à 
l’étude d’une classe particulière de mouvement de la particule relativiste, où la notion de réciprocité 
est étendue aux points eux-mêmes, la donnée de l’un permettant de passer à l’autre, et récipro- 
quement. On élargit ainsi le mouvement de Weyssenhoff, où les deux points réciproques sont 
confondus. 


Abstract. — À relativistic particle without external field, obeying to the laws of conservation 
of linear momentum and angular momentum, is usually described by a skew tensor A4, its cova- 


riant derivative A4, and a constant linear momentum p,. Inorder to materialize the physical 
fact of internal rotation, we introduce in the theory a second point particle without supplementary 
assumptions : the motion of a material point-particle can be described as the relative motion of 
two particles, relative motion which is summarized by a skew tensor w,g, of the same nature as 
the rotation pseudpvector of a rigid body in classical mechanics. So one finds a formula relating 
the unitary velocities of the two point-particles, but og depends still on 6 arbitrary parameters. 
By obvious physical considerations of reciprocity between these two point-particles (reciprocity 
of momenta and velocities), we can make these degrees of freedom disappear and the rotation 
reduce to a simple expression. We then restrict ourselves to the study of a particular class of 
motion for a relativistic particle, the idea of reciprocity being extended to point-particles them- 
selves. If one point-particle is given, the second is consequently found, and reciprocally. So one 
generalizes the so-called Weyssenhoff”s motion, in which the two reciprocal points have the same 
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space-time coordinates. 


Nous nous proposons de définir dans cet article 
des quantités physiques capables de décrire une 
rotation interne d’une particule relativiste infi- 
niment petite ; puis nous tenterons d’élargir ces 
définitions à une particule ayant une extension 
finie, à un fluide relativiste. On sait que ce dernier 
problème est difficile : la particule relativiste n’est 
pas un corps solide et on ne peut concevoir un 
nombre limité de degrés de liberté interne. Mais 
pour la particule ponctuelle cette difficulté dispa- 
raît et on peut doter cette dernière d’une structure 
interne très complexe sans que le problème de la 
disposition spatiale des tourbillons locaux appa- 
raisse. Cependant, comme l’ont montré D. Bohm et 
J. P. Vigier [1], un mouvement de rotation ayant 
une réalité physique est obligatoirement lié à la 
présence de deux points (centre de masse et centre 
de matière), même si ce mouvement peut se décrire 
ponctuellement. C’est de cette idée que nous parti- 
rons pour élargir le concept de rotation interne à 
un fluide relativiste, Nous montrerons d’abord que 


les descriptions déjà faites de rotations internes 
dans des particules ponctuelles sont en réalité des . 
mouvements à deux points, puis nous écrirons une 
condition physique liant ces deux points. Enfin, 
étendant notre définition bi-ponctuelle d’un mou- 
vement de rotation interne, nous supposerons qu'un 
fluide en équilibre est tel que la condition physique 
liant les deux points s’étend à toutes les « molé- 
cules » du fluide, prise deux à deux et suivies le 
long de leur trajet d'espace temps. 

Cet article reste entièrement classique, dans le 
cadre de la relativiste restreinte. Nous utilisons les 
coordonnées de Minkowski, avec € — 1, 4 — il. 


*k 
RAR 


I. Définition du rotateur relativiste. Son origine 
physique. — Il est classique de décrire un fluide 
relativiste par la donnée d’un tenseur énergie- 
impulsion symétrique 7°,,, obéissant aux équations 
du mouvement 9, 7,, — 0, Un tel tenseur permet 
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de définir un «centre de masse » animé d’un mouve- 


ment rectiligne uniforme, et une impulsion p} 


constante du mouvement. On pose 
à = — if Tu dv (1) 


Mous = — à [Ur — Ya) Tas— (re — 9) Pol do (2 
x, est la variable d'intégration, Y, est un point de 
l’espace-temps où nous calculons AC,,,(Y). Les inté- 
grations précédentes sont faites sur un hyperplan 
t — constante, mais 1l est facile de démontrer que 
les quantités p, et AC,, ne dépendent pas de la 
surface d'intégration, à cause de la loi de conser. 
vation 2, 7,, = 0 (cf. Müller [2]). La même loi de 
conservation conduit à 


En 0 k (3) 
Jeu = Pu Ur Py Pu (4) 
où À désigne la dérivée le long de la ligne de courant 


suivie par Ÿ,, À = dAfJdt et v, — dY,/dt, avec 
la relation supplémentaire vu, = — 1. 

On voit qu'à chaque point ÿ, d’un fluide rela- 
tiviste on peut associer une vitesse ’,, une impul- 
sion p,, la même pour tous les points ; un tenseur 
antisymétrique Mu généralisation du mouvement 
cinétique de la mécanique classique. Mais on peut 
encore associer à chaque point Ÿ, d’autres quan- 
tités importantes, [3] qui sont 


Av py = mm Ry Vs Du Dre (5) 
gs Ps mSy A RES (1/2) Euvxr Aux (6) 
Aus D = My ER Re (7) 
ous Op Sur (8) 
R, et S, sont directement liés au centre de masse 


que nous voulons définir. À, désigne le vecteur 
joignant le point Y, au lieu de ce centre de masse, 
l’axe central, qui est une droite parallèle à l’impul- 
sion Py- S1 C, désigne le centre de masse, on montre 
que l’on'a GC, =:Y, —-h;; Ilsufit d'écrire en 
effet : 


Ayy Py = — ip [ (ru Tva — 29 Ty) dv 


im [ Ty an) Fin vs) f Ta do. (9) 


Soit en utilisant (1) 
ACyy Py = — ipy [tu Tya — Xy Tya) A0 


Tes Yu (Po Pa) + Pu(Po Ya): (10) 


Cette équation s’interprète facilement : appe- 
lons I, l’'hyperplan perpendiculaire à p4 et passant 
par Y.. Dans cet hyperplan Pa Ya est une cons- 
tante, Il en est de même de la quantité sous le 
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signe À qui ne fait pas intervenir le point Y,. 


Posons maintenant À, — Y, = — C,. il vient, 
en utilisant (10) 
Ry Lei Yu RATES. ipu [ (ru T'ya — 19 Tua) dv 

+.pu(Pa Ya) = — Cu (11) 


C, est done le même pour tous les points Y, de 
1K hyperplan I,. Lorsque IT,se déplace parallèlement 
à lui-même, on voit facilement que €, décrit une 
droite (l'axe central) qui est le lieu du centre de masse 
de D. Bohmet J. P. Vigier [1]. R, apparaît comme 
le vecteur d’origine le centre de masse, d'extrémité | 
le point courant, dans cet hyperplan. 

On montre de plus que S,,au même titre que ps, 
est une constante du mouvement [3]. 

Si 5, et R,,de par leur définition, sont des quan- 
tités globales, qui décrivent une propriété unique 
du fluide, la rotation autour de l’axe central, les 
quantités 4, et s, sont en revanche des quantités 
locales liées en particulier à la vitesse unitaire. On 
peut d’ailleurs lier 4, ets, locaux à leurs homologues 
globaux sous des formules qui ressemblent aux 
théorèmes de Koenig. On utilise pour cela l’expres- 


sion de À, 


11e Ru = MES Pv = m'°[vy — (mm?) ru] (12) 


d’où 


AGuv PE mu — Ru) (13) 


Me À = Su — (m/m') Sue (14) 


On peut interpréter s, comme le moment ciné- 
tique pr en un point quelconque de la masse 


fluide, =>; re comme le moment cinétique pris aux 
centre “ la SRE Ho, Rest en effet, dans le 
repère Il,, le moment cinétique de rotation du 
point YŸ, tournant autour de C,, avec la masse m’, 
Suivant Takabayasi [41, qui nomme impulsion la 
quantité p,, nous appellerons m’ la masse au repos 
apparente, m la masse au repos vraie, 

Il est possible de définir un autre point que le 
centre de masse C,, qui donne une idée du mouve- 
ment tourbillonnaire de la masse fluide. C’est le 
centre de densité de matière (1), à partir de la 
donnée d’un flux de courant 7, conservatif, 
du Ju — 0. Mais le centre de densité de matière n est 
défini que dans un repère bien particulier, ce n’est 
pas une donnée covariante. 


ve! (av) — far tan. 


Cette définition est analogue à celle du centre de 


Masse 
Cu | [8 do | = fr NW dD 


(15) 


(16) 
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. Les intégrations sont faites dans le repère IT,. 

Mais alors qu’on peut donner par (11) une défi- 
mition covariante du centre de masse (on démontre 
facilement l’équivalence entre les équations (11) et 
(16)), il n’est pas possible de le faire pour le centre 
de densité de matière. Son sens physique est cepen- 
dant précis. Le mouvement moyen de la masse 
fluide ne peut-être décrit par un seul point, car 
l'observateur placé dans le repère propre de ce 
point n’aurait aucun moyen de connaître son mou- 
vement de rotation. Au contraire, s’il en voit un 
autre, il sera sensible à ses variationx angulaires et 
en distance. Un cas bien étudié jusqu'alors a été le 
mouvement de Weyssenhoff, qui correspond à la 
condition supplémentaire 4, — 0. Dans ce cas le 
centre de masse et le centre matière sont animés 
lun autour de l’autre du même mouvement de 
rotation pure, avec même vitesse angulaire. 

Nous en arrivons maintenant à la définition d’un 
rotateur relativiste, qui découle tout naturellement 
des données précédentes. C’est une particule ponc- 
tuelle animée d’une impulsion constante (conser- 
vation de l’énergie impulsion), et d’un moment 
cinétique AC, tel que AC, = p, v,— p, v, (conser- 
vation du moment cinétique). Cette particule est 
encore définie soit par la donnée de p4, R,, Su Dur 
soit par la donnée de p,,,1,,s,, v,. Donc par quatre 
quantités vectorielles, dont une longueur 4, ou À, : 
on a ainsi la donnée de deux points, et nous allons 
montrer que le mouvement de l’un par rapport à 
l’autre est caractéristique des propriétés dites in- 
ternes du premier point, soit p, et s. 


II. Description du mouvement relatif de deux 
points. — Avant de passer à l’étude que nous avons 
définie dans le $ I, nous allons établir une formule 
permettant de lier entre eux les mouvements de 
deux points sur leurs lignes de courant respectives. 
Nous allons pour celà introduire un tenseur anti- 
symétrique @,e, qui est la généralisation relativiste 
des vitesses angulaires, comme nous le verrons plus 
Join. 

Nous attachons à chaque position de la trajec- 
toire du point matériel envisagé un ensemble de 
quatre vecteurs orthonormés (vierbein) b® où (E) 
désigne l’un des quatre vecteurs (£ — 1, 2, 3, 4), 
et p sa pieme projection sur les axes du laboratoire. 
Les 16 quantités b® ainsi définies sont liées par 
10 relations d’orthonormalisation. 

De pe = dyy D bu = den. (1) 

On a donc affaire à 6 paramètres arbitraires fonc- 
tion de +, le temps propre du point matériel par- 
courant sa ligne de courant. Cette variation du 
repère mobile (que par la suite nous confondrons 
avec le repère propre du point en alignant b£et v,,) 
peut-être caractérisée par les différentes vitesses 
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cinématiques E (b®)). Appelons «,, la projection 


de la vitesse du y ième axe du laboratoire [# sur le 
v ième, vitesse prise évidemment par rapport au 
repère mobile. On a 

1 = 6 b® (2) 


d’où en dérivant par rapport au temps propre 
des D®, et en gardant les axes b#) fixes 


le Les Le AI 9 
d'et u (3) 
ou 
Ouy = ble) b® pt) PAUL) (4) 
ô re es 
Ouy = à) LE => (ie 1E) kr, be) be) (5) 


t 


à cause des relations d’orthonormalisation 
(E) DE) 
pie pe — du 
On pourrait tout aussi bien utiliser les vitesses 


angulaires du repère mobile par rapport au repère 
du laboratoire ; on trouverait alors 


DHErEr PRpÈTe (6) 


Considérons maintenant un autre point maté- 
riel X décrivant une autre trajectoire que la précé- 
dente, avec le temps propre 0. Nous allons décrire 
le mouvement de X dans le repère b®. », désignera 
la vitesse unitaire de M, w, celle de X. On se donne 
en plus une fonction arbitraire + — /(0) telle qu'à 
un repère b® correspondent deux points bien dé- 
finis M et X sur leur ligne de courant respective. 
Ces lignes de courant étant supposées des trajec- 


“toires réelles, il suffira par exemple de prendre M 


et X au même temps physique £. 
Posons MX=K=K, Fe On a 


X = M + MX. (?) 


Dérivant sous cette forme par rapport au temps 
propre +, on fait apparaître les vitesse unitaires v, 
et w,. I vient 

LATE C d , 
de —= di” Ie = —— + — MX (8) 
cM 


MX = Ky It Ke DS 142. © (9) 


—= Uy Iv 
Nous notons par le signe prime les projections 
sur les axes mobiles, et nous utilisons les formules 
de transformation du repère mobile au repère fixe, 
et inversement, qui sont : 
rt E y 
Ke = DE Ky 


Ras Ke be (10) 


d’où le calcul suivant : 


dMX d pen ve) Ge » J() SE 
= CA LP) 6 = (ui KE + &,0Ÿ ie) Fe 


— (2 É Ke + Ov k,) Ie. (11) 


4 


438 
Comparant (8), (9) et (11) et revenant à la nota- 
tion tensorielle, il vient : 


do 
dr 


d 


We = Vy, + œuv Æy + RE 


(12) 

Tels que nous avons fixés les w,,, c’est-à-dire 
absolument indépendants du contexte physique, 
cette formule n’est pas une définition de la rotation 
instantanée d’uneligne par rapport à l’autre. Cela ne 
sera que s’il est possible d’imposer à w,,, six condi- 
tions physiques simples. Pour le corps solide de la 
mécanique classique, l’on sait que l’on a une for- 
mule analogue à (12), mais où le dernier terme est 
nul. Cherchons dès maintenant la signification de 
l’annulation de ce dernier terme. Partons de l’ex- 
pression K; — DE K, et dérivons par rapport à 7, 
puis saturons par TER il vient : 


ke DO = Ky + op Kv. (13) 
La condition cherchée est donc 
ky: — Qu ÆY. (14) 


Elle signifie que la distance X,, K,, est constante. 
Réciproquement, si cette condition est réalisée 
(X, K, = 0), il existe w,, tel que le dernier terme 
de (12) s’annulle. Remarquons simplement pour le 
moment que cette condition apparente l’ensemble 
des deux points à un corps solide. 


III. Étude de w,, dans le cas général. — Le 
point courant M est défini par les quantités p,,,s,, 
tv. Nous y adjoignons un second rotateur rela- 
tiviste simple, le point X, de vitesse w,,, tel que 


A, = Y, —à (1) 
r sera le temps propre de M(Y,), 0 celui de X. 
En utilisant IT (12) il vient : 


d8 de 
FEU = y — Opy L — RE) te LE. (2) 


La dérivation de (1) donne aussi 


d0 


= My = Vy — ty (3) 
qui entraine, en saturant par v,,p, et en élevant 
au carré, les formules suivantes, que nous aurons à 


utiliser 


(JO dr Em ENS me -L Ve ta) 
10 

—— Pa We = M + le Pa = æ(d0 dr) (5) 
US EC (6) 
— (dO dr)? = 21 96 be Etote (7) 

Nous poserons en outre 
1 
er A LE — Ag. (8) 


dr 
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Introduisons la distance de X à l’axe central, ce 

sera, d’après (1) 
1B = Py — Projulta) (9) 


où Proj,(t) représente la projection de #, sur un 
plan LE nie di à p,. On a,en posant pe = MU 
(HU = — 


Proju (te) = Ouv + Wu Us) ty = tu + (ua to) tu MO) 


AD = —= Py ES lu << (ue la)Uy. (11) ‘ 

Si on définit pour le point adjoint X un tenseur » 
antisymétrique M,£ qui rapporte ce point au même . 
axe central que M, on aura donc 


Mis = Dr 62 — pp 1 a Fe ee en Se Pre (12) 
Soit en utilisant la valeur de Fe (tenseur anti- 


symétrique de M rapporté au même axe central) 


Mg = Mg — (Pa 1@ — PR ta). (13) 


Par la suite, nous considérons que, physi- 
quement, p, et5, représentent les propriétés com- 
munes aux points X et M. Ainsi nous appellerons 
masses au repos vraies les quantités m = — p, ve, 
u — — p,w,. Ce seront les masses dont nous « 
affecterons les points X et M lorsque nous les 
considérons comme deux points conjugués au sens « 
de Yukawa (théorie de la bilocalité). 

Les 16 paramètres b® sont liés par 10 rela- 
tions bE 5 = 5, On peut donc fixer 6 conditions 
qui doivent contribuer à donner un sens physique 


à (2). On est d’abord conduit à aligner D avec la 4 


vitesse propre de Y. 
D'où 3 conditions 


CE = — y (14) 


Ov tb &,, Ont alors la forme suivante, si on pose 
On Vi = Py 


Ouv = Un Vy — Du Uy — LEpvex Ve OX (15) 
Ouv = y Éy = y Pos LEuvxX Ve UN (16) 
Nous sommes conduits, en utilisant des résultats … 
de mécanique classique, à considérer que Fénergie 
de rotation de la masse fluide est le produit d’un 
terme de vitesse angulaire par un terme de moment 


CRE i PENSE 
cinétique 1 — 5 Panies Ici nous prendrons pour 


DE 2 . 
1 5 Ag (ou Mg), pour © : ©,8. Nous intro- 


durons donc deux énergies de rotation 


[g AI il à 
1 1 = à ACp Cap (17) 
Lg AI 1 / 

1 = L Mg Oop. 


Nous avons déjà étudié [5] des termes d'énergie 


: 
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analogues, dans le cas du mouvement de Weys- 
senhoff. Nous avions alors trouvé T — m2. Ici 
deux raisons vont nous conduire à adopter la même 
valeur pour T,;. On remarque d’abord que X a 
exactement un mouvement à la Weyssenhoff dans 
le repère propre de A (il est facile de montrer 
que M,g vg — 0). D'une façon générale on peut 
done associer à tout point un point qui ait autour 
du premier un mouvement satisfaisant aux équa- 
tions de Weyssenhoff. D’autre part nous supposons 
que les points X et Y sont conjugués, la vitesse de 
lun étant l'impulsion de l’autre. 7, représente alors 
énergie du point Y, de vitesse v,, d’impulsion 
UN, (5) , Ta l'énergie de X, de vitesse w, d0/dr, 
d'impulsion mv,, d’où deux conditions supplé- 
mentaires (nous ne ferons pas intervenir le fac- 
teur v, w, d0/dr, pour rendre les expressions plus 
simples). 


T;, = (1/2) ud0/dr Ta Yi PA (18) 
Pour calculer 7; et 7,, nous utilisons l'expression 
de AC,£ en fonction de p,, D, Sy. 
l Huy — LEpvxA Sx VX + ne Du — {y y) (19) 
ACyy = Sy Vu — Su Vy + LMEnyxà le VX. 


Il vient alors, pour T'; et 7, 


| 


: ‘| . : - 
T; = % ÂACp Wap — 2 [la Va + 9x Sa] 
(20) 
L4 n : 
me : Mg Ou = 5 la Va + Ex SA] + 3 uv Pv lu. 


D'où la valeur de +; s; en fonction des données ; 
ainsi que de ©,, p, ty 

en 81 —= u(d0 dr) — mt Van + (Pa INUx) É (21) 
(22) 


Ouy tv Pa = Mm— U{(d0[dr) = Pa de 


Reprenons l’expression (2) pour voir ce qu’en- 


traînent pour la quantité À, — — bé + t les 


T 
hypothèses précédentes. Il suffit de saturer (2) res- 
pectivement par v,, Py, l, Ce qui donne : : 
(d0 Jdr) we Ve + 1 + op tv Vu = 
(d0[dr) we Px — Va Pa + Ouv tv Pu = Au Pu 
(d0 dr) wy to Æ 


Au Vu 
(23) 
Au lu. 


Soit encore, en utilisant les équations (3), (4), 


(15), (22) 
porter 0 (4) 


Ag lu = — lg la 
On est alors tout naturellement conduit à poser 
comme 6€ condition 
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condition qui conduit à étendre la symétrie imposée 
aux deux points X et Y par les équations (20). À, 
est un vecteur du genre espace qui est perpendi- 
culaire à toutes les vitesses données, p,, 0, We. 

La dernière relation permet de trouver l’expres- 


sion de oo, É,4, 


(d0 [d=) Wu tr cr Ur lu + Ouv % ly = 0 


Opy y lu = Lg lo. 


Nous allons commencer par calculer le vecteur 
Ouyly — Q,, qui est entièrement défini par les 
relations (23) et (27), soit : 


Qy Va = ta Va 
Qu Pa = le Pa 
(ox le —= la lœ 
OPA 0: 


Les trois premières de ces équations permettent 
d'écrire Q, sous la forme d’un vecteur égal à £, 
œ 
plus un vecteur perpendiculaire à D, Pa, le SOit 


@ 2 Lœ + im Q Eve Vy Py a (29) 


La dernière relation Q,+t,— 0 permettant de 
fixer l’unique inconnue Q 
QUE 0e ta té im O'euive DhPvi5 ta OU] 
mi, se calcule facilement à partir des expressions 
I (4) et I (7) 


Mo = — te + MUplIx dx + 1) — Po 


(31) 


- LEpa/y'v Da’ Su’ y. 


Si on porte cette expression dans (30), seul sub- 
siste le dernier terme de (31), à cause de l’ortho- 


gonalité de (Q, — t) QUO AD AR OTÉE 


0 —= lo la + Q Epatiy Epx'u'v” Lœ Üx! Un Su’ Py Uy’. (32) 


On symbolise ssuv: ours Par de, qui est égal 
à O si les groupes auv, «'4"v’ sont globalement diffé- 
rents, à + 1 si on passe de l’un à l’autre de ces 
groupes par une permutation paire ou impaire. On 
aura donc à considérer, avec les signes adéquats, 
les permutations 


uv — av” e (ta Da) (Cu Su) (Ps D) = 0 
œuv = L'a'v” à (la Sa) (Eu du) (Pv Ds) = 0 
Œuv = VX + (x Sæ) (Un bu) (Py ds) 

au = va Es (lo Va) (bu Su) (Py 0v) = 0 
uv + v'ax’u” 2e (la La) (Vu Du) (pv Sv) = 0 


auv = «VU Se (x Da) (Du Vu) (pv Sy). 
Soit : 


Q | (to ve) (Po si) — (Te A) (Pa Da) Lis tulu. (33) 
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Pour calculer 9, nous utiliserons la valeur de Q; 
qui le contient, soit 


Port — VA (D La) — LExapy VE Pyto (34) 


et nous égalons cette quantité à (29). Il suffit de 
multiplier par ie, Sy V, pour faire apparaître 9: 
que l’on cherche. Én effet l'on a, par un calcul 
semblable au précédent 


EhaBy Sauve VB y la Su Vy = lof Pa Sa) — (la Sa) Ph (35) 
où 4 S4 est donnée par (21). 
L’équation qui donne directement ; s'écrit alors 
LEpvrp Sie Uy On — Exapy Exuve VB Pr la Su Vu 
= lp(Pa Sa) — Ppllo Sa). (36) 
Soit, en exprimant Lejuvo Su V, en fonction de A(;,, 
Sas Va 


10e — Mo Va + Mia Ve) (x = lo Pa Sa) — Pplla Sa). (37) 


Reste à utiliser l'expression (21) de w,s,, et on 
obtient une relation simple entre @, et Q; 
(ta Sa) Pp — Abex 1 + u(dO dr) t 


Le calcul de AC,1 Q; en fonction des données ini- 
tiales est un peu plus laborieux. On a 


(38) 


AC,x Qù = Àrta D, — Milo VA + LEpauv Su Vu) 
—,rn|Da ta) to + iQ papy Su Vv (39) 
Epauv Su Vy à = Eau n Su Uy 
+ im Q eus Exaëy Vo PB - Su Uv. (40) 


Le second terme du second membre de (40) se 
calcule encore de la même façon que plus haut, il 
vient 


im (À Euv Exapy Va PR # Su Vy 
— 1m O[to(Pa Sa) — Pelte Sx) + VolPa Sa) (tx 0x) + mix sx] 
(41) 
LEpruv Su Vy Où — Ai à — Le { ml to) — MO {Pa Sax) (to Vo) 
— mm? À 14 So — m Op £a) — Polte Sa)] + mto(vx in). 
(42) 


Pour simplifier cette équation, on peut utiliser la 
valeur de (, qui s’écrit encore suivant (44), compte 
tenu de ce que, en saturant (31) par s,, on obtient 
(43). 


la £a) + lle So) = — Po Sa (43) 

Oflta Vx) (Pa 8x) — Pa Sa — Mltx Sa)) = — tata (44) 
d’où 

MD À lo ta — Q[(Pa Sa) (la Va) + Ma Sa) } = MPa Sa) Ve. 


(45) 

L’équation (39) se simplifie alors beaucoup, et 
il reste 

ACL Où = Ah D — MPa S Sa) V9 — mt sise Sa) — Polta Sa]. 

(45) 


Nous allons maintenant tirer de cette équation a | 
une expression plus simple que (33), pour Q 


ga v — 0 entraîne, en repartant de (38), et en el 


prenant pour A6, Q la valeur ci-dessus 


se) + MA da) (Pa Sa) — Mia Sul. 
(46) 


0 — Mmle ta) + M(Po 


Utilisons encore une fois (43) pour simplifier Je 
facteur de Q il vient : 


0 — ma la) + MPa Sa) 


+ mOi Va) (Pa Sa) + la Sœx) + Pa Sal. 
Si on utilise Pexpression (33) de Q il vient : 
DEEE 


(48) | 


IV. Cas de la bi-ponctualité réciproque. — Dans 
le paragraphe précédent nous avons utilisé pour 
fixer les valeurs de l’énergie, la notion de réciprocité 
entre-les deux points X et M. En fait, nous n’avons | 
guère poussé cette question, supposant simplement 
qu’elle se limitait à l’énergie (à l’impulsion). Mais, 
tels que nous avons les avons fixé dans le para- 
graphe précédent, si la donnée de M donne X, 
inversement, la donnée de X, considéré lui aussi 
comme un rotateur relativiste (p,, w,, M,,w,, 


M,,w,), ne donne pas M. Il suffit d’ailleurs, pour « 
s’en assurer, de remarque que f,, le vecteur jo1- 
gnant X à MW, n’est pas perpendiculaire à w,, donc 
n’est pas parallèle à AZ,, w,. 

Une bonne définition du rotateur relativiste 


€pur », et qui mène à l’étude d’une classe de mouve- « 


ments, qui généralisent immédiatement le mouve- = 
ment de Weyssenhoff, est donc de considérer les 
deux points comme réciproques en tant que rota- 
teurs relativistes (p} — pé, si— s£, 1} — — 1%) ce 
qui, joint à la condition de réciprocité énergétique 
comme nous l’avons introduite plus haut, a un sens 
physique bien clair. Le mouvement apparait en 
effet comme stable, puisque l’itération du procédé 
qui fait passer de M en X ramène en M. On pour- 
rait peut-être étudier des itérations particulières qui 
mèneraient à un point d’accumulation, ou redon- 
neraient cycliquement le point de départ. Mais dans 
ce cas on imagine mal quel serait le processus d’ité- 
ration de l’énergie. Nous signalons que l’un des cas 
envisagé dans ce paragraphe est justement le mou- 
vement que décrit T. Takabayasi [4] pour le lagran- 
gien à variables spinorielles : 


RE Le AN Ci: 
+R +: En 1} (0 


Si si — si, 1% — — t}, comme nous en faisons 
maintenant 1 hypothèse, ‘alors on à 


0. 


We SES We 


(47) 4 
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Ces deux conditions, portées dans III (3), donnent 
tu Sa = 0 (2) 
la ta = 0. (3) 


On en déduit d’abord que la distance 


AM = VU te 


est une constante du mouvement. D’autre part, il 
est facile de montrer que, dans le cas général, on a 


Ste = 0 (4) 
à cause de 


Su — (Sx da) Vu + Êm Epvxn tx Un v. (5) 


On en conclut que s, tiest aussi une constante du 
mouvement. Nous avions introduit le vecteur À, 
au paragraphe précédent 


(Av — A, pu = Au =-0) 
Nous avions en outre At, — —t, 1. lei cette 


quantité est nulle, et on en déduit qu’il existe une 
relation linéaire entre p,, te, Ve, Wa. 


; Pu —= Aty + Bvy + Cwy. (6) 

En saturant d’abord par v, et w,, on calcule B 
et C d’où 

Pu = Alu + n 


TT (mwy — uvy). (7) 
Si maintenant on multiplie par s,, on en tire la 
valeur de À, soit Pa Sa, Mais on a aussi, en multi- 
pliant HI (31)pars, 
Pa Sa = —— la Sa (8) 
d’où la valeur de p, 


mt LL 


— — jhty < Vu + : 
Pu CAIN ET 7 LACS RERS TEEN 


Wu 
ep Cf 
D Te 


Remarquons qu’en saturant (9) et IIT (31) par 4, 
au lieu de $,, on obtient 


Fate = — Maty à partir de (9) (10) 


à partir de III (31) 
(11) 


Paula — —— la le + LEpvxn Sx vx Dy lu 


Le second terme du second membre de (11) est 
Su $9, à cause de (5), d’où 


Sa Se — 0, équation qui fait pendant à tt = 0,» (19) 
Reprenons l'identité IIT (33), où Q — — 1, c’est 
(13) 


(£œ Va) Pa SA = (la Sa) PA x, 


MOUVEMENTS INTERNES D'UNE PARTICULE RELATIVISTE AA 


que nous comparons à 
PA SA = (lo Sa) Px UN Fe Me Sy. (14) 


Ici comme £, s, = 0, d’après (2), il vient, suivant 
les cas 


A) Si 
Pa da = — M # 0 et Pa Se 0 
alors 
ta De — Le Un = A. 
B) Si 
crie ROUE AU EL 


on devra avoir, soit r — 0, soit 4,5, — 0. Nous 
allons donner rapidement quelques indications sur 
chacun de ces types de mouvement. 


A) Comme v,t, — — 1, si nous portons cette 
valeur dans IIT (3) il vient 


(40 /d+) wa Ve = 0. (15) 


Ou bien d0/dr — 0, et alors X est un point fixe 
de repère du laboratoire, et M tourne autour de ce 
point fixe à une distance constante ; ou bien 
Ve Wa — 0 et les deux points tournent l’un autour 
de l’autre de telle sorte que M atteint la vitesse de 
la lumière dans le repère propre de X, et récipro- 
quement. On a alors 


Pu = — y + My + UWy. 


Dans le cas où d@/dr — 0, les équations du 
paragraphe précédent donnent lieu à des calculs 
simples, et on trouve 


Uno Su == Phiy TROIE Ay — 


(16) 


Opv [M —" Une 


B) Ici deux cas se présentent encore, selon que 
la Sa égale zéro ou diffère de zéro. Dans le premier 
cas, on en tire d’abord p,s, — 0. Ces valeurs, 
portées dans III (45), conduisent à AC, 41 = 0: 

Mais plus intéressant est le cas où t, 8, — 0, ce 
qui conduit tout naturellement à poser 74 — 0. 
Nous allons montrer que t, et s, sont alors paral- 
lèles. 

Remarquons tout d’abord que r# — 0 entraine 


Pa Sx = Ù Pa lœ = 0 (17) 


En) Pa e = 0, 


Si nous nous plaçons alors dans le repère p; = 0, 
nous voyons que nous avons affaire à quatre vec- 


teurs Se So les la Qui se trouvent entiérement dans 
un hyperplan à trois dimensions. Comme, d’après 
les équations du début du paragraphe, Sa et le Sont 
tous deux perpendiculaires à &£, et {,, on en conclut 


que s, et {, sont colinéaires. D’où 
Su — Aty. (18) 
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À est une constante du mouvement. En multi- 
pliant (18) par s,, et en se souvenant que 5, s, et do la 
sont des constantes du mouvement, on en déduit 
la même propriété pour À. On peut donc écrire 


Si nous utilisons maintenant les valeurs de &, et 


{, données plus haut par IIT (31) et IV (5), et que 
nous les portons dans (11), il vient la valeur de À, 
soit : 


Sy = mi 
& He (20) 
Su — Ml 
Comme 7% — 0 la valeur de p,, dans le cas envi- 
sagé est une combinaison de v, et w,, c’est 


1 
Pu AREA (mDy — uwy) 
JTE y — uvy) d'après (9) (2) 
ES (ve PA mwy Up eh es 


A, étant un vecteur perpendiculaire à Par Var Wp 
A, est 101 nul, puisqu’ii existe une combinaison li- 
néare entre ces trois vecteurs. On peut donc écrire 


(d0 dr) Wu —= Vu ee CONS ly. (22) 


Le cas ci-dessus étudié est justement le mouve- 
ment de Takabayasi [4], où les hypothèses de base 
étaient mn — p, s, — 0. Nous avons vu ici comment, 
grâce à la notion de réciprocité, ces hypothèses 
s’introduisent d’elles-mêmes. 
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Nous avons done vu que la notion physique de 
mouvement bi-ponctuel fait apparaître commo- 
dément un certain nombre de mouvements simples. 
Il fallait évidemment que la donnée du premier 
rotateur relativiste permette de définir le second 
sans ambiguité, et c’est pourquoi nous avons sup- 
posé que le tenseur J{0,, à lui seul permettait cette 
définition, qu’il représentait donc à lui seul le mou- 
vement interne d’une particule. Nous avons alors 
introduit l’autre extrémité du vecteur {,, qui joue 
un rôle analogue au premier rotateur relativiste 
(la particule donnée). 

Comme nous avons vu que A, décrit la rotation 
locale, nous pouvons supposer un champ de JC, 
c’est-à-dire un fluide de particules, et adjoindre à ce 
fluide l’ensemble des particules réciproques. On 
pourra alors, par exemple, écrire que le mouvement 
global du fluide est stable quand l’ensemble des 
lignes de courant réciproques se confond avec l’en- 


. semble des lignes de courant directes. Ce serait une 


sorte de condition de Weyssenhoff appliquée à un 
fluide. 

On peut enfin espérer, dans le cadre des théories 
de particules qu’on élabore à l’Institut Henri- 
Poincaré, que ces quelques mouvements classiques 
aideront à définir des critères de classification ou 
sont des descriptions classiques de mouvements in- 
ternes de particules de la classification. 


Manuscrit reçu le 19 novembre 1960. 
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ÉLARGISSEMENT ET DÉPLACEMENT DES RAIES DE RÉSONANCE MAGNÉTIQUE 
CAUSÉS PAR UNE EXCITATION OPTIQUE * 


Par J. P. BARRAT, 


Laboratoire d’Optique, Faculté des Sciences, Caen. 


et GC. COHEN-TANNOUDJI, 
Laboratoire de Physique de l'E. N. $., 24 ,rue Lhomond, Paris (5€). 


Résumé. — L'effet global des processus d’absorption et d'émission de photons de résonance 
optique sur l’état fondamental d’un ensemble d’atomes est étudié à partir d'équations d’évolution 
établies précédemment. On met ces équations sous une forme plus simple grâce à une « approxi- 
mation séculaire », toujours valable dans les cas pratiques. On tient compte également de l’effet 
éventuel d’un champ de radiofréquence. Certains résultats sont établis, tels que : la description 
du pompage optique comme un processus de relaxation dans l’état fondamental : la possibilité 
d’une conservation partielle de la cohérence au cours du cycle de pompage ; la prévision de certains 
déplacements des raies de résonance magnétique sous influence de l° excitation 0 optique ; l’expres- 
sion ve signaux de détection optique à partir des éléments de la matrice densité dans l° état fonda- 
mental. 


Abstract. — The total effect of absorption and emission of resonance radiation on the ground 
state of an ensemble of atoms is derived from equations established in an earlier paper. These 
equations are put in a simpler form by means of a ‘‘ secular approximation ?” which is satisfied 
in all practical cases. The effect of a radiofrequency field is taken into account. Some of the 
results obtained are : description of optical pumping as a relaxation process in the ground state ; 
possibility of a partial conservation of coherence during the pumping cycle ; prediction of some 
shifts of the magnetic resonance lines under the influence of the optical excitation ; expression 
for the optically detected signals as a function of the elements of the density matrix in the ground 
state. 


IV. ÉVOLUTION GLOBALE 
DE L'ÉTAT FONDAMENTAL. 


constantes de même ordre de grandeur multipliées 
par des exponentielles eff (r entier). 
2. RÉSOLUTION DES ÉQUATIONS (IV-1) PAR 


A. Équations d’évolution. — 1. ÉTABLISSEMENT 
DES ÉQUATIONS. — Pour obtenir l’évolution globale 
de la matrice densité 9 représentant l’ensemble des 
atomes dans l’état fondamental, nous ajoutons les 
effets de l'excitation optique, décrits par (II, 13), et 
ceux de la retombée, décrits par (IIT, 10), obtenant : 


d 1 ] &) ài( RNA 
1 u—{4”)0 4 
MP Éa + AE" 2 Ag’ Pu’u’ € f 
2 3 


L MT iu/—p)o,t 
> Cr: Le ia£) D Ayry Puur € EAU 
2Llp u” 


sp. 
nn 


Du" 


Quru/" 
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Bin Lun 


il(u—1v)—(ù "pot 
e = 
i(u — x) we — 


(RE — u") of] 


(IV-1) 


I" + 


Le système d'équations différentielles linéaires 
auquel obéissent les p,4 à pour coeflicients soit des 


1 
constantes d'ordre de grandeur 7 °u AE, soit des 
? 


(*) Cette étude fait suite à un premier article paru dans 
le précédent numéro du Journal de Physique [1]. 


APPROXIMATIONS SUCCESSIVES. APPROXIMATION SÉ- 
CULAIRE. — Nous montrons en appendice I que le 
système obtenu (I1V-2) en ne gardant que les coeffi- 
clients constants du second membre des équations 
(IV-1), 


d L AE" 
qu = [(o7g + EE) Aux 
— iAE") Al pus 
“te Es =? Ay'u| Puu 
UE Du AV 
: | rt  (IV-2) 
j Th ul! LUN Oe Tor) ia 


DATANT 


a pour solution une solution approchée de (IV-L), 
lorsque les condition (IV-5) 


of [AE > 
Of Th D { 


sont satisfaites. L’erreur commise sur la solution 
si l’on remplace (IV-1) par (IV-2) est de l’ordre 
de 1/or T,. Cette approximation, dite « approxi- 
mation séculaire », sera faite dans tout ce qui suit. 
Physiquement, les termes exponentiels ef#;t oseil- 
lent trop rapidement (à la fréquence de Larmor 


(IV-3) 


AAA 


@, [2x de l’état fondamental) pour qu’ils puissent 
avoir un effet appréciable sur l’évolution de p qui 
n’est notable que pour des temps de l’ordre de 
AJAE", T, > 2r/or. 


3. JUSTIFICATION DES APPROXIMATIONS FAITES. 
— Nous venons de voir que p, varie avec une 
vitesse caractérisée par les temps 7, et 1/AZ. 

Or nous avons supposé dans le calcul de l’expres- 
sion (II-4) que ps variait très lentement sur des 
intervalles de temps de l’ordre de 1 /A. Cette hypo- 
thèse est justifiée puisque d’après (1-1) et ([V-3) 


1IAE 1 lot LT» VIAE. 


D'autre part, dans le calcul de lexpression 
(III, 10) nous avons admis que p,, variait lente- 
ment devant la durée de vie + — 1/[°' de l’état 
excité. Nous devons donc compléter les conditions 
de validité (I1V-3) par 


FJAE’" 51 
LT 


Les conditions ([V-3) et (IV-4) expriment physi- 
quement que le déplacement AE’ et la largeur 1 IT 
d’origine optique des raies de résonance magné- 
tique sont petits devant la fréquence «7/27 de ces 
raies (fréquence de Larmor de l’état fondamental), 
et devant la largeur naturelle du niveau excité. 
Elles portent en! fait sur la fréquence «,/2r et sur 
le flux lumineux de pompage. Elles sont toujours 
remplies dans les expériences usuelles. 

Il faut y ajouter la condition (1-1) À > «e, y, 
VF. On pourrait s’affranchir de la restriction À > &, 
w, en compliquant l'écriture des équations. 


(IV-4) 


4, MATRICE DENSITÉ DANS L'ÉTAT EXCITÉ. — 
Les approximations (IV-3), qui impliquent aussi 
pratiquement © 7, > 1 permettent d'intégrer 
(11-18), en cherchant directement la solution forcée 
de l’équation avec second membre et en négligeant 
la variation de p,,,, devant l’exponentielle 

Qilum—m") SL EE ot 
On trouve 
Emme’ | ettr—-mo —{u—u)olt 
< mJlex,.Dlu > < w’le,, . Dim’ > 
LP + i(m—m) oe — (u—p) of] 


Ce résultat, joint à (III-4), permet de calculer 
les intensités de la lumière de fluorescence réémises 
dans les différentes directions. Nous avons égale- 
ment vu ($ II) comment calculer la quantité de 
lumière absorbée. Notre formalisme permet donc 
de décrire l’évolution de l’état fondamental sous 
l'influence du pompage optique, et de calculer 
d'autre part les signaux de détection optique. 


(LV-5) 


B. Interprétation physique. — 1. SrRucTURE 
DES ÉQUATIONS. — Les équations (IV-2) montrent 
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qu’ un élément p,, de la matrice densité est couplé 1 


à tous les autres éléments p,» ur tels sur 
u—p=u eu. 
Elles se divisent donc en groupes distincts, 


21 + 1. équations couplées pour p — pu" — 0, 
21 pour u — pu’ — 1, et ainsi de suite. Elles sont, 


à cet égard, très ‘analogues aux équations de rela- | 


xation de Bloch-Ayant-Kubo-Tomita [2] décrivant 
l’évolution d’un ensemble de spins soumis à une 

interaction aléatoire avec le réseau, qui varie assez » 
vite pour que l’approximation de rétrécissement : 
par le mouvement soit valable. L’approximation 

séculaire est également utilisée dans ce cas, la. 
condition or 7, > 1 s'écrit alors wo T > 1 où T 
est le temps de relaxation. 


2. ÉVOLUTION DES POPULATIONS. — Les popu- - 
lations des niveaux y sont les p,,, qui satisfont à 
deu he) Auu 
dt OT, To À Bu! bu (IV-6) 
On pose 
pl, 
Puy — = ra I< me, Dlu >[? 
To mn 
lCu(F,m; m—u,u)f? (IV-7) 


1 "2 
7 1<ule.-Dim >|? est la probabilité par unité 


de temps de passage du niveau y’ au niveau m 
par absorption d’un photon ; 
[Car(F, m;, M-—Kt, u)f? 

est la probabilité relative de passage d’un niveau m 
déterminé à Fun des niveaux x par émission spon- 
tanée. P,_,, représente donc la probabilité par 
unité de temps d’un passage du niveau pu’ au ni- 
veau par résonance optique. D'ailleurs, on sait 


queY [Cu(F, m; m—u,u)? = 1 et l’on peut 
1m 

écrire : 

A 

+ — FE | mles, .Dlu >|? 


= = 2 IS m lex .Diu >|? |Car(F, m; m—w,w)|? 


D mu 


= D Puy (IV-8) 
[Ps 


(1V-6) s'écrit alors : 


deu 


di -È Pau) Qu +2 (Puy) Qu'u (IV-9) 
L'interprétation en est alors évidente : la varia- 
tion de la population Puy du niveau y est égale, par 
unité de temps, à tout ce qui arrive des autres 
niveaux ’, diminué de tout ce qui part vers les 
autres niveaux y’. Ceci est à rapprocher des-équa- 
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tions obtenues en théorie de la relaxation. La prin- 
cipale différence est que, dans notre cas, 


Pin a Pins 

au moins en général. Si bien que l’on peut obtenir 
des différences de populations appréciables entre 
les niveaux par suite de l’absorption et de la réémis- 
sion de photons. C’est le principe même du pom- 
page optique. La matrice des (27 + 1) quantités 

u—>uw n'est autre que la matrice de pompage 
optique déjà utilisée par divers auteurs [3]. 

Remarquons que bien entendu 


d de, 
gs lrace p — &È Dur — 0; 


c’est-à-dire que la somme des populations de l’état 
fondamental est constante. 


3. ÉVOLUTION DE LA COHÉRENCE. CAS PARTI- 
GULIER DE 2 NIVEAUX. — Dans ce cas, qui est par 
exemple celui de l’isotope 1#H9, la matrice p ne 
comporte qu'un seul élément non diagonal, p1 1 

3 


. . . 2 
(que nous ne distinguons pas de p_ 11 qui en est 
22 


complexe conjugué), puisqu'il n’y a que 2 sous- 


niveaux Zeeman (1 =7, b= + ;) L'équa- 


tion (IV-2) s’écrit alors : 


CERN NE 
TRES 
É 2 2 2 9 Th 
& nn me 
TT 22 5e 
DT 
Die 
+ re Ë IV 
Ty T'+ we — of) ii (IV-10) 
D'où 
si + ie) — to) 
PRESS (4) = p1 1 (Lo) € Cr (1V-11) 
2 2 202 
Avec 
at 3 4 1PT 2 
Non ee 2 pet T2 
# F2 TR IT lee — 0) 
(IV-12a) 
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1h 2 f | 
; 4 =" ee 22 2 De of 
e— AE (as re) Te RE L los age 
| (LV-12b) 
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a) Interprétation de T, : T, est la constante de 
temps avec laquelle le module de p1_1décroît. 


2022 
1/7, est la probabilité par unité de temps de dispa- 
ee de la cohérence. 1/T, comprend 2 termes, 
e{°r, 


est la demi-somme des inverses des durées de vie 
[T1 = A1 [Ty AIT _1= A_1_1/7T, des niveaux 


2 22 2 22 
u — + 1/2 et u — — 1/2. Il représente la dispa- 
rition de la cohérence sous l’effet de l’excitation 
optique qui arrache les atomes à la superposition 
cohérente dans laquelle ils se trouvent. On peut 
dire qu’il représente la contribution de la rela- 
xation longitudinale à la relaxation transversale. 

Le 2° terme : 


représente la fraction de la cohérence qui s’est 
conservée au cours du pompage et qui se retrouve, 
après absorption et émission d’un photon, dans 
l’état fondamental. Nous montrons dans l’appen- 
dice IT qu’il est toujours négatif, et inférieur en 
module au 1% terme. D’après l’expression (III-11) 
de Bt#, cet effet n’existe que s’il existe 2 transi- 
tions de même nature (7x, 6, ou 6), donc paral- 
lèles dans le diagramme de Grotrian partant res- 
pectivement des niveaux 1/2 et — 1/2 (fig. 1a). 
Ainsi il n’existe pas pour la composante hyperfine 
dela raie 2 537 À partant duniveau F — 1/2 du ®Hg 


en excitation o. (fig. 1b), mais il existe pour la com- 


posante hyperfine partant du niveau F— 3/2 (fig. 1c), 
La variation de. cet effet avec le champ magné- 
tique suit une loi analogue à celle de Peffet 
Hanle [4]. On doit pouvoir l’observer expérimen- 
talement en constatant une diminution de la contri- 
bution de l’intensité lumineuse à la largeur des 
raies de résonance magnétique lorsqu’on diminue le 
champ magnétique. 


b) Znterprétation de e : La matrice densité o(é) en 
représentation de Schr°dinger est ljée à p(t), ma- 
trice densité en représentation d'interaction, par 


6 (1) — e—iJCot p(t) eidCot, 
de eo DE 3 
ee 2 2 2 
F= hs ® ® (l ENS: 
4 s4 
EU PAUL [Ets 
Lea LS se 2 
2 2 
Fire, de 
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Donc 


et dans le cas examiné : 


(7 +iag+is) to) 
% 


(to) € (IV-13) 


ans 
E 
£ représente donc un changement de la différence 
d’énergie entre les sous-niveaux Zeeman. D’après 
(IV-126), e comprend aussi 2 termes. 
Le 1% terme AE’ {A:11— À _ 1 1)représente le 


22 2 
déplacement dû à la différence des selfs-énergies 
AE'Ayiet AE’A 1 _1 des 2 sous-niveaux + 1/2 *. 


22 sus 
Nous avons déjà donné plus haut l'interprétation 
physique de AE" ($ II, 4). 

Remarquons que notre calcul s’applique à une 
composante hyperfine particulière de la raie de 
résonance. Dans le cas de l’isotope Hp, il s’agit 
par exemple de la transition X,, (ou #42) reliant le 
niveau hyperfin unique F = 1/2 du niveau fonda- 
mental 614$, au niveau hyperfin F# — 1/2 (ou 
F — 3/2) du niveau excité 6%P,. Nous supposons le 
centre k, de la raie excitatrice suffisamment proche 
de Æo1 et koo Comme il s’agit d’une self-énergie, 
nous devons tenir compte également des transi- 
tions virtuelles vers l’autre niveau hyperfin 
F = 3/2 (ou F — 1/2) du niveau de résonance et 
vers les autres niveaux excités. Considérons l’effet 
des transitions virtuelles vers les niveaux hyper- 
fins ke d’un autre niveau excité. 4, — kr varie 
alors infiniment peu d’une composante hyperfine à 
l’autre. D’après (IT, 8, b) il en est donc de même 
pour AE”. Il est facile alors de voir que par suite 
des règles de somme sur les 4,,, l'effet de ces 
transitions virtuelles est de déplacer les 2 sous- 
niveaux — 1/2 de la même quantité. 

Le déplacement de la raie de résonance magné- 
tique est donc sensible uniquement aux transitions 
virtuelles vers l’ensemble des niveaux hyperfins du 
niveau excité proche de k.. 

Le 2° terme représente un déplacement d’énergie 
fondamentalement différent du précédent Il est 
lié aux transitions réelles vers le niveau de réso- 
nance. Il n’existe que si une fraction de cohérence 
se conserve au cours du pompage Lors du temps 
r = 1/[ passé en moyenne dans l’état excicé, cette 
fraction de cohérence « voit » un effet Zeeman 
et non plus «y. D’où un déphasage de l’ordre de 
(oe — 6x)/l qu'il faut multiplier par 


F2 


D? Eos = œ)? 


ec! 
NW 
Di Die 


(*) Cet effet a été récemment mis en évidence par l’un de 
nous dans le cas du pompage optique du *#?Hg [5]. 
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fraction de la cohérence qui se conserve Si le | 
champ magnétique est.grand, l’effet est petit parce | 
qu’une faible fraction de la cohérence se conserve. | 


1 NA 
No 


S'il est petit, l'effet l’est aussi parce que les effets 


Zeeman et les déphasages le sont. Il y a optimum 


lorsque we — y? TT. Sur le Hg, on peut ima- ne | 
giner des expériences où le déplacement, dû à cet... 


effet, de la raie de résonance magnétique serait de 
l'ordre de 1/2 de la largeur due à l’excitation 
optique, donc facilement observable. 


4. ÉVOLUTION DE LA COHÉRENCE DANS LE CAS 
GÉNÉRAL. — Dans ce Cas, on a : he 


d ’ 1 : 
je (pl à) 
L Burur 18 Qu" 
AUDE Ty LT + tu — p°) (we — of) 
CESR Lu” ï 
u'=u"=u—t 
(IV-14) 


A[Tour et uw ont des définitions analogues 
à (IV-12). Leur interprétation a déjà été donnée. 
Les termes nouveaux représentent un transfert 
partiel vers le couple (4, 4’) de la cohérence exis- 
tant dans le couple (4”, 4”), par l'intermédiaire de 


l’état excité. Le système auquel satisfont les ps, . 


pour u —wu” fixé, est très analogue au système 
obtenu en théorie de la relaxation. La matrice asso- 
ciée à ce système [matrice d’élargissement [6]] a 
ses valeurs propres à partie réelle négative (cf. 
appendice II). Le pompage fait donc disparaître 
la cohérence comme un processus de relaxation. 


V. ÉVOLUTION GLOBALE 
EN PRÉSENCE D’UN CHAMP 
DE RADIOFRÉQUENCE. 


1. Position du problème. — Nous avons jusqu’à 
présent étudié la variation dans le temps, sous 
l'influence d’une excitation optique, d’un ensemble 
d’atomes dans l’état fondamental, représentés par 
une matrice densité p(t), connaissant la matrice p(4,) 
non diagonale a priori (c’est-à-dire qu'il peut y 


avoir Cohérence à l’instant t,). Nous allons main-. 


tenant étudier l’évolution de l’état fondamental 


sous l’action simultanée de l’excitation optique et 
d’un champ de radiofréquence H;, perpendiculaire 
au champ statique créant l’effet Zeeman, tournant 
à la fréquence © voisme de of. Nous négligerons 
Paction de ce champ sur l’état excité, soit que ce 
soit très différent de ©, soit que la durée de vie 
Tr — 1" 1 de l’état excité soit trop courte pour qu’il 
puisse subir une action sensible de la part de H,. 
L’Hamiltonien représentant l’action de H, est : 
JORF — EE + 


g—iot L JT eiot] (V-1) 


Now? 


où y est le rapport gyromagnétique de l’état fonda- 
mental, | son moment cinétique. 

Indiquons dès maintenant que les équations fi- 
nales que nous obtenons dans ce cas sont analogues 
aux équations (IV-2), les seules différences étant : 

a) p(t), matrice densité en représentation d’inter- 
action, est remplacée par p(t), matrice densité dans 
le référentiel À tournant à la vitesse angulaire © 
autour du champ statique : 


p (t) — eilo—0;)Tat e(t) e—i@— 0) Ta (V-2) 


b) il apparaît au second membre un terme 


‘es il, pee 
où sa 
LR = (of — ©) 73 + YH, lo (V-3) 
représente l’eflet du champ statique et du 


champ H, dans le référentiel À. 

c) «x doit, dans les équations, être remplacé 
par ©. 

Pour établir ces résultats, nous étudierons succes- 
sivement le processus d'absorption (cf. $ IJ) et 
celui de réémission (cf. $ III). : 


2. Équations du problème. — Au lieu de passer 
en représentation d'interaction pour l’état fonda- 
mental et l’état excité, nous effectuerons la trans- 
formation qui consiste à passer dans le référentiel 
tournant À pour l’état fondamental, et en repré- 
sentation d'interaction pour ce qui concerne l’état 
excité et le rayonnement. On pose : 


ay — pinot UM 
Em: —k; == ciE—k;+ mo )t Am ki (V-4) 
du: Ris KA = ENUO TEEN Gui Like 


Cette transformation remplace 7 par J et les 
équations de départ (1-9) deviennent, avec cet 


- Hamiltonien supplémentaire : 


. ; S ER 7 ss 
| 14y — » Ju Au’ +- »> << LU, Ki kolJCr|m 2 Em :—k; 
Lu” i 


mt 


| 
(V-Ea) 
ibm:—k; — >ine mlrlu’, Love ay 
Lu’ 
( : re 
Ne re mlibriu”,K, À > ap'iki:Kka (VE) 
U'kA 
iau” 5—k;: k A — a Ju’ au het, KA 
u’ 
| PRIT UT ASC r|m > bmi-k;. (V-Ec) 
mm 


Les éléments de matrice de 3, sont obtenus en 
remplaçant © par © dans les expressions (1-10) 
de ceux de J7. On a posé KL, = < pK" >. 
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3. Processus d'absorption. — Comme au $ IT, 
on intègre (V-5c) et on reporte dans (V-5b) obtenant, 


ibm —k; = Die m|#rlu’, ki, do > au 
u' 


s t v 
Sr 18 di < mir) -kN> 


u°kA 
um 
<u’le-ie-tilu > < uk, XSr()m > Em: kif). (V-6) 


Le calcul s’effectue alors de la manière classique 
indiquée au $ 11. Si l’on explicite les éléments de 


matrice de 7, on voit apparaître une exponentielle 
ele Etm'e,—uo—%] EE : comme k, est très supérieur 
aux éléments de matrice de $, à m'a et à u&, 
on peut remplacer cette exponentielle par ef (fn 
ce quirevient à remplacer < u’le-it-0}y > par 
< plu > = 0, Il reste alors 


Ce résultat était prévisible : la radiofréquence agis- 
sant sur l’état fondamental ne modifie ni la durée 
de vie, ni la self-énergie de l’état excité, On a 
donc : 


| Mines | 
ibm ik; = — 1Ë iAF) bi ki 


HE < mlXrlu ki ko > au (V-1) 
me 


L’élimination de b»:_K, ; entre (V-7) et (V-5a) se 
fait alors de manière identique à ce qui à été fait 
au $ IT, puisque n'intervient pas. Elle conduit à 


du. BR 
— On — 1 41) / 
di Puu [ Pluu 
1 Sur e Eu 
(ir + AE) TE Au burn cltu—u"ot 
7 me 
1 ste M y ee. 
SES 2T, — 1ÂE D App Puu’ eitt re °- ot, (N -8) 
Dr mn é 
On trouvera de même 
d fl ; ; ; 
di Pmmt = — L'Emm + T, D CORTE nel 
D uu’ 
< mes Dlu > < u’le0DÎm’ > euu. (V- 


4. Processus de réémission. — On généralise faci- 
lement le caleul du $ (III, 1). On intègre (V-5c), 
on prend le complexe conjugué de l'expression obte- 
nue et on le multiplie membre à membre par (V-5c). 
On retrouve les mêmes simplifications que dans le 
calcul fait ci-dessus au $ TIT, en raison du fait que 
émission spontanée ne dépend pas de la radio- 
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fréquence dans l’état fondamental, On obtient 
ainsi : 


d'ou re 

Ce = = — it plu 

LH OP ZX  eiu-ute,—ot Cir(F,m;m—4,u) 
mm’ 


m—m ht 


Cur(F, m', m — 4, w) pmm’(t). (V-10) 


L’élimination de pm entre (V-9) et (V-10) donne 
alors 


dt Pyu’ Le “1 
UE — — ÿ[Rpluue 
di [ Put 
T Burger eiltu=ur—(u"—vot Blé, 
+ (V-11) 


Tone D + lu — 4) &e — (4 — pol 


5. Équations d’évolution complètes. — En ajou- 
tant (V-8) et (V-11) (en ne gardant bien entendu 
qu’une fois le terme contenant #&), on obtient 
l’équation d'évolution globale de6. On arrive aux 
équations (V-12) 


d ’ ; 5 
“ue. = —i[x,phy 
4 0 / y RES 
= [5 (Auu + Auu) + 1AE'(Auu — Au] Puu’ 
p 

3 a Hs 

M TT UPS à EL LÉ 110 un 
BR Dr placo 


en utilisant toujours l’approximation séculaire qui 
impose la condition de validité supplémentaire 
(V-13) : YA, € «w. (Les éléments de matrice de 
étant de l’ordre de y.) 

Cette condition exprime que la largeur due à la 
radiofréquence des raies de résonance magnétique 
est petite devant la fréquence de ces raies. 

Le système (V-12) montre bien que le pompage 
optique peut être globalement considéré comme 
un processus de relaxation. Il montre que les termes 
discutés en (1V-3b) interviennent bien pour dé- 
placer les raies de résonance magnétique. Il devrait 
être complété par l’introduction (plus ou moins 
phénoménologique) de termes décrivant la rela- 
xation proprement dite (désorientation des atomes 
par collision contre les parois, les molécules d’un 
gaz étranger, etc). La résolution du sys- 
tème (V-12) permet de calculer les signaux de 
détection optique. 


6. Signaux de détection optique. — a) Lumière 
absorbée par unité de temps ; elle s'obtient aisément 
en prenant la trace de (V-8) 


Ù 1 A p/\t “ : ” 
LS Cr ar ) > App’ Puy eitu—u"ot 
tr. up 


+ complexe conjugué,  (V-14) 
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On voit qu'une modulation peut apparaître aux k 11 
fréquences rw/2r (r entier) (expériences de fais- | 1 


ceau croisé [7]). 


b) Lumière de fluorescence émise par unité de … | 


temps, dans une direction et avec un état de polari- 
sation donné. On l’obtient en calculant 
Lr = 


aa 
_— duikskÀ d'u; &kA] 
a Der 


ET Z eitm—-mot <ule1 .D|m > pmm < m'lex DIR > 


(V-15) 
AVEC Pmm donné par la généralisation de (IV-5) 
Da 2 = 2 Euu’(t) eillm—m)og—(u—14)01 
< mlex Dlu > < p'lexo.D|m' > (V-16) 


D + of(m — mm) © — (up — y) 0] 


CONCLUSION 


En résumé, nous pouvons classer les résultats 
de cette étude en 3 parties : 

1. Vue de l’état fondamental, et en présence ou 
non de radiofréquence, l’excitation optique par la 
raie de résonance peut être décrite mathémati- 
quement et physiquement comme un processus de 
relaxation qui modifie la répartition des popu- 
lations et détruit la cohérence. Cette destruction 
n’est cependant pas totale au cours d’un cycle de 
pompage : Une fraction de la cohérence peut se 
conserver au cours du cycle et ceci d'autant mieux 
que l'écart Zeeman des raies optiques est petit 
devant la largeur naturelle du niveau excité. 

2. Les signaux de détection optique : absorption 
ou fluorescence, s’expriment entièrement en fonc- 
tion des éléments de la matrice densité dans l’état 
fondamental. Les expressions exactes de ces 
signaux sont établies. 

3. Des déplacements de fréquence des raies de 
résonance magnétique peuvent être provoqués par 
Pexeitation optique. Le mécanisme de ces dépla- 
cements relève de 2 processus physiques entière- 
ment différents : La self-énergie des sous-niveaux 
Zeeman en présence du rayonnement excitateur 
qui peut différer d’un sous-niveau à l’autre ; la 
conservation de la cohérence au cours du cycle de 
pompage qui permet de ramener dans l'effet 
Zeeman de l’état fondamental une partie de leffet 
Zeeman de l’état excité. 


APPENDICE I 


JUSTIFICATION 
DE L’APPROXIMATION SÉCULAIRE 
(cf. $ IV-2) 
1 — Pour simplifier, considérons d’abord l’équa- 
tion 
dxfdt = (a + beicit) x, (1) 


A 0 


N° 7 
Sa solution rigoureuse est : 


b tort 
a+ 0 "t 
€ = To E ter (2) 


où #, est une constante, 


On peut écrire : 


Sep A a rée 
æ = toe# (x) X 7 einosé, (3) 
n=0 


1 Si |b/où| € 1, on peut limiter le développement 


2 = roet (4) 

qui est la solution de 
dæfdt = ax, (5) 
obtenue en négligeant le terme oscillant dans (1). 


Ce résultat peut être retrouvé autrement. On 
résout 


dæfdt = (a + Ab eioit) x (6) 
par un développement en puissances de À et on 
fait ensuite À = 1. A l’ordre 0 : 

dæ® Jdt — ax), 
d’où 
IE T en, 
A l’ordre 1 : 
dr Jdt = ax + beiogt x(0), 
d’où 


xt) — ES eiogt eat 
1Of 


ete. On retrouve ainsi tous les termes du dévelop- 
pement (3), dont on ne garde que le 1% (ordre 0) 
Physiquement, le 2° terme du second membre 

de (1) oscille trop rapidement par rapport à son 
amplitude pour avoir des effets notables. 

- 2. — Le système (IV-1) s'écrit sous une forme 
générale 

IX n 
de. 
les coefficients À, B, C, ont pour ordre de gran- 
deur 1/T,. On résoudra de même par un dévelop- 
pement en puissances de à le système 


= >> (Anm + Bam eiOtt + Cum er +...) Xm. (7) 
m 


dXn Jde = D Aym Am + XD (Bam EtStt 
mn 
+ Cum eiogt + ...)Xm (8) 
Le système 
mn 


AA dt — 3} Anm A9 (9) 
m 


obtenu à l’ordre 0 n’est autre que celui de l«appro- 
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ximation séculaire ». Sa solution générale est de la 
forme : 


x = ps K, Xrn ePit (10) | 
Tr 


où les «, et les 8, sont des constantes dépendant 
des Aum,les K, des constantes déterminées par les 
conditions initiales. Les 8, ont pour ordre de gran- 
deur 1/7,. On peut toujours supposer les œ, 
d'ordre de grandeur unité par une normalisation 
convenable. A l’ordre (1) on obtient : 


d x 
dé 


= D Anm Xh) + D (Bnm eos 
m m 


Con edo sr LE ar ERA) 


On en déduit pour les X® des solutions de la 
forme : 


sa SE > Qnr e,t pE MM, el(B,—8,)+kio,lt dt’ (12) 
Tel 0 

où les ME ont pour ordre de grandeur XJT, 

(Æ : ordre de grandeur de X,), k étant un nombre 

entier. On voit que le rapport des ordres de gran- 


(1) (0) SA 
deur de X® et XX est RE MEET 


Il en sera de même pour les rapports des termes 
d'ordressupérieur, AXES XO IX FIST 
Si © T, > 1, on pourra se limiter aux ordres les 
plus bas, et en pratique on fera une erreur d’ordre 
de grandeur 1 /o$ 7, en ne gardant que les termes 
d'ordre 0 (approximation séculaire). 


APPENDICE II 


Nous voulons démontrer que les valeurs propres 
de la matrice de relaxation Am définie dans le 
précédent appendice et associée au système sécu- 
laire ([V-2), ont leur partie réelle négative ou nulle. 

D’après la remarque faite dans (IV, B, 1) cette 
matrice se subdivise en plusieurs sous-matrices 


correspondant aux différentes valeurs de uw — 1”. 


En particulier pour & — un" — 0, nous obtenons la 
matrice de pompage optique : 


deéuu [At = — (2 Put) Fuu + > Purœu qu'u. (1) 


Montrons tout d’abord que les valeurs propres 
de la matrice de pompage optique ont leur partie 
réelle négative ou nulle. 

Nous devons démontrer qu'aucune des quan- 
tités p,(t) ne tend en module vers l'infini quand # 
tend vers l’infini. 

A l'instant initial] t — 1, toutes les quan- 
tités pyu(to) sont positives ou nulles (populations). 
D'autre part, à partir de (1), on établit facilement 
que : 


; 2 PESTE Puy’) Puu — 0. 


DA Ë Eur 


(TA 
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Si nous démontrons que les p,,(t) restent tou- 
jours positifs ou nuls au cours du temps, comme 
leur somme est constante, nous aurons démontré 
qu'ils ne peuvent tendre vers l’infini. 

Désignons par #, le premier instant après #, où 
l’une des quantités p,, (supposées toutes initia- 
lement positives) s’annule : p,u(£,) = 0. Il s’ensuit 
d’après (1) que 


d | S 
_ — Dot ’u’(t1). 
F Fun tt si w—>u) Eu’u’(t1) 


Comme les p,,(11) sont positifs par hypothèse 
ainsi que les P,.,.,, d’après ([V-7) on en déduit que 


n cui) 
— 0. 
[5 Fuu A > 


Puy ne pourra donc jamais devenir strictement 
négatif. Cette démonstration se généralise faci- 
lement au cas où plusieurs 0, s’annulent en même 
temps pour la première fois. 

Pour une valeur bien déterminée et non nulle 
de u — pu’, la matrice 4,, est la matrice d’élargis- 
sement. 

D’après la définition même des éléments d’une 
matrice densité, la solution rigoureuse des équa- 
tions d'évolution (I, 9) satisfait à 


leuu’(t)| < Ve) uw. (2) 


Mais nous avons démontré que moyennant les 
conditions (IV-3) et (IV-4) la solution rigoureuse 
de (1-9) était très voisine de celle des équations 
séculaires (IV-2), l’erreur relative étant de l’ordre 
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ne) 
de 1/or T,. L’inégalité (2) est donc valable égale: 1 
ment pour la solution de ([V-2). Comme nous (fe 
venons de démontrer que les p,,(t) ne tendent pas 15 
vers l’infini, il s’ensuit qu’il en est de même pour, | 
les |p,,| et que les valeurs propres des matrices 
d’élargissement ont leur partie réelle négative ou. 
nulle. 


Démontrons enfin que dans le cas de 2 niveaux fn 


($ IV, B, 3), la quantité 1/T; est positive. Physi- … 
quement ceci signifie que dans le cas de 2 niveaux, 
il y a une «restitution » partielle de la cohérence au. 
cours de la retombée et non une destruction supplé-" 
mentaire 


Il suffit de démontrer que B## est positif. 


< mle>, Dlu > < u'le, D|m’ 2 


CuifF,m;m—u,u) Cu(F, m'; m— u’, us 


Mais d’après l'expression (III-5) des <yle,DIm>, 
il vient en utilisant le fait que dans la sommation - | 
surmetm,m—uw—=m —4 


His 


[CulF,m;m—u,u)]? 


[CirlF, m ; HE u”, HA äâ 
[Su = 0jex Dim = m—u >1-0|? > 0° 


Manuscrit reçu le 20 février 1961. 
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LA DILATATION THERMIQUE DU MILIEU CRISTALLIN 


Par J, LAVAL, 


Collège de France, Paris, 


Les facteurs de la dilatation 


Résumé. — Les atomes, en oscillant, exercent les uns sur les autres des forces répulsives qui ont 
pour résultantes les tensions thermiques. Quand la température s’élève, les tensions thermiques 


s’amplifient et dilatent le milieu cristallin. 


Cette première étude a pour objet les tensions thermiques qui ont lieu dans les cristaux cubiques 
formés par les métaux. On suppose que chaque oscillation fondamentale de l’agitation thermique 
conserve la même énergie moyenne qu’un oscillateur harmonique, et l’on évalue les tensions ther- 
miques en fonction de la température dans le milieu cristallin libre de toute contrainte. 


Abstract. — By oscillating the atoms exert upon one another repulsive forces which have 
thermal strains as resultants. When the temperature rises, the thermal strains increase and dilate 


the crystalline medium. 


The subject of this first studyis the thermal strains that take place in the cubic crystals of 
metals. It is supposed that each fundamental oscillation of the thermal agitation keeps the same 
average energy as that of an harmonic oscillator ; and thermal strains are estimated with reference 
to temperature in a crystalline medium free from any constraint, 


A. Introduction. —Je me référerai dans cette 
étude aux deux mémoires intitulés « L’agitation 
thermique des atomes dans le milieu cristallin » [1] 
et « Les tensions thermiques dans le milieu cris- 
tallin » [2], et je représenterai les grandeurs par les 
symboles dont j'ai fait usage dans les mêmes 
mémoires. 

Selon la théorie classique de l’état cristallin 
(Debye, Max Born), les forces de rappel développées 
par l’agitation thermique des atomes obéiraient 
à la loi de Hooke, et l'agitation thermique qui a lieu 
dans un cristal serait analysable en oscillations 
harmoniques indépendantes les unes des autres. 

Cette conception est en désaccord avec les faits. 
Les oscillations thermiques du milieu cristallin ne 
développent pas des forces de rappel proportion- 
nelles aux élongations, même aux températures 
proches du zéro absolu. Néanmoins, l’agitation 
thermique d’un cristal peut être résolue en compo- 
santes harmoniques. Mais ces composantes ne sont 
pas indépendantes. Elles comprennent des oscil- 
lations fondamentales, qui se confondraient avec 
les oscillations harmoniques considérées par la 
théorie classique si les forces de rappel devenaient 
conformes à la loi de Hooke, et des oscillations 
secondaires produites par les interférences entre les 
oscillations fondamentales [1]. D'autre part, les 
atomes, en oscillant, se repoussent mutuellement. 
Les tensions thermiques sont les résultantes de ces 
forces répulsives. Elles sont fonctions des périodes 


du milieu cristallin, des dérivées troisièmes, qua- 
trièmes, cinquièmes.. de son énergie potentielle par. 
rapport aux translations qui transportent les ato- 
mes les uns sur les autres, enfin des propriétés des 
oscillations thermiques : leur vecteur d’onde, leur 
fréquence, leur amplitude, la direction de leur élon- 
gation [2]. Si la structure atomique du milieu cris- 
tallin est simple et sa symétrie élevée, il est possible 
de déterminer, par la mécanique quantique, l’éner- 
gie potentielle d’un. cristal à la température du zéro 
absolu, c’est-à-dire : lorsque l’agitation thermique 
est réduite au minimum [3], [4], L'énergie poten- 
tielle varie avec la dilatation du milieu cristallin, 
donc avec la température de ce milieu libre de 
toute contrainte. Mais, si l’on connaît l’énergie 
potentielle au zéro absolu, et la dilatation entre le 
zéro absolu et la température @, on peut calculer 
l’énergie potentielle et les coefficients d’élasticité 
NAS [5] du cristal porté à la température O (1). 
Une élévation AO de la température accroît les 
tensions thermiques, et par conséquent la dilata- 
tion. Les augmentations des tensions AT, et celles 


(:) Je donne, dans le mémoire « L’agitation thermique 
des atomes dans le milieu cristallin» [1] les expressions des 

our vie cs mpq 
dérivées secondes CYE, troisièmes Di RICE 
potentielle en fonction de la dilatation et des dérivées secon- 
des, troisièmes, ... de la même énergie, mais qui se rap- 
portent au milieu cristallin en état inerte, soit sensibilement 
au Zéro absolu de la température. 


de l'énergie 
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des dilatations élémentaires, Anges, («, B, y, d = 1, 
2, 3) sont liées par la relation : 


ATœY -2 Vars Arps. 

J'ai établi des formules qui donnent les accrois- 
sement des tensions ATY[2]. D’où l’on déduit les 
dilatations Amgs, puis l’énergie potentielle du cris- 
tal à la température ® + AO. Mais les coefficients 
d’élasticité et les tensions thermiques varient 
quand le milieu cristallin se dilate. Le calcul de la 
dilatation ne peut être fait que point par point. 
Une élévation minime AO de la température pro- 
duit une dilatation minime, qui peut être calculée 
avec une précision d'autant plus grande que l’élé- 
vation de la température est plus petite, en prenant 
pour données les accroissements des tensions ther- 
miques et les coefficients d’élasticité qui se rap- 
portent à l’énergie potentielle et à l’agitation ther- 
mique du cristal à la température initiale @. La 
dilatation connue, on calcule l’énergie potentielle 
et l’on détermine le régime de l’agitation thermique 
à la température finale @ + AO. Puis. on répète 
la même opération qui donne cette énergie et défi- 
nit ce régime après une nouvelle élévation AO!’ de 
la température. Aïnsi, de proche en proche, on 
décèle l’évolution de l’énergie potentielle et celle de 
l’agitation thermique entre le zéro absolu de la 
température et le point de fusion. 


B. Les facteurs de la dilatation. — Je consa- 
crerai cette étude à la dilatation des cristaux cubi- 
ques dont la maille élémentaire renferme un seul 
atome. Tous ces cristaux sont formés par des 
métaux. Les uns ont un réseau de cubes centrés 
(Li, Na, K, Fe,, Mo, W, ...) les autres un réseau 
de cubes à faces centrées (AI, Ni, Cu, Ag, Pt, 
Au, ...). Toutes leurs oscillations fondamentales 
sont du genre acoustique. D’autre part, en raison 
de leur haute symétrie, ils sont le siège d’une ten- 
sion thermique qui est isotrope ; et leur dilatation 
a la même propriété. Cela allège le calcul de la 
dilatation. Néanmoins, la méthode que je suivrai 
pour définir la dilatation de ces cristaux cubiques 
s’applique à la dilatation de tous les cristaux ; 
mais les calculs deviennent de plus en plus lourds 
quand la symétrie du milieu cristallin s’abaisse et 
quand les atomes qui constituent le motif cristal- 
lin se multiplient. 

Je supposerai le milieu cristallin à température 
uniforme, libre de toute contrainte exercée par le 
milieu extérieur, et, pour fixer les idées, j’envisa- 
gerai un cristal formé de N atomes, qui est done 
le siège de 3N oscillation fondamentales. J’en- 
tends par oscillation fondamentale l’ensemble des 
oscillations progressives ayant la même fréquence, 
pilotées par le même vecteur d’onde, faites par 
tous les atomes du cristal. Ainsi définie, une oscil- 
lation fondamentale sera notée Sr, l’indice S rap- 
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pelant son vecteur d’onde fondamental S, l’in- 
dice + sa fréquence v4.. Les vecteurs d’onde fonda- 


mentaux, au nombre de W, sont inscriptibles dans | 


une première zone de L. Brillouin, l’origine au . 
centre de la zone, les extrémités formant un réseau 
de densité uniforme. Chacun d’eux pilote trois 
oscillations fondamentales, qui sont rectilignes et 
orthogonales, l’une sensiblement longitudinale 
(+ — 1), les deux autres sensiblement transver- | 
sales (r = 2, 3). 

Je désignerai la température absolue par O, la 
constante de Boltzmann par k, le volume de la 
maille élémentaire par v, la masse d’un atome par | 
u, le coefficient de dilatation linéaire par « et le 
coefficient de dilatation cubique par Ê 


Lu 
ARS | 


Souvent, je prendrai pour variable 


ño = 
t=;S (© —2Ty). 

L'énergie potentielle envisagée sera toujours 
celle du cristal. Développée en série de Taylor par 
rapport aux déplacements u”, u?, ... accomplis 
par les atomes à partir de leurs positions moyennes, 
cette énergie a pour expression : 


1 
WE) — [ue FH Cr EN 


= 
4 
=: mp af 


L 
+ sp. da DÉS ur GEUS 
©: mpa afy 


frs û 


où VW, est, à la température O, l’énergie potentielle 
du cristal inerte, les atomes au repos dans leurs 
positions moyennes. Et les dérivées de cette éner- 
gie, C#_$, D ..., dont il sera fait état, seront 
celles qui se rapportent aux coordonnées des 
atomes en positions moyennes : 


Lavileanvan 
da—=mMm—}?, 
Din = D 
ee [: 2 2 VW 
da drf day | (x — m ; x? — p ; xa — q), 
UE (2) 


Les composantes harmoniques de l’agitation 
thermique ne sont pas indépendantes les unes des 
autres [1]. Les oscillations fondamentales sont 
liées entre elles par les dérivées troisièmes et sui- 
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vantes de l'énergie potentielle ; et les oscillations 
secondaires résultent des interférences entre les 
oscillations fondamentales. L'énergie libre de l’agi- 
tation thermique n’est donc pas celle d’une assem- 
blée d’oscillations harmoniques indépendantes. 
Néanmoins, je supposerai que chaque oscillation 
fondamentale conserve sensiblement, à toute tem- 
pérature, la même énergie moyenne, W, qu’un 
oscillateur harmonique de fréquence égale (1) : 


W = ño re + : 
(1e) 1 2 (3) 


Cette hypothèse rend compte, sensiblement, des 
chaleurs spécifiques mesurées. Aux températures 
basses et moyennes, c’est-à-dire tant que la tempé- 


MARS , 2 ù 
rature est inférieure à Op, si le cristal fond à la 


température Or, l’énergie totale des oscillations 
secondaires reste minime par rapport à l’énergie 
totale des oscillations fondamentales ; et les cha- 
leurs spécifiques mesurées se confondent avec celles 
qui sont prévues par la théorie classique (Debye, 
Max Born), laquelle suppose l’agitation thermique 
. décomposable en oscillations harmoniques indé- 
pendantes et affecte à chacune d’elles l'énergie 
moyenne d’un oscillateur harmonique. Aux tempé- 
ratures élevées, proches du point de fusion, l’éner- 
gie moyenne des oscillations secondaires entre en 
compte. Elle est déterminée par les énergies 
moyennes des oscillations fondamentales. Et, si 
Von attribue à chaque 6scillation fondamentale 
l’énergie kO, les chaleurs spécifiques calculées, 
compte tenu de la part prise à l’énergie thermique 
par les électrons de conductibilité et, s’il y a lieu, 
par le magnétisme, se confondent, à peu près, avec 
celles qui sont mesurées sur les cristaux cubiques 
envisagés, portés à haute température. 

Outre l’énergie moyenne, W(3), des oscillations 
fondamentales, le calcul de la dilatation fait état 
des facteurs suivants : 

1) La fonction d’Einstein (?) : 


1 /0W x? et 
ol ee () 


(:) Je fais cette hypothèse pour les raisons suivantes, 
Jusqu'ici, seule a été quantifiée, en toute rigueur, énergie 
d’un oscillateur harmonique. L’agitation thermique du 
milieu cristallin est un mouvement complexe. Son énergie 
est définie par une somme de termes nombreux. Et la quan- 
tification de cette énergie ne peut être abordée que par le 
calcul des perturbations, poussé jusqu’à des ordres élevés. 
Ce calcul est fort long, son résultat incertain. De sorte que 
l'énergie libre de l’agitation thermique, et partant l'énergie 
moyenne d’une oscillation fondamentale, échappe présen- 
tement à tout calcul précis. 

(?) Selon le symbolisme en usage, l’indice w signifie que 


la dérivée (5) se rapporte à une pulsation & constante. 
o 
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Quand la température se rapproche du zéro 
absolu, la fonction d’Einstein s’annule. Lorsque 
la température s’élève, elle croît constamment, et 
tend vers l’unité dès que x tombe au-dessous de 
l’unité. 

2) La fonction : 


p_ wo 
Her li AE 

led llErrat) ee 

BV Tee AS (5) 


Égale à 4/2 au zéro absolu de la température, la 
fonction F augmente aussi constamment quand la 
température s’élève et, lorsque x décroit au-dessous 


de lunité, elle tend one 


ho 
3) La dérivée 
DÉS RE 
(ER FE io Ÿ 
Lies enfer 2) er 26r2] 
6 (et — 1)° 3 (6) 


Quand la température tend vers le zéro absolu, 
la fonction y s’annule comme x e*. Elle passe par 
un maximum, égal à 2,69, pour x —4,5, puis 
quand x décroît à partir de 4, 5, elle finit par tendre 
vers zéro comme 2%,0 (fig. 1). 


() 4,5 10 15 


HTC ME 


4) La dérivée (négative) 


(as) 
do \w?/@ 


Cette dérivée s’annule quand la température 
tombe au zéro absolu et, lorsque la température 
s'élève, elle diminue constamment et tend vers 


2 


wo? 


aÿ ç® (e 


+1) 
wo? (e% — {hs S- (7) 
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o) Le coefficient 
CnemeEsele y = 3 TD 


= — 5 DC (ui)e + 208 at afl;a = mp, (8) 


cl et cl? sont les coefficients d’élasticité et 
le coefficient de compressibilité, isothermes, tous 
rapportés aux constantes de rappel atomiques 
2, CE (contrevariantes, égales aux dérivées 
secondes de l’énergie potentielle (1-2)). 

6) « Les facteurs de structures » gs., des oscil- 
lations fondamentales. Soient : 2/ la hauteur de la 
maille cubique, et !1, l2, 18, 3 vecteurs dirigés sui- 
vant ses arêtes et de même module, {. Les transla- 
tions du réseau cristallin, celles qui transportent un 
atome sur un autre, peuvent être exprimées : 


m= ml + ml? + mal. 


Les nombres entiers, M1, Mo, Ma, sont assujettis 
aux conditions suivantes : 

19 Si le réseau est cubique centré, M1, Me, Ma, 
sont tous les trois pairs, ou tous les trois impairs. 

20 Si le réseau est cubique à faces centrées, m1, 
M, Ma, Sont tous les trois pairs, ou l’un d’eux est pair 
et les deux autres impairs (1). 


X4 


pf1o0] 


O 


Pro? 


Les axes de coordonnées étant dirigés suivant les 
axes de symétrie quaternaire, envisageons le cristal 
divisé en deux parties À et B par un plan P (100) ; 
et calculons la tension thermique principale T4 [2], 

> 
de direction et de sens Ox, (fig. 2), (donc normale 
aux plans 100), exercée par le demi-cristal B sur 
le demi-cristal A. Pour ce faire, il convient de som- 
mer, selon la méthode que j’ai définie [2], les forces 
répulsives exercées par les strates 100 véritables, 
d’épaisseur {, contenues dans la partie B, sur les 
strates 100, de même épaisseur, contenues dans la 
partie À. On obtient, par ce procédé, les facteurs gsr 
(positifs) qui se rapportent à la maille élémen- 


(1) Bien entendu, pour définir toutes les translations du 
réseau, il est nécessaire de prendre, au lieu d’un ou deux 
nombres pairs, un ou deux Zéros. 
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taire (2). Et que le réseau soit cubique centré, ou | 
cubique à faces centrées, ces facteurs sont donnés 
par la même expression : 


gr = — À cpsr Cy.sr D | D'ÉYsin?, x$a 
EY (2 
1 
2 J >) D18,* J ab AE > Die Jab | a} Ps (9) 
b=a jol> la] 
B, y = 1,2,3, 
où 
Xav = sin?.7r$Sa + sin?.7rSb — sin?.7rS(b — a) (10) 
GO ER EE PE TA 
et 
a—m—p; b—m—q; 


p et q étant, comme m, des translations du 
réseau cristallin. 
Enfin je ferai état des relations : 


2 DA D CO 

de LE EU — PU, 

dw °duw dv Feet 
—= à = (2 , 

co dv Ge dv 

dg  dg dv dg 
Et = 40] ; 

de dv de ” dv 

ac de dv ae 

— = —.— —= AU —. 

de dv de dv 


Une élévation de la température accroit les ten- 
sions thermiques, donc dilate le milieu cristallin, 
et par là modifie sa structure électronique ({). La 


(2) J'ai défini les facteurs gs- (notés £G,-) quise rappor- 
tent à la maille cubique d’un cirstal dont le réseau est formé 
de cubes centrés ([2], pp. 584, 585, 586) : 


BST = D CR.Sr Cy,ST (EÏ,gy.5r SE SE,Gv.S7 “0 VE ey.sc). 
BY 

J’ai rapporté ces facteurs à la maille cubique afin de don- 
ner une méthode générale pour évaluer les tensions qui ont 
lieu dans le milieu cristallin. Mais l'expression obtenue ne se 
prête pas à des calculs simples. De plus, pour les cristaux 
envisagés ici, les facteurs gs+, rapportés aux mailles cubi- 
ques, prennent deux formes différentes, l’une propre au 
réseau de cubes centrés, l’autre au réseau de cubes à faces 
centrées. 

(5) Les constantes de rappel, DIBY, DIPY, qui figurent 
dans la formule (9) sont égales aux dérivées troisièmes de 
l’énergie potentielle (1-2) : 


Di8r — Diby, Difr = DIX (p + @). 


Elles sont posées contrevariantes, et dans un système 
orthonormé, elles ont les mêmes modules que les constantes 
de rappel (covariantes) : 

D D jk ñ ? 
16Y 167 
qui leur correspondent dans les formules 6 et 7 [2] (j, 
k — 1,2) mais sont de signes opposés. 

(*) En particulier, les électrons qui se trouvent au som- 
met de la bande de Fermi sont portés par l’agitation ther- 
mique à des niveaux d'énergie plus élevés que les niveaux 
où ils étaient primitivement. 
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. dilatation et la modification de la structure élec- 
. tronique réduisent l’énergie potentielle du cristal 
en grandeur absolue ; partant, elles diminuent le 
coeflicient © (8), et affaiblissent dans l’ensemble 
les pulsations &,4. et les facteurs g4. Les dérivés 
do do dg d€ 


totales =, ; de: dei 


effets. 


tiennent compte de tous ces 


C. Les tensions thermiques. — Il est commode 
d'exprimer chaque tension thermique par une 
. série [2]. Et dans les cristaux cubiques envisagés, 
la tension thermique, 7, isotrope, définie posi- 
tive (gs. > 0) a, de la sorte, pour expression : 


DÉS DE (12) 


où le terme 7, (1)est de degré 27 par rapport aux 


amplitudes quadratiques moyennes 4,4, des oscil- 
- lations fondamentales [2]. 
- La tension principale T, [2] peut être formulée 


Ti; (13) 
uU 
il ia 
== SRE 14 
Q N 5 Fe $7 Le 
SLPÉS NE 19,3, 
et 
[4] Ws- £s- 
2157 de 


Afin de simplifier les formules, j’écrirai de même, 
entre deux crochets, affectés de l’indice Sr, des 
expressions où l’énergie W, la pulsation ©, le fac- 
teur g, et les dérivées de ces grandeurs, se rappor- 
tent à la même oscillation fondamentale S7. La 


T: [ : ] T. 
5 . 5 


compte SN termes. Donc, sauf si le cristal est 
extrêmement petit, à une température donnée, le 
facteur Q (14) est constant, caractéristique de 
l'espèce cristalline. 
Du zéro absolu à la température ©, la dilatation 
cubique des cristaux considérés s’élève à 


® 
(6) = f Eco) 4@. 


Donc, si la demi-arête de la maille cubique est 
l, au zéro absolu, elle devient, à la température © 


(16) 


ë(@) 
lo) te 3 À 
eb 
ô(@) 
a@) = de 3 
Cie Mo nee 210472 nr, (17) 


Quand la température s’élève, l’énergie poten- 


{) La tension étant isotrope, nous avons : 
HU TNT = Tr. 
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tielle du cristal Wio, (1) et ses dérivées (2) décrois- 
sent en grandeur absolue ; donc les facteurs gs (9) 
et les pulsations w$, diminuent aussi. Mais, au 
cours de toute la dilatation les dérivées 


d Log.v 


d Log. CC 
d Log.v 


CAD 


conservent à peu près la même grandeur (?). À 
Si les dérivées de l’énergie potentielle sont C&, 

D, ..., au zéro absolu, à la température O, elles 

ont sensiblement pour expression : 

—(g28 + n4800))00), 

Ro et id Pt 

—(paÊT +név8e;)0(e) 

e 


serolene"e. 6 nus ele eu etre tele eme on ele ele + 6 eee sister on else tie nr eneheie 


pET, mir, étant des coefficients 


constants, et 
[n°130 «€ [pa°, [na°"|D(e) < lgev|, le 


La dilatation thermique des cristaux pris en 
compte est isotrope. Donc à température crois- 
sante, les vecteurs d’onde S gardent tous les mêmes 
directions dans le réseau cristallin, et leur modules 
restent dans des rapports constants. De ce fait, les 
facteurs e**"Sa, sin?.7Sa, ab (10) sont invariables. 
En outre, les « vecteurs propres » c$. de compo- 
santes Ce.s- (x —1,2,3) ne changent pas considéra- 
blement de direction dans le réseau cristallin (#). 
A ce changement près, négligeable, les facteurs 
£s- (9) sont seulement fonctions des dérivées troi- 
sièmces DT, DT, et de l’arête 21 de la maille 
cubique (par l'intermédiaire des translations at) (17). 


EST = — À ChgrCy, Sr 


1 N 
it 1 DÉDE 


> Dire sin?.rSa 
a 
l eue 
Lo outre +mea)a 
US »] Dire e D se 4 ab 
+ ab 
bi =: 
[l NI N 
— (ar; ngre)s 
EP SRDI À xab |Ga, (19) 


ab 
RAEATA 


(*) Ce fait sera établi plus loin (formules (64) et (66)). 

(*) Les coefficients œ%6, 128, paêr, naPT, sont, à 
quelques rares exceplions près, positifs. 

{*) Les « vecteurs propres » cs des oscillations fonda- 
mentales se confondent sensiblement avec les vecteurs 
propres des matrices de Fourier (20). Et, si à température 
croissante, les dérivées secondes Ca (1-2) diminuaient uni- 
formément (toutes dans la même proportion) les vecteurs 
propres des matrices de Fourier seraient constants ; donc 
les vecteurs propres cs; ne varieraient pas sensiblement 
Quand la température s'élève, la décroissance des dérivées 
secondes Cab n’est pas uniforme ; mais les cristaux cubi- 
ques formés par les métaux ayant une symétrie élevée, les 
vecteurs propres €s- gardent néanmoins à toute tempéra- 
ture, à peu près la même direction dans le réseau cristallin. 
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D'autre part, les oscillations fondamentales sont 
entretenues par des forces de rappel qui restent, 
même si la température est élevée, à peu près 
linéaire par rapport aux dérivées secondes C% 
(a —m—p). De la sorte, les carrés w%, de leurs 
pulsations se confondent sensiblement avec les 
valeurs propres de la matrice de Fourier [y] 
ayant pour éléments : 


vb = 15 CS eienSa — ee DH cr sin?.7Sa, (20) 
a 4 
donc 
2 \ —(pu8 +-ne8S)S 
o°?sr — TPEADS Ca, St Cp,sr D (ous e (R : Ta ) 
EL af a 


sin2.nsa. (21) 


Lorsque les distances |a|, [b| augmentent, les 
dérivées C%, D, DT, diminuent vite en gran- 
deur absolue. Partant, les pulsations w$. (21) et les 
facteurs gs. (19) sont presque entièrement déter- 
minés par les dérivées secondes C%, et par les 
dérivées troisièmes D£%*, DST, qui portent sur des 
atomes voisins. Quand la température s'élève, ces 
dérivées ne décroissent pas unifermément. Mais 
du zéro absolu au point de fusion, leurs décroissan- 
ces relatives 

[Cl ICter . 1DÉ —ID8]. | 
Poe 


tombent rarement au-dessous de 1/2 et restent 
donc, dans l’ensemble, du même ordre de gran- 
deur. En outre, chaque carré de pulsation, «%., 
est fonction de toutes les dérivées secondes C%, 
chaque facteur g$ de toutes les dérivées troi- 
sièmes DST, D (1). Donc si les pulsations 5. et 
les facteurs gs. prennent les grandeurs oo et 
£sro au Zéro absolu, et à la température les 
grandeurs ose) et £suo) quand S et + varient, à 
toute température, par un effet de moyenne, les 
deux affaiblissements 


ERONE Ge (6), (22) 
EST 
2 
HO 5 (0), (23) 
OST 


n’oscillent que faiblement. 
Et, à température croissante, l’évolution de tout 

facteur 

£8r'0) 


dre w8-(e) 


Ste qsrl9) 
87 rsrl®) ; 


AA) 


(24) 


(*) Car à toute dérivée DIPY correspondent deux dérivées 
Def, DE Ye égales : 
DPY = DEPV — DAY 
telles que les translations a’, a”, se déduisent de la transla- 
tion a par des rotations de et de . autour d’un axe de 


symétrie ternaire. Les dérivées Dig ont la même propriété. 
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comme celle du facteur moyen 


il gsr(®) ; 


ue 0 nee ST 
210) 3N $r sr 9) 


est déterminée par le rapport entre la diminution 


d'ensemble, en grandeur absolue, des produits … 


D'Y. at, DÉ.a1, et la diminution d’ensemble, de 
même en grandeur absolue, des dérivées secondes 


%, la première diminution étant définie par les. 


affaiblissements se) (22), la seconde par les 
affaiblissements rsue) (23). 
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(5). 


à 


Pour une même élévation de la température, la 


décroissance relative, en grandeur absolue, des … 
dérivées troisièmes D#°T, DéÉY, excède dans l’ensem- 


ble celle des dérivées secondes C#(2). Mais l’aug- 


mentation des translations at (17) rétablit approxi- | 


mativement l’égalité entre la diminution relative 
des produits DT. at, DT. at, et celle des dérivées 
ce 

En conséquence, du zéro absolu de la tempéra- 
ture au point de fusion, le facteur moyen ëçe, (25) 


conserve à peu près la même grandeur, caractéris- … | 


tique de l’espèce cristalline ; et la variation rela- 
tive 
lEsrol — Ésrto) | 
lÉSro| 


de tout facteur &,. reste légère (de quelques cen- 
tièmes) (#). On a donc grossièrement : 


d Log.«sr 14 Log.gsr 
d bog.v - 2 d Log.v 


(?) Les énergies potentielles mutuelles des électrons de 
conductibilité, celle des ions positifs, celles des électrons 


de conductibilité et des ions positifs, sont en gros conformes 


à la loi de Coulomb. Même les énergies d'échange entre les 
électrons obéissent à cette loi ; et les énergies de Van der 
Waals résultent d'interactions coulombiennes. Ainsi, l’éner- 
gie potentielle d’un cristal relève principalement de la loi 
de Coulomb. En conséquence, lorsque le milieu cristallin 
se dilate, les dérivées de son énergie potentielle (1-2) 
décroissent, en grandeur absolue, d'autant plus rapidement 
que leur ordre est plus élevé. 
(8) Si les égalités 
PUB — pe, V8 = M, 
et GREY = QÊT = ps TAT = TOY = 


étaient satisfaites, on aurait : 


o 
Éce — ÉÊte e— Lea (113)- pa+(ns—1n2)D]8 
oO 
(2 . 2 To 
et les égalités 
1 
P3 RE Tu Oo Ms 02 
entraînerait la relation : 
2 ST 
Er — 5" — constante. 
o 
STo 


C’est la décroissance non uniforme des dérivées Dib: 
DIfY et C£ qui fait que les différents facteurs £s- ne varient 


pas de même avec la dilatation 5s-, donc avec la tempé- 
rature ©, mais leurs variations restent minimes. 


(26)» S 


M x 7 


Soit V» (£) la fréquence moyenne des oscilla- 
tions fondamentales, qui ont lieu dans le cristal à 


LA F . . . 
la température moyenne. Si ces oscillations ont 


des fréquences distribuées selon une densité p(v) 


entre les deux limites w et vr ( < vr), 


et L 
Or 1 ŸL 
va (5) = ge etui va (27) 
J’évaluerai les tensions secondaires To, Ta, ... 
(formule 12) au moyen des expressions : 
b 
nee" ie 
T, = * ie (28) 
Pi W: (29) 
ÉTÉ uv m> 


où b et c sont des coefficients constants, positifs, 
propres à l’espèce cristalline, et où l’énergie 
moyenne W, est celle d’un oscillateur harmonique 


4 (] La 
dont la fréquence, v» (2) , est constante (indépen- 


dante de la température (3-4)). 

Ces expressions des tensions secondaires sont 
approximatives ; et, parmi les approximations 
faites, deux d’entre elles sont particulièrement 
importantes. 

1° Aux basses températures, chaque oscillation 
fondamentale S+ prend en moyenne dans le temps, 
une énergie Ws. (3) qui lui est propre ; car cette 
- énergie est fonction de sa pulsation w,., laquelle 
diminue lorsque le milieu cristallin se dilate. A ces 
énergies individuelles W4, est substituée dans les 
expressions (28) et (29) une énergie moyenne, W», 
fonction uniquement de la température puisqu'elle 


 serapporte à une fréquence constante v» (5) (27). 


Donc si la température du cristal est basse 

(o<7 
4) 

erreur sur la grandeur des tensions secondaires. 
Mais aux basses températures, les tensions secon- 
daires restent minimes par rapport à la tension 
principale 7, et une lourde erreur sur leur gran- 
deur n’affecte que légèrement la tension globale 
D lle dénrs at 0) 

(12) et celle des vées set oz (): 
20 Si la température du cristal est élevée 


le cette substitution entraîne une forte 


(1) La tension globale Test fonction de la température 
et de la pression. Mais le cristal étant supposé sous une pres- 
i Il t az 
sion nulle, 45 et 463: 
- uniquement de la température. 


sont des dérivées totales, fonction 
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(0) > 30), selon notre hypothèse, toutes les 


oscillations fondamentales ont sensiblement la 
même énergie moyenne, #O, qui est indépendante 
de leur fréquence. Donc remplacer les énergies 
individuelles W,. par l'énergie moyenne W,, n’en- 
traine pas non plus une erreur notable sur la gran- 
deur de la tension globale et sur celle des dérivées 
GE os EU 
do’ d@? 

La tension 7, est fonction linéaire des facteurs 


/ 


o # _) 
sr, sr" — es (30) 
OST 0) ,S7" 
la tension 7, des facteurs 
0 Ver ren 
Es, s'r,8"r" — RARE (31) 


sr W°$/r D°$/r" 


Et les coefficients b et c sont respectivement les 
grandeurs moyennes, prises sur toutes les oseil- 
lations fondamentales, des facteurs Es.s et 
Es, sr, 87 

Les facteurs gs se et Less. Sont comme les 
facteurs gs. (9) proportionnels à l’arête 27 de la 
maille cubique. Et, les oscillations secondaires 
négligées, chaque facteur gs. varie linéaire- 
ment avec les dérivées cinquièmes de lénergie 
potentielle W4(®) (1), chaque facteur gse,s,s". 
varie linéairement avec les dérivées septièmes de 
la même énergie. D’autre part, les produits &$. 
Œ et ©$ © ©ès sont sensiblement des 
fonctions quadratiques et cubiques des dérivées 
secondes C% (2-18) (1). En majeure partie l’énergie 
potentielle du milieu cristallin obéit à la loi de 
Coulomb. Par conséquent, pour une même éléva- 
tion de la température, les dérivées cinquièmes de 
l’énergie potentielle décroissent dans l’ensemble, 
en grandeur absolue, à peu près dans la même pro- 
portion que les carrés des dérivées secondes ; et, 
dans l’ensemble aussi, en grandeur absolue, les 
dérivées septièmes diminuent un peu moins que les 
cubes des dérivées secondes. Donc, à température 
croissante, l’arête de la maille cubique se dilatant, 
les facteurs Es... et surtout les facteurs Es. 5,81 
augmentent dans l’ensemble, et partant les coefi- 
cients b et c grandissent. Mais du zéro absolu de la 
température au point de fusion, la dilatation cubi- 
que Ô(O) (16) ne s’élève qu’à quelques centièmes ; 
il est donc à présumer que dans le même intervalle 
de température les coefficients b et c ne croissent 
que dans une proportion minime par rapport à leurs 
grandeurs au zéro absolu. 

Même quand le cristal est à une température 
élevée, à peine inférieure à la température de 


(1) Je suppose les oscillations fondamentales entretenues 
par des forces de rappel linéaires par rapport aux dérivées 
secondes et les vecteurs propres cs: (de composantes 
CSI = 2 NS) CODSIANLS: 
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fusion, les tensions secondaires n’entrent dans la 
tension globale que pour une faible proportion (de 
l’ordre de 10—1). Donc, même si les facteurs Es..5, 
Esrsx,s augmentent considérablement quand le 
milieu cristallin se dilate, cela n’entraine qu’un 
accroissement minime de la tension globale. 

Mais le coefficient de dilatation linéaire « est, 


2 A pre L daT . 
en gros, proportionnel à la dérivée 10’ et sa varla- 


d 
tion en fonction de la température est principale- 
ment déterminée par la dérivée —— 10% 

La comparaison des coefficients de dilatation 
définis par cette théorie et des coefficients d’élas- 
ticité mesurés, révèle d’abord que les tensions 
secondaires 7, T3, ..., s’exercent dans le même 
sens que la tension principale T;, puis qu'aux tem- 
pératures élevées proches du point de fusion, la 
somme positive 


ds de 
40 410: 


reste plus petite mais est du même.ordre de gran- 
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deur que la dérivée = ee et que la somme 


10 J 
d° Ts, d°Te 
d®? d0? ? 


positive aussi, devient, en valeur absolue, beau- 
Fes laquelle 
décroîit et finit par être négative quand la tempé-. 
rature est voisine du point de fusion. 

En conséquence, si les facteurs Es.s+ (30) 
Es. srs"r" (31) augmentent notablement quand la 
température s’élève, cette augmentation entraine 
aux températures élevées, un accroissement sen-. 
sible du coefficient de dilatation. 

Même aux températures élevées, la tension TE 
reste minime par rapport à la tension T;, et la” 


coup plus grande que la dérivée 


tension T; reste minime par rapport à la tension | 


T;. La série (12) qui exprime la tension globale” 
converge rapidement. Et pour calculer le coeffi-" 
cient de dilatation je ne retiendrai que les tensions. 

secondaires T, et T3. 


Manuscrit reçu le 11 avril 1961. 
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REMARQUES 
SUR L'UTILISATION DES POUDRES 
DANS LES ÉMULSIONS NUCLÉAIRES 


Par Cao xuan Cauan et NGuyËN Bic Nuu, 
Office de l'Énergie Atomique, Saïgon 


Une des techniques récentes dans l’étude des réac- 
tions nucléaires par les émulsions nucléaires consiste 
à incorporer des poudres de la matière à étudier entre 
deux couches d’émulsion et à observer les étoiles de 
- désintégration issues de ces grains [1, 2]. 
Mais cette méthode introduit deux nouvelles dif- 
ficultés : 
- 1. On ne voit pas l’origine des désintégrations. En 
_effet les particules chargées issues de la désintégration 
ne laissent des traces observables qu’une fois sorties des 
grains. Si ces particules n’ont pas assez d'énergie pour 
s'échapper du grain, le phénomène est perdu. Cet effet 
entraîne donc une imprécision dans l'interprétation 
. des résultats. Il serait intéressant de pouvoir faire une 
estimation théorique de l’importance de ce phénomène 
- et de faire les corrections nécessaires dans la statistique 
finale. 
2. Durant son parcours à l’intérieur du grain, la 
particule peut subir des chocs contre les noyaux du 
. milieu, Ces chocs provoquent des déviations des tra- 
jectoires et entraînent des perturbations dans l’étude 
L des distributions angulaires. 
Évidemment les corrections sont d’autant plus 
faibles que le diamètre des grains est plus petit (cor- 
- rections sur les parcours aussi bien que sur les angles). 
Mais on est limité pour réduire le diamètre des grains, 
… lorsque les sections efficaces sont faibles. Si le diamètre 
_ moyen est trop petit, très peu de désintégrations se 
produiront dans les grains. 
Dans l’étude présente, nous nous proposons de faire 
une analyse théorique du premier de ces deux effets. 


fee 


HiGeate 


Pour simplifier le problème nous supposerons que 
les grains sont des sphères de rayon À. 
Soit une particule &« émise avec une énergie E à 
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l’intérieur du grain. Pour que cette particule soit 
observable, nous devons avoir : 


FAP ARS (1) 


AE étant la perte d'énergie de la particule dans la 
sphère. À 
En étant l’énergie minima pour que la trace de la 
. particule « soit observable dans l’émulsion, 
Nous prendrons pour 


Em 


les ions °Be Em 


— 0,2 MeV environ 
0,4 MeV environ. 


les œ 
D’après la relation de Bethe [3], on peut écrire : 


knet z? NZ 
= . Ln 
mv ? 


2mv°? 


KZ @) 


(E) 


Nous supposerons négligeable la perte d’énergie 
dans le grain de sorte que f(Æ£) ne varie pratiquement 
pas pendant ce parcours 
ze : charge de la particule. 


: numéro atomique des atomes constituants du grain. 
: masse électronique. 

: vitesse de la particule. 

: nombre d’atomes par cm° du milieu traversé. 
L’inégalité (1) peut s’écrire : 


E — f(E) At > Em 


ER TION 


ou 


FE) ê 
B est défini comme le plus grand parcours permis pour 
la particule dans la sphère, c’est-à-dire que si ce par- 
cours est plus grand que B la particule ne peut pas 
laisser une trace observable. 

Supposons que la particule « soit émise au point M. 
O est le centre de la sphère. Posons 


ONE 
La particule & suit la trajectoire MI pour sortir 
du grain 
ES 
MISE IOM—=0; 
Nous avons 


(Ax)? = R2 + 7? — 2Rr cos 0. (4) 


De (3) et (4) nous tirons : 
R? Er —9Rr cos 0 < B2 
ou 


RE rÈ BE 


2Rr de 


cos 0 > 
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Une discussion simple montre que la probabilité 
pour que cos Ô soit supérieure à X est 


D 1 si K <—1 
PS He ANS (4) 


FE étant ee sur le cercle trigonométrique défini 
par cos ® — À (fig. 2). 


Tee 


D'autre part, la probabilité p’ pour que le point M 
soit à l’intérieur d’une couronne sphérique de rayons 
r et r + dr est : 


krr? dr 3r? dr 
dp — — 


Ur 


_ La probabilité pour que la particule « soit observée 
est le produit pp’. En faisant l'intégration sur toute 
la sphère nous aurons : 


DE. FR 
2 
= DA, mi R? + r2 = pp? 
ë 2Rr 


d’où 


CIRE at 3 [ er 1 (5) 
8R? 16 8K?| J(E) 16° 


‘On distingue deux valeurs de F: Ejnet En pour 
lesquelles P = 0 et P —1. 

Les corrections doivent être faites lorsque 
Enin < E < Epax. Si Ê < Enim la particule ne peut 
pas laisser de trace observable dans l’émulsion tandis 
qu’on peut toujours l’observer si £ > Fax 

Deux facteurs entrent en jeu dans notre problème : 

1) L'énergie de la particule «x, Æ au moment de 
l'émission. 2) Le rayon À? du grain incorporé dans 
J’émulsion. En général l’énergie Æ est imposée par les 
conditions expérimentales. Dans ce cas, la relation (5) 
permet de choisir le rayon À? de manière à éviter ces 
corrections et de se mettre ainsi dans les conditions 
d'observation les plus favorables. 

Comme exemple nous avons construit des courbes 
montrant les variations de la probabilité P avec l’éner- 
gie Æ pour différents rayons /? pour les particules « 
et les ions °Be dans le cas des grains de diamant [1]. 
Pour construire ces courbes, nous supposerons que la 
constante k de la formule (2) est k — 11,2 eV. En outre 
l'énergie minima £i pour que la trace soit observable 
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est supposée égale à.0,2 MeB pour la particule « et. 
égale à 0,4 MeV pour l'ion ‘Be (fig. 3). 


Particule 


IPS SUIS 


F1G. 3. — Variation de P en fonction de l’énergie. 
P probabilité d'observation de la particule. 
E énergie de la particule au moment de l’émission. 
R-rayon des grains. 


Cette méthode peut être bien entendu utilisée pour | 
d’autres particules chargées et dans le cas où plusieurs - : 
atomes de nombre atomique différent sont présents 
dans le grain. 


Qu'il nous soit permis d’exprimer ici notre vive 
reconnaissance à M. R. CHasrez, Professeur à la Fa- 
culté des Sciences de Bordeaux, pour les précieux 
conseils dont nous avons bénéficiés au cours de ce 
travail. 


Lettre reçue le 26 mai 1961. 
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STRUCTURE HYPERFINE DU PRASÉODYME 
DANS LES SOLIDES 


Par Yves Ayanr et Élie BELORIZKY, 


Laboratoire d’Électrostatique et de Physique du Métal 
Laboratoire de Physique Générale 


Le praséodyme est l’un des métaux de terres rares 
possédant un isotope unique doué de spin, donc fa- 
vorable à l'étude expérimentale de structure hyperfine. 
de noyau Pr possède les caractéristiques suivantes : 

= 5/2, u = 3,8 MBN ; Q — —0,054.10% cm°). 

Dans cet article, nous allons successivement étudier 
les cas d'ions Prit et Pr+ dans un site cubique pur. 


.)}, Communication au 22 col- = à 


Nas 


1° Cas de l'ion Prtt, — La valence 4 du Pr est 
connue ; il semble qu'on puisse réaliser lion Prit 
comme inclusion dans la thorine. L’ion Prit est un 
système (4j)! dont le fondamental est ?F,,,, lequel se 
décompose dans un champ cubique en 2 niveaux appar- 
tenant aux représentations l, et I ; c’est le niveau 
quadruple F, qui est le plus bas. (Lorsque l’environ- 
ee est un cube de charges négatives, ce qui est le 
cas. 

En tenant compte du spin nucléaire, le niveau fon- 
damental appartient à la représentation [y X D, 
(27 étant la représentation usuelle du groupe des rota- 
tions, mais restreinte évidemment, ici, au groupe du 
cube). On calcule aisément la réduction : 


me X Die = a JR + 1BE — 2h _. SD . Sie 


Nous voyons ainsi que le multiplet hyperfin (du 
fondamental [°;) comprend 2 niveaux simples, 2 ni- 
veaux doubles et 6 niveaux triples. 

On traite le fondamental l', comme un niveau qua- 
druple d’un atome libre, auquel on associe un moment 
cinétique fictif J', tel que J’ — 3/2. On définit les 
états [Mr > en fonction des [M; > (où My est le 
rai nombre quantique magnétique associé à J) par : 


ST SIGNE 59" Aie} TT 73h :> 
| 142>=f2> HP tr 10 > 
a —=(5j6) 2} SR > + je) [3 > 


Et on vérifie que les composantes J, du moment 
cinétique J ont des éléments de matrice entre < M'r:| 
et My > données par la relation : 


Jr = © 


Pour obtenir les niveaux hyperfins issus du niveau 
fondamental électronique L,, il suffit de diagonaliser 
le couplage hyperfin À — A I..J. La constante À a 
été estimée à 0,043 cm1 d’après une étude sur le fon- 
damental 47,,, de Pr 2%), Le plus simple est de travail- 
ler dans la base [Mr > [My > où Mr est le nombre 
quantique magnétique nucléaire, en remarquant que 
seuls existent les éléments de matrice obéissant à la 
règle de sélection. 


A(M7r + My’) = 0, + 4. 


En réduisant la matrice primitive d’ordre 24, 
on trouve une sous matrice correspondant à 
M; + My: = 0, + 4, d'ordre 6 mais symétrique par 
rapport à son centre, ce qui ramène à des problèmes 
d'ordre 3 ; de même on trouve une sous matrice Ccor- 
respondant à My + My — + 2 ayant les mêmes 
propriétés. On trouve pour 12 £.J les valeurs propres 
suivantes : 


D 33 TU 60,1, 56,1: 
I, 64,87 ; 24,68 ; — 60,55 
DORA Te 47 15; 


Ci-contre le diagramme des niveaux (fig. 1). 
L’écartement total du multiplet est de 0,46 cmt, 


Remarques : 1) Le terme quadrupolaire, très faible, 
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(ordre de grandeur 1074 cm1) n’introduit aucune 
décomposition supplémentaire ; il a par conséquent 
été négligé dans cette étude. 

2) Effet Zeeman : Au 1e ordre, il est nul pour les 
niveaux [',, [,, l, et isotrope pour les niveaux l, 
et l'; ; cela est important car il serait possible de faire 
des expériences de résonance paramagnétique sur une 
poudre, à condition d'utiliser un champ magnétique 
assez faible, 


20 Cas de l’ion Pr°+. — Le niveau fondamental 
de l'ion Prè+ libre (4f)? est °/7,, qui en potentiel cu- 
bique pur se décompose en [, + P, + 1, + F5; les 
valeurs de ces niveaux sont proportionnelles à 14 
(P,) ; 7 (l,) ; 203) — 13(l;) dans l'hypothèse vrai- 
semblable où le potentiel cristallin du 6° ordre est 
négligeable — le signe correspondant toujours à 
8 charges — au sommet d’un cube). 

Le fondamental de l’ion est ainsi un niveau (élec- 
tronique) triple l',, auquel on associe, comme plus haut, 
un moment cinétique fictif J’ avec J' = 1. Il convient 
de poser : 


FA = MBPS > MR) 
| 0 = f1/2}/2] 2 > — (1/2) 222 
EM = Nes cr sen 


avec les mêmes conventions que plus haut. Dans la 
base [My > on trouve aisément qie J ro (1). 


Au premier ordre du caleul des perturbations, nous 
avons donc, pour le couplage hyperfin magnétique, 
à traiter la matrice À I1..J. En remplaçant J selon (1), 
on obtient un couplage H — A'I.J' (2) (A' a été es- 
timé à 0,092 cm1 de la même façon que A(1) ; (2) se 


70 


r, 
60 
50 T: 
40 Ts 
30 
20 Le 
10 De 
0 
-10 
- 20 
-30 T, 
-40 i 
_50 Te 
- 60 T 
-70 r' 


TGS 
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résout aisément en introduisant F' — I. + J' et on 
obtient 3 niveaux F' = 7/2, 5/2, 3/2. 

En fait ces niveaux sont respectivement associés aux 
représentations. 


Lo Dojo = Te + Lo + Ve: Vo Dia 
DE NS A 1 TV; Dyyo = Ts. 


Donc un calcul plus poussé faisant intervenir le 
2ème ordre du couplage magnétique et le {er ordre du 
couplage quadrupolaire, fera apparaître une décompo- 
sition supplémentaire. À priori ces deux effets se tradui- 
sent-par un hamiltonien effectif Æ' ajouté à Æ et qui 
sera de la forme : 


“H=alJ +6 (Id) + ce +de JEe 
k 


Seul le terme en d joue un rôle : a, b, c, ne donnent 
lieu à aucune décomposition des niveaux /", mais seule- 
ment à un très léger déplacement. 


Tous calculs faits, la contribution à d du terme ma- 
12 


5] ème . 27 A \ 4 
gnétique au 2ème ordre est : 80 AZ °Ù AE est l’écarte- 


ment des deux plus bas niveaux cristallins. 

La contribution du terme quadrupolaire, est de 

35 7 
— 107$ em 1, Le calcul de la décomposition se 
présente ainsi ; on passe en base |F, Mr > et il 
suffit de former les sous matrices 
< F',M'r|H|F", Mr > 

Si F' — 7/2 et 5/2, ces sous matrices sont nécessaire- 
ment du type bien connu «/(F") + 8 


3 
2F' +1 


HONTE LE 4 Fe Terre 

Les matrices Æ7(F") sont les mêmes que celles qu’on 
rencontre pour représenter les couplages cristallins 
du 4ème ordre en symétrie cubique ; il est aisé d'obtenir 
leurs valeurs propres. En calculant directement des 
éléments de matrice simples, on trouve : 


Pour F' — 7/2 | 
ML eus 
ÉTAT Ta 
Pour F7=15/2 
“ A 34 
LÉ EU E 


et on utilise le fait que les matrices Z(F") ont respecti- 
vement les valeurs propres : 


HC/? k2(T$) ; 6(T) ; 
H6/? 6(Ts) ; —42(T;) 


Cela fournit la décomposition des niveaux Æ° = 7/2 
et 5/2. Le niveau F” — 3/2 ne subit aucune décompo- 
sition. 

(La décomposition supplémentaire n’est pas à 
l'échelle sur la figure 2.) 


— 54(T,;) 


FTrGa2: 


Conelusion. — Nous étudions actuellement le Pr 
en présence de champ cristallin non parfaitement 
cubique, le terme supplémentaire pouvant être tétra- 
gonal, trigonal, ou même comme dans le cas des gal- 
lates, totalement dépourvu de symétrie. 


Lettre reçue le 26 mai 1961: 
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COMMUNICATION A LA SOCIÉTÉ FRANÇAISE DE PHYSIQUE 


L'ÉMISSION PHOTOÉLECTRIQUE : 
EFFET DE SURFACE, EFFET DE VOLUME 


Par P. VERNIER, 
Faculté des Sciences de Dijon 


Malgré de nombreux travaux, l’origine des élec- 
trons émis par les solides n’est pas encore éclaircie. 
Nous définirons d’abord le problème dans un schéma 
corpusCulaire : 

Dans un solide, un électron subit l’action d’un pho- 
ton qui le fait passer d’un état Æ à un état !, Dans le 
-premier état, l’électron ne peut traverser la surface du 
solide, dans le second il a une probabilité p, de le faire. 
Gette interaction doit conserver la quantité de mou- 
vement et l’énergie entre le photon absorbé, l’électron 
excité et le solide. L’action du solide est nécessaire, 
car l’énergie échangée entre un électron libre et un 
photon lumière dans l’effet Compton est négligeable, 


… On peut distinguer trois types de transitions : 


1) Effet de surface : La conservation de la quantité 
de mouvement est assurée grâce à l’action sur l’élec- 
tron du champ électrique superficiel qui enferme les 
électrons non excités dans le solide. 

2) Effet de volume : La conservation de la quantité 
. de mouvement est assurée grâce à l’action du champ 
périodique qui règne dans le solide. 

3) Effet de défaut : On peut encore envisager une 
action des champs perturbateurs dus aux impuretés 
et autres défauts, superficiels et profonds. Les carac- 
tères de l’émission photoélectrique qui en résulte sont 
fondamentalement différents des précédents. C’est 
pourquoi, nous le classerons à part sous le nom d’effet 


£ de défaut. 


L'effet de surface peut être seul envisagé dans l’hy- 
pothèse où les électrons sont libres dans un métal 
parfait. À partir du modèle de Sommerfeld on peut 
ainsi calculer la répartition spectrale du rendement 
quantique et la répartition énergétique des électrons 
émis. L'accord est relativement satisfaisant près du 
seuil, mais, d’après Mayer [2] une théorie fondée sur 
un effet de volume pourrait donner des résultats très 
voisins. 


Le point de vue ondulatoire, — Pour tenir compte 
de la polarisation de la lumière, il convient de décrire 
celle-ci par une onde. Une théorie mixte, ondulatoire 


| pour la lumière et corpusculaire pour les électrons peut 


donner des résultats mais nous adopterons ici une 
description purement ondulatoire. 

L'état électronique d’un volume V de solide peut 
être décrit par un ensemble fini de fonctions d’ondes, 
orthogonales entre elles et normées sur le volume V : 


.Ekt 
ok!r) = uk{r) exp (- 27] h } 


À l’intérieur d’un cristal parfait, og(r) peut être mise 
sous la forme : 


kr) = uk(r) exp 277 kr (2) 
où ug(r) est une fonction possédant la triple périodicité 
du cristal et k un vecteur pris dans la première zone 


de Brillouin. Dans le voisinage de la surface d’un cristal 
parfait ok(r) peut être mis sous la forme : 


kr) = Uk‘) Xp 277 kr (3) 


k; représente la composante du vecteur k parallèle à 
la surface et v4,(r) est une fonction possédant la double 


périodicité de la face de sortie du cristal. 

Dans le voisinage d’une impureté ug et 4, perdent 
tout caractère périodique. 

L’onde 44 ne s’annule pas aux limites du volume V. 
En fait, pour les états électroniques normalement oc- 
cupés dans le solide, elle présente vers l’extérieur du 
solide une forme évanescente à laquelle aucun courant 
n’est associé. 

L’influence d’un champ électromagnétique variable 
peut être considérée comme une faible perturbation 
apportée à l'équation du Schrôdinger dont les 44 sont 
les solutions non perturbées. Pour représenter l’état 
électronique du solide, on peut, dans une première 
approximation, remplacer chaque 44 par la somme : 

dk + CkIKt) GI (&) 
où dy représente une fonction d'onde, solution de l’équa- 
tion de Schrôdinger non perturbée, mais correspondant 
à un état électronique d’énergie Æ, trop élevé pour 
être occupé normalement. Si on tenait compte de 
toutes les actions subies par les électrons, notamment 
de celles qui sont responsables de l’effet Joule dans les 
conducteurs qui ferment le circuit photoélectrique, 
c(t) serait une constante, En composant les densités 
de courant associées à tous les termes cy(t), pour tous 
les vecteurs k et 1, à l'extérieur du solide, on obtien- 
drait un courant photoélectrique stationnaire. 

Il est plus commode de ne faire intervenir que le 
champ électromagnétique de la lumière. Si on cherche 
la contribution d’un volume V à l’effet photoélec- 
trique, quand lé potentiel vecteur peut, dans le volume 
V, être mis sous la forme : 


À —= (ar) cos 2vt (5) 
on trouve : 


€ y XP 2j (VkI — v)t 
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avec : 
Hu — | sa graû qu 
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représente la probabilité pour qu’un électron dans 
l’état 4 passe dans l’état (y pendant l’unité de temps. 

Si la densité de courant j, associée à dy n’est pas 
nulle à l’extérieur du solide, on peut associer, à une 
transition vers l’état dy, une probabilité d'émission 
pi-pi peut être calculé en comparant le flux de ÿ; à 
travers une surface entourant la photocathode et à 
travers une section bien choisie du volume V. Le cou- 
rant photoélectrique Z peut être obtenu en faisant 
la double somme : 


T1 =eZ rk pl. (8) 
kl 


Ce courant est la seule quantité observable, avec la 
répartition géométrique et énergétique des électrons 
émis. Il n’y a donc rien d'étonnant à ce que xx pour 
un couple kl donné, varie avec le temps et dépende des 
dimensions du volume V d’une façon analogue. 


Les règles de sélection. — cyy prend des valeurs 
beaucoup plus importantes sik et L satisfont certaines 
relations, qu’on peut appeler règles de sélection. En 
raison du volume fini du cristal et du temps d’obser- 
vation fini, les transitions à considérer dans l’évalua- 
tionde Z sont celles qui satisfont approximativement 
ces règles, qui correspondent aux lois de conservation 
de l’énergie et de la quantité de mouvement dans le 
schéma corpusculaire. Les règles de sélection diffèrent 
suivant l’effet envisagé : 


A) EFFET DE SURFACE. — On envisage les contri- 
butions, à l'intégrale Hki, du gradient moyen de 
qui ne diffère de zéro qu’au voisinage immédiat de la 
surface. En utilisant l’expression(2) on trouve : 


Hi = d vj,*(r) @ grad vk,(r) exp 27j(k; — l;) rdr. (9) 
Y 


La règle de sélection pour l’effet de surface est donc : 


ou 


La (10) 


Si on peut négliger les variations périodiques de 
vx((r), (cas des électrons libres), grad vy,(r) ne pré- 
sente pas de composante parallèle à la surface du 
solide. Il résulte alors de 6, 7, 8, que le courant photo- 
électrique est proportionnel à da’, a, représentant la 
composante normale de l’amplitude à du potentiel 
vecteur, 


B) EFFET DE VOLUME. — On envisage la contri- 


bution à {4 des variations périodiques de uk, d’après 
l'expression 2, 


Hi À u*(r) a grad wy/r) exp 2j(k— ljrdr (11) 
v 
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Cet effet est nul si l’approximation des électrons … À 
Jibres est valable (ux — Cte dans le solide). - 
La règle de sélection est ici plus sévère que dans | 
l'effet de surface : 


VE = V TOR 


KE 


ar €, Sel eos 


mais l'intégrale Æ4y a des contributions provenant | de. 
tout le volume du cristal et plus seulement de la zone | 
superficielle. 

Un effet vectoriel est possible, parce que la proba- 
bilité d'émission p, de l’électron excité dépend de law 
direction du vecteur Let que la direction de grad wx(r) 
dépend de la direction de k= 1. Mais la COMPOSARES 
normale a, n’est plus forcément la seule active, ni. 
même la plus active. 4 


C) EFFET DE DÉFAUT. — Les défauts internes et" | 
superficiels provoquent des perturbations dans 
uk(r) et vx(r) et détruisent leur périodicité. Il en ré-" 
sulte une contribution à A4. Les règles de sélection 
et les effets correspondants sont différents de ceux 
trouvés dans les deux premiers effets. | 


Conclusion, — La polarisation de la lumière est un. Z 
test important sur l’origine des électrons dans l’émis- 
sion photoélectrique. Les résultats expérimentaux déjà 
obtenus pour certains métaux indiquent un aceroisse-" 
ment spectaculaire de l’émission, quand l’angle de la. 
vibration lumineuse avec la surface augmente. Ceci 
pourrait faire penser à une contribution dominante, 
de l’effet de surface. Toutefois l’argument n’est pas 
absolument décisif, puisque cet effet vectoriel n’est 
pas exclusif de l'effet de surface. En fait, d’autres. 
éléments expérimentaux font penser à un effet de, 1 
volume [2]. | 

Dans les semiconducteurs, on admet généralement … 
que l'effet de volume domine et il serait intéressant den 
développer les études expérimentales sur l'effet 
vectoriel, : 


Section Bourgogne-Franche-Comté. 
Séance du 28 avril 1961 à Besançon. 
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